龙格库塔的文字介绍 龙格-库塔（Runge-Kutta）法 到目前为止，我们已经学习了多步法，例如：亚当斯-巴什福思（Adams-Bashorth）法，亚当斯-莫尔顿（Adams-Monlton）法，都是常微分方程的积分 。
什么是欧拉方法（Euler's method）？ ？www.zhihu.com 简单来说，隐式欧拉这里 是已知的，这里的 才是未知量，F是函数。我们需要求得当F=0时，y究竟应该是多少，也就是根是多少。怎么办呢？先猜一个数，然后。
用C#编写一段代码，实现欧拉格式和龙格库塔格式。这里有一段C语言的代码，怎么改写成C#？ 直接粘过去就行了…只需要把2.0改成2或者(float)2.0就行了
四阶RK(龙格库塔)方法求初值问题的数值解 f=inline('-y+x+1'，'x'，'y')；微分方程的右边项dx=0.05；x方向步长xleft=0；区域的左边界xright=10；区域的右边界xx=xleft：dx：xright；一系列离散的点n=length(xx)；点的个数y0=1；(1)欧拉法Euler=y0；for i=2：nEuler(i)=Euler(i-1)+dx*f(xx(i-1)，Euler(i-1))；end(2)龙格库塔法RK=y0；for i=2：nk1=f(xx(i-1)，RK(i-1))；k2=f(xx(i-1)+dx/2，RK(i-1)+k1*dx/2)；k3=f(xx(i-1)+dx/2，RK(i-1)+k2*dx/2)；k4=f(xx(i-1)+dx，RK(i-1)+k3*dx)；RK(i)=RK(i-1)+dx*(k1+2*k2+2*k3+k4)/6；endEuler和Rk法结果比较plot(xx，Euler，xx，RK)hold on精确解用作图syms xrightsolve=dsolve('Dy=-y+x+1'，'y(0)=1'，'x')；求出解析解rightdata=subs(rightsolve，xx)；将xx代入解析解，得到解析解对应的数值plot(xx，rightdata，'r*')legend('Euler'，'Runge-Kutta'，'analytic')
