untitled2.m

利用MATLAB进行求解常微分方程组的四阶R-K算法法

clear;
%程序C8-1  二阶R-K法
M=input('输入M(M~=0)的值:');
K=input('输入K的值:');
B=input('输入B的值:');
y=[1,0]';%initial condition
b=K/M;a=B/M;
n=0;t=0;h=0.05;
y_rec(:,1)=y;t_rec(1)=0;
while t<=0.4
    n=n+1;
    K1=h*f_sm(y,t,a,b);
    K2=h*f_sm(y+K1,t+h,a,b);
    y=y+0.5*(K1+K2);
    t_rec(n+1)=n*h;
    y_rec(:,n+1)=y;
    t=t+h;
end
y=y_rec;
t=t_rec;
plot(t,y(1,:),'-',t,y(2,:),':')
y(1,:)
y(2,:)