本文档利用 C 语言编写求解常微分方程的方法,包括改进欧拉方法求解常微分方 程的初值问题,以及四阶龙格-库塔法求解常微分方程的初值问题。程序均采用 txt 文 件输入输出方式。 /* 改进欧拉方法求解常微分方程的初值问题 */ #include #include #include #define f(x,y) (y-2*x/y) void main() { FILE *fp1; /*依次输入初值、步长、下限、上限*/ float x,h,l,f0,f1,f2; if((fp1=fopen(“in.txt“,“r“))==NULL) { printf(“Can't open this file!\n“); exit(0); } fscanf(fp1,“%f%f%f%f“, fclose(fp1); FILE *fp2; if((fp2=fopen(“out.txt“,“w“))==NULL) { printf(“Can't open this file!\n“); exit(0); } while(x #include #include float f(float x,float y) /*定义一个函数*/ { float f; f=(y-2*x/y); return f; } void main() { int n,N; /*定义各变量*/ float x0,y0,h,x1,y1,l; float k1,k2,k3,k4; FILE *fp1; /*依次输入初值、步长、下限、上限*/ if((fp1=fopen(“in.txt“,“r“))==NULL) { printf(“Can't open this file!\n“); exit(0); } fscanf(fp1,“%f%f%f%f“, fclose(fp1); FILE *fp2; if((fp2=fopen(“out.txt“,“w“))==NULL) { printf(“Can't open this file!\n“); exit(0); } for(n=1;x0l;n++) /*四阶龙格- 库塔法*/ { x1=x0+h; k1=f(x0,y0); k2=f(x0+h/2,y0+h/2*k1); k3=f(x0+h/2,y0+h/2*k2); k4=f(x1,y0+h*k3); y1=y0+h/6*(k1+2*k2+2*k3+k4); fprintf(fp2,“x%d=%3.1f,y%d=%f\n“,n,x1,n,y1); x0=x1; y0=y1; } fclose(fp2); }
展开阅读全文