数值分析欧拉法、改进的欧拉法、4阶龙格库塔法 MATLAB算法

数值解与理论解对比可知，四阶龙格-库塔法的精度已经很高，用它来解一般常微分方程已经足够了。 有程序运行说明

function y = DEEuler（f haby0varvec）
%f:一阶常微分方程的一般表达式的右端函数
%h:积分步长
%a:自变量取值下限
%b:自变量取值上限
%y0:函数初值
%varvec:常微分方程的变量组
format long;
N = （b-a）/h;
y = zeros（N+11）;
x = a:h:b;
y（1） = y0;
for i=2:N+1
    y（i） = y（i-1）+h*（-y（i-1）+x（i-1）+1）;
end
format short;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        319  2009-05-19 19:38  5\DEEuler.m

     文件        501  2009-05-19 19:54  5\DELGKT4_lungkuta.m

     文件        423  2009-05-19 19:44  5\DEModifEuler.m

     文件      61440  2010-11-04 23:07  5\Given the initial.doc

     目录          0  2010-11-04 23:07  5

----------- ---------  ---------- -----  ----

                62683                    5