解析解
y'' - y '+ytsint = t , t在0到1之间
y(0) = 0.4,y'(0) = 0.2 求数值解并画出解函数的图形
主函数
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
不知道你这个问题是否已经解决, 如果还没有解决的话: