| 查看: 1661 | 回复: 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); } 请问是什么问题，怎么解决。我在运行其他程序计算积分时没有这种问题 |
回复此楼» 猜你喜欢
能否退出参与的面上项目解除限项
已经有18人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有3人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
-
我面上完蛋了
已经有9人回复
-
2026国自然函评费到账
已经有18人回复
-
建议基金发布提前给出明确的时间点
已经有15人回复
-
范进中举一文的中心思想
已经有7人回复
-
朋友圈看到的
已经有9人回复
-
让我中一个面上吧！
已经有16人回复
2楼2010-05-13 13:35:38
muqiuyan
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2963.5
- 散金: 355
- 红花: 1
- 帖子: 162
- 在线: 64.1小时
- 虫号: 829827
- 注册: 2009-08-18
- 性别: MM
- 专业: 能源化工
3楼2010-05-13 14:48:15
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
holmescn
金虫 (正式写手)
- 程序强帖: 37
- 应助: 1 (幼儿园)
- 金币: 1918.8
- 散金: 275
- 红花: 1
- 帖子: 699
- 在线: 102.6小时
- 虫号: 913482
- 注册: 2009-11-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2010-05-14 11:05:23
muqiuyan
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2963.5
- 散金: 355
- 红花: 1
- 帖子: 162
- 在线: 64.1小时
- 虫号: 829827
- 注册: 2009-08-18
- 性别: MM
- 专业: 能源化工
7楼2010-05-15 10:04:33
muqiuyan
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2963.5
- 散金: 355
- 红花: 1
- 帖子: 162
- 在线: 64.1小时
- 虫号: 829827
- 注册: 2009-08-18
- 性别: MM
- 专业: 能源化工
8楼2010-05-15 10:17:14
wenbo985
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 653.3
- 帖子: 14
- 在线: 3.5小时
- 虫号: 566193
- 注册: 2008-05-30
- 性别: GG
- 专业: 粒子物理学和场论
9楼2010-05-15 14:07:35
holmescn
金虫 (正式写手)
- 程序强帖: 37
- 应助: 1 (幼儿园)
- 金币: 1918.8
- 散金: 275
- 红花: 1
- 帖子: 699
- 在线: 102.6小时
- 虫号: 913482
- 注册: 2009-11-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2010-05-15 14:47:56
20