C语言编写求解常微分方程初值问题.doc

下载文档

10
0
约1.75千字
约 3页
2017-05-05 发布于四川
举报
版权申诉
保障服务

C语言编写求解常微分方程初值问题.doc

1、本文档共3页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

C语言编写求解常微分方程初值问题

本文档利用C语言编写求解常微分方程的方法，包括改进欧拉方法求解常微分方程的初值问题，以及四阶龙格－库塔法求解常微分方程的初值问题。程序均采用txt文件输入输出方式。 /* 改进欧拉方法求解常微分方程的初值问题*/ #includemath.h #includestdio.h #includestdlib.h #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(Cant open this file!\n); exit(0); } fscanf(fp1,%f%f%f%f,f0,h,x,l); fclose(fp1); FILE *fp2; if((fp2=fopen(out.txt,w))==NULL) { printf(Cant open this file!\n); exit(0); } while(x=l) /*改进欧拉法算法*/ { f1=f0+h*f(x,f0); x=x+h; f2=f0+h*f(x,f1); f0=0.5*(f1+f2); fprintf(fp2,x=%f,f=%f\n,x,f0); /*输出运算结果*/ } fclose(fp2); } /* 四阶龙格－库塔法求解常微分方程的初值问题*/ #includestdio.h #includemath.h #includestdlib.h 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(Cant open this file!\n); exit(0); } fscanf(fp1,%f%f%f%f,y0,h,x0,l); fclose(fp1); FILE *fp2; if((fp2=fopen(out.txt,w))==NULL) { printf(Cant 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); }