2017-11-27から1日間の記事一覧
メモ。導出はTaylor展開してtにh,2hを代入して解きます。 4次精度の微分はこう。 (f(x - 2h) - 8 f(x - h) + 8 f(x + h) - f(x + 2h))/(12h) + O(h⁴) function V = deriv4(Etime,Epos) em2 = Epos(1:end-4); em1= Epos(2:end-3); %Epos(3:end-2) ep1 = Epos(…
メモ。導出はTaylor展開してtにh,2hを代入して解きます。 4次精度の微分はこう。 (f(x - 2h) - 8 f(x - h) + 8 f(x + h) - f(x + 2h))/(12h) + O(h⁴) function V = deriv4(Etime,Epos) em2 = Epos(1:end-4); em1= Epos(2:end-3); %Epos(3:end-2) ep1 = Epos(…