二阶龙格-库塔(共5页).docx

精选优质文档-倾情为你奉上#include#include /*n表示几等分，n+1表示他输出的个数*/ int RungeKutta(double y0,double a,double b,int n,double *x,double *y,double (*function)(double,double) double h=(b-a)/n,k1,k2;int i; x0=a; y0=y0; for(i=0;in;i+) xi+1=xi+h; k1=function(xi,yi); k2=function(xi+h/2,yi+h*k1/2);yi+1=yi+h*k2; return 1; double function(double x,double y) return y-2*x/y; int main() int i;