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

>>阅读原文<<