main.m

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

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