四阶龙格-库塔法求解常微分方程的初值问题

算法原理和程序框图

龙格—库塔法是一种求其准确解y(x)在一系列点x_i处y（x_i）的近似值y_i的方法，y_i称为数值解。经典的四阶龙格库塔法方程如下：

y'=f(t,y),y(t₀)=y₀输出按如下求解y_n+1=y_n+h(k₁+2k₂+3k₃+4k₄)/6其中                       k₁=f(t_n,y_n)

k₂=f(t_n+h/2,y_n+hk₁/2)

k3=f(t_n+h/2,y_n+hk₂/2)

k4=f(t_n+h,y_n+hk₃)

4.2程序使用说明

本程序单独编写了一阶初值问题单独求解和高阶常微分方程的求解问题。常微分方程的阶数，x的上下限和步长直接输入就好。常微分方程则在系统提示input后再输入，并且按照y(1,1)=y,y（2,1）=y’,y(3,1)=y’’的格式输入,然后输入初始条件f(x),f’(x),f’’(x)等计算结果。最后求得所要的函数值。

4.