Question

Matlab用四阶龙格库塔法求解不可微分方程组的初值

u'=x+u+v u(0)=0 0<x<1
v'=-x+u v(0)=1
w'=u+w w(0)=1

Answer 1

fun=@(x,uvw)([x+uvw(1)+uvw(2);-x+uvw(1);uvw(1)+uvw(3)]);

[x,uvw]=ode45(fun,[0,1],[0 1 1]);
plot(x,uvw)
legend('u','v','w')
%下面是相图
figure
plot3(uvw(:,1),uvw(:,2),uvw(:,3))

Answer 2

你刚刚问过同样类型的一道题，做法相同的。