main.m

用四阶古典RK方法解初值问题： 取h=1/8。每隔8步打印出数值解与真解的值（u(t)=(t^2)/2-t）

clear;
[T,U]=kutta(0,3,0.125,0);
plot(T,orig(T),'b.-',T,U,'m*-');
legend('精确解','h=1/8时的数值解');
xlabel('自变量');
ylabel('解');