求助:有关解高阶微分方程
哪位高手解答下:
初学MATLAB,要用四阶龙格库塔函数解一个二阶微分方程
(注:不调用ode45函数)
程序运行总是出错
请高手指点修改下:
方程是:y''+4y'+6y=1
y(0)=0,y'(0)=1,0<=x<=3,
程序代码如下:
syms x y z h d n m ff f1 f2
f1=sym('dy=z')
f2=sym('dz=-4*z-6*y+1')
x=0:0.1:2;h=0.1;
y(0)=0;z(0)=1;
for n=1:length(x)-1
k1=-4*z(n)-6*y(n)+1;
k2=-4*[z(n)+h*k1/2]-6*[y(n)+h*z(n)/2]+1;
k3=-4*[z(n)+h*k2/2]-6*[y(n)+h*z(n)/2+h^2*k1/4]+1;
k4=-4*[z(n)+h*k3]-6*[y(n)+h*z(n)+h^2*k2/2]+1;
z(n+1)=z(n)+h*(k1+2*k2+2*k3+k4)/6;
y(n+1)=y(n)+h*z(n)+h^2*(k1+k2+k3)/6;
end
错误提示:
??? Subscript indices must either be real positive integers or logicals.
Error in ==> Untitled34 at 6
y(0)=0;z(0)=1;
不知道这里如何修改,请朋友们帮忙,谢谢拉!!!!