求解ODE（常微分方程）的C程序（数值解法）.docVIP

求解ODE(常微分方程)的C程序(数值解法) 将(Ieuler Rkutta gear IMRK63 Adambm Hamming VSRK4)放入TC2.0 的Include文件夹中后,执行TEXT程序即可 TEXT Ieuler Rkutta gear IMRK63 Adambm Hamming VSRK4 /*主程序TEXT*/ double Y_initial=3.0; /*初始条件y(0)，作为外部变量*/ int N=11; /*节点数N，也作为外部变量*/ #include stdio.h #include math.h #include stdlib.h/*以上3个均为系统提供的库函数头文件*/ #include IEULER.C #include RKUTTA.C #include GEAR.C #include IMRK63.C #include ADAMBM.C #include HAMMING.C #include VSRK4.C/*以上7个均为程序开发人员（我）自定义的函数，即7种算法*/ double f( double x, double y ) /*Dy/Dx＝f（x，y）中的f（x，y）*/ { double f; f=x; return(f); } void g( double x,