Runge-Kutta-c-language-program.rar_runge kutta_龙格－库塔

#include<stdio.h> #include<math.h> double fun(double x,double y) { double f; f=100*(sin(x)-y); return(f); } main() { double y[400],K[4]={0,0,0,0},x[400],h,j; int i,k; printf("输入步长 h\n"); scanf("%lf",&h); x[0]=0;y[0]=0; j=3/h; for(k=0;k<j;k++) { x[k+1]=h+x[k]; K[0]=h*fun(x[k],y[k]); K[1]=h*fun(x[k]+h/2,y[k]+K[0]/2); K[2]=h*fun(x[k]+h/2,y[k]+K[1]/2); K[3]=h*fun(x[k]+h,y[k]+K[2]); y[k+1]=y[k]+(K[0]+2*K[1]+2*K[2]+K[3])/6; } for(i=0;i<=k;i++) { printf("x[%d]=%lf,y[%d]=%lf\n",i,x[i],i,y[i]); } }