c语言编写常微分方程例题,c语言编写求解常微分方程初值问题

本文档利用 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); }

展开阅读全文

  • 0
    点赞
  • 0
    评论
  • 0
    收藏
  • 一键三连
    一键三连
  • 扫一扫,分享海报

表情包
插入表情
评论将由博主筛选后显示,对所有人可见 | 还能输入1000个字符
©️2021 CSDN 皮肤主题: 数字20 设计师:CSDN官方博客 返回首页
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值