| 查看: 1670 | 回复: 14 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
muqiuyan木虫 (小有名气)
|
[交流]
【求助】16位 MS-DOS 子系统 已有7人参与
|
||
|
求助 我在用C语言运行一个积分程序时，出现 16位 MS-DOS 子系统 提示 程序如下：#include "math.h" #include "stdio.h" double simps1(a,b,f) double a,b,(*f)(); { int n=100,k; double t,t2,s,s2,h,x; h = (b-a)/n; /* 求得初始步长*/ t = ((*f)(a)+(*f)(b))/2.0; /* 计算初始的积分值*/ for(k=1; k x = a+k*h; t = t+(*f)(x); /* 累次求和*/ } t = t*h; /* 尽量减少乘法次数*/ s = t; do { t2 = 0.0; for(k=0; k x = a+(k+0.5)*h; /* 累加计算*/ t2 = t2+(*f)(x); } t2 = (t+h*t2)/2.0; s2 = (4.0*t2-t)/3.0; /* 计算新的积分值*/ t = t2; s = s2; /* 更新积分值*/ h = h/2.0; /* 更新步长*/ n = 2*n; /* 更新划分数*/ }while((fabs(s2-s)>1e-8)); return(s); } double cal_c(double zz,double D,double r) { double seigerma=3.2536e-10,rad=0.300e-10,k=1.380658e-23,ipuseilong=152.68*k,T=273.15; double c,derlta_4,derlta_5,derlta_10,derlta_11,oumiga_0,oumiga_1,oumiga; double f(double); derlta_4=1/4*(pow((1-r/D-rad/D),-4)-pow((1+r/D-rad/D),-4)); derlta_5=1/5*(pow((1-r/D-rad/D),-5)-pow((1+r/D-rad/D),-5)); derlta_10=1/10*(pow((1-r/D-rad/D),-10)-pow((1+r/D-rad/D),-10)); derlta_11=1/11*(pow((1-r/D-rad/D),-11)-pow((1+r/D-rad/D),-11)); oumiga_0=pow(seigerma,10)*(derlta_10+rad/D*derlta_11)/pow(D,11)/r; oumiga_1=pow(seigerma,6)*(derlta_4+rad/D*derlta_5)/pow(D,5)/r; oumiga=2*zz*ipuseilong*(oumiga_0-oumiga_1); c=simps1(0,D,f); /* 用变步长辛卜生求积法积分*/ c=4*3.1415/k/T*c; printf("simp: c=%2.5f\n",c); return(c); } main() { int i; double zz[2]={20,24},D[2]={3.95e-10,4.33e-10}; double c[2],r; for(i=0;i=1;i++) { c=cal_c(zz,D,r); printf("c[%d]=%.4f\n",i,c); } } #include "math.h" double f(r) double r; /* 被积分的函数*/ { double oumiga,k,T; double y; y=exp(-oumiga/k/T)*r*r; return(y); } 请问是什么问题，怎么解决。我在运行其他程序计算积分时没有这种问题 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
muqiuyan
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2963.5
- 散金: 355
- 红花: 1
- 帖子: 162
- 在线: 64.1小时
- 虫号: 829827
- 注册: 2009-08-18
- 性别: MM
- 专业: 能源化工
3楼2010-05-13 14:48:15
2楼2010-05-13 13:35:38
holmescn
金虫 (正式写手)
- 程序强帖: 37
- 应助: 1 (幼儿园)
- 金币: 1918.8
- 散金: 275
- 红花: 1
- 帖子: 699
- 在线: 102.6小时
- 虫号: 913482
- 注册: 2009-11-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2010-05-13 23:48:13
muqiuyan
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2963.5
- 散金: 355
- 红花: 1
- 帖子: 162
- 在线: 64.1小时
- 虫号: 829827
- 注册: 2009-08-18
- 性别: MM
- 专业: 能源化工
5楼2010-05-14 08:47:53
