zjode4.m

来自「基于matlab开发软件 的选择输入类型的数值计算仿真程序 迭代方法包括 」· M 代码 · 共 13 行

function [x1,y1]=zjode4(a,b,c,d,u,h,n,x0)
x=x0;
y=c*x;
for i=1:n
    k1=a*x(:,i)+b*u;
    k2=(a*(x(:,i)+0.5*h*k1)+b*u);
    k3=(a*(x(:,i)+0.5*h*k2)+b*u);
    k4=(a*(x(:,i)+h*k3)+b*u);
    x(:,i+1)=x(:,i)+h/6*(k1+2*k2+2*k3+k4);
    y(i+1)=c*x(:,i+1);
end
x1=x;
y1=y;