| 查看: 1362 | 回复: 14 | ||
[交流]
【求助】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); } 请问是什么问题，怎么解决。我在运行其他程序计算积分时没有这种问题 |
回复此楼
jjdg
版主 (知名作家)
- 程序强帖: 9
- 应助: 223 (大学生)
- 贵宾: 3.926
- 金币: 51005.2
- 红花: 238
- 帖子: 7227
- 在线: 474.2小时
- 虫号: 27000
- 注册: 2003-10-21
- 专业: 药物代谢与药物动力学
- 管辖: 程序语言
2楼2010-05-13 13:35:38
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
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
7楼2010-05-15 10:04:33
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
