f.m

来自「用四阶古典RK方法解初值问题： 取h=1/8。每隔8步打印出数值解与真解」· M 代码 · 共 2 行

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