Runge_Kutta法解初值问题

西交大计算方法上机可用。RK-4 该算法用Runge-Kutta 法解初值问题。a,b—求解区间的左右端点，y 0 —初 值，步长选用h。依次由公式计算出K1，K2,K3，K4 的值并计算出出每步的y(n)。 同时每算依次就更新依次t(n)的取值，最终输出t i 和y

%%Runge Kutta
clear;
h=1/128;   %步长 1/81/321/128
n=1/h;
%初值
y=-1;y_2=3;y_3=2;
y=[y;y_2;y_3];
i=1;
for t=0:h:1-h
   t1=t;
   y1=y;
   f1=[0 1 0 0;0 0 1 0;-1 1 1 2]*[y1;t1]+[0 0 -3]‘;
   K1=h*f1;    %计算k1
   t2=t+h/2;
   y2=y+1/2*K1;   
   f2=[0 1 0 0;0 0 1 0;-1 1 1 2]*[y2;t2]+[0 0 -3]‘;
   K2=h*f2;    %计算k2
   t3=t+h/2;
   y3=y+1/2*K2;
   f3=[0 1 0 0;0 0 1 0;-1 1 1 2]*[y3;t3]+[0 0 -