| 查看: 314 | 回复: 0 | ||
冉冉升起93新虫 (初入文坛)
|
[求助]
不能求出解析解的方程求数值解,新手求指导
|
|
#include<stdio.h> #include<math.h> #define step 0.01 void fm(); double fs(double); double ki; FILE *fp1; main() { fp1=fopen("fig2a.txt","w"  ;ki=0.0; fm(); fclose(fp1); fp1=fopen("fig2b.txt","w" ;ki=1000.0; fm(); fclose(fp1); fp1=fopen("fig2c.txt","w" ;ki=3000.0; fm(); fclose(fp1); fp1=fopen("fig2d.txt","w" ;ki=6000.0; fm(); fclose(fp1); } //////////////////////////////////////////////////////////////////// void fm() { float a; double t=0.0; double a0=0.01; long int i; fprintf(fp1,"t\t a\n",t,a); for(i=0;i<=150;i++)//对时间t的循环 { //a=a0+fs(a0)*step; double result; double ks=0.000001, n0=1000,i0=3000,k=200; double A,B,C,res1,res2,res3,res4,res,m,mn,a; //for(i=0;i<1000;i++)//a { A=(k/ki)*((ki+i0)/k*(1+k*n0)-i0/(1+k*ki+k*n0)); B=(1+i0/ki)/(1+k*n0); C=-i0/(ki*(1+k*ki+k*n0)); res1=pow((1+k*a),(-A/k)); res2=pow((n0-a),B); res3=pow((n0-a+ki),C); res4=exp(-ks*t-m); res=res1*res2*res3; mn=res-res4; if(mn<0.00001) { fprintf(fp1,"%f\t %f\n",t,a); } a=a+0.01; } t=t+step; a0=a; |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
两块石头
已经有8人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有14人回复
-
广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！
已经有7人回复
