Euler and Runge Kutta methods

向前Euler方法code function Forward_Euler(fun,u_0,a,b,N) %% 用 Forward-Euler法来求区间[a,b]的常微分数值解,N为区间[a,b]分割的份数 %u'(t)=f(t,u); %u(a)=u_0; %迭代格式 U_k+1=U_k+delta_t*f(t_k,U_k); %显式格式 % step1 将区间[a,b]等分为N份 delta

>>阅读原文<<