f.m

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

function r=f(t,u)
r=10*u+11*t-5*(t.^2)-1;