常微分方程的数值解法
odefun = @(x,y) y*tan(x) + sec(x);
[x,y] = ode45(odefun, [0,1], pi/2);
plot(x,y,'r-')
hold on
fplot(@(x)(x+pi/2)./(cos(x)),[0,1],'b--')
legend('数值解','精确解')
可见数值解和精确解几乎重合
odefun = @(x,y) y*tan(x) + sec(x);
[x,y] = ode45(odefun, [0,1], pi/2);
plot(x,y,'r-')
hold on
fplot(@(x)(x+pi/2)./(cos(x)),[0,1],'b--')
legend('数值解','精确解')
可见数值解和精确解几乎重合