Matlab用四阶龙格库塔法求解不可微分方程组的初值
问题描述:
Matlab用四阶龙格库塔法求解不可微分方程组的初值
u'=x+u+v u(0)=0 0
答
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')%下面是相图figureplot3(uvw(:,1),uvw(:,2),uvw(:,3))