TP03.rar_runge-kutta matlab

clc clear all t=0; h=0.5; tf=5; y=[0;1]; Y=[y] T=[t]; for t=0:h:tf-h; k1=h*fun2(y,t); k2=h*fun2(y+(k1/2),t+(h/2)); k3=h*fun2(y+(k2/2),t+(h/2)); k4=h*fun2(y+k3,t+h); y=y+(1/6)*(k1+2*k2+2*k3+k4); T=[T t+h]; Y=[Y y]; end Ya=(3/4).*T.^3+(3/2).*T.^2+T t1=0; T1=[t1]; y1=[0;1]; Y1=[y1]; h=0.1; for t1=0:h:tf-h; y1=y1+h*fun2(y1,t1); Y1=[Y1 y1]; T1=[T1 t1+h]; end plot(T,Y,'g',T1,Y1,'r')