matlab函数用法怎么用速度

Question

解析解
y'' - y '+ytsint = t ， t在0到1之间
y(0) = 0.4,y'(0) = 0.2 求数值解并画出解函数的图形

Answer 1

主函数
clc,clear
t = [0:0.01:1]';
x0 = [0.4 0.2]';
[t,x] = ode45('eqq',t,x0);
plot(t,x);

函数文件
function  xdot = eqq(t,x)
xdot = [x(2);
       x(2)-x(1)*t*sin(t)+t;]
end