zjode4.m

基于matlab开发软件 的选择输入类型的数值计算仿真程序 迭代方法包括 1可自启动的4阶龙格库塔法 2需要使用其他方法启动的阿达姆斯法

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;