| 查看: 430 | 回复: 9 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
lansesanji银虫 (小有名气)
|
[交流]
c语言求助
|
||
| 各位精通C语言的大虾，我有个问题，我的程序定义了将近10个数组，但是数组的维数如果定义大了就不能正常运行，维数如果小的话就没有问题，到底是怎么回事啊？帮帮忙啊，不胜感激啊…… |
回复此楼» 猜你喜欢
293求调剂
已经有14人回复
-
0703化学336分求调剂
已经有6人回复
-
304求调剂
已经有9人回复
-
281求调剂（0805）
已经有8人回复
-
环境领域全国重点实验室招收博士1-2名
已经有3人回复
-
材料专硕306英一数二
已经有10人回复
-
301求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有7人回复
-
302求调剂
已经有6人回复
-
26博士申请
已经有3人回复
lansesanji
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1287
- 散金: 349
- 红花: 1
- 帖子: 124
- 在线: 21.6小时
- 虫号: 606371
- 注册: 2008-09-18
- 专业: 油气储运工程
7楼2009-12-23 22:42:32
shenao_07 
木虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 1581.3
- 散金: 726
- 红花: 2
- 帖子: 1350
- 在线: 146.9小时
- 虫号: 926332
- 注册: 2009-12-12
- 性别: GG
- 专业: 能源化工
2楼2009-12-23 17:54:05
sun_sun_sun
金虫 (小有名气)
- 博学EPI: 1
- 应助: 0 (幼儿园)
- 金币: 777.8
- 帖子: 135
- 在线: 28.5小时
- 虫号: 874679
- 注册: 2009-10-16
- 性别: GG
- 专业: 电气
3楼2009-12-23 18:23:53
lansesanji
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1287
- 散金: 349
- 红花: 1
- 帖子: 124
- 在线: 21.6小时
- 虫号: 606371
- 注册: 2008-09-18
- 专业: 油气储运工程
现在贴上源程序，大家帮帮忙了！
|
/* Note:Your choice is C IDE */ #include #include #define k 6 #define numb 67 void main() { double T[numb][numb],detar,detax,aw,ae,an,as,ap,xw,xe,rn,rs,e1=0,e2=0,b[numb][numb],a[numb][numb]={0}; double r[numb],x[numb],t[numb][numb]; double T0[7][numb][numb]={0},a0[7][numb][numb]={0}; int i,j,m,n,Nk; Nk=pow(2,k)+1; for(i=0;i for(j=0;j for(j=1;j<=Nk;j++) T[j][0]=400; for(i=1;i } for(i=0;i<=Nk;i++) {for(j=0;j<=Nk;j++) T0[1][j]=T[j]; } detar=1.0/Nk; detax=1.0/Nk; for(i=0;i<=Nk;i++) r=x=detar*i; do{ e1=0; e2=0; for(i=0;i<=Nk;i++) {for(j=0;j<=Nk;j++) t[j]=T[j]; } for(n=1;n<7;n++) { for(m=0;m<3;m++) {for(i=1;i xe=x[i+1]-x; rn=r[j+1]-r[j]; rs=r[j]-r[j-1]; ap=2/xw/xe+2/rs/rn+1/r[j]/rn; aw=2/(xw*xw+xe*xw); ae=2/(xe*xe+xe*xw); as=2/(rs*rs+rs*rn); an=2/(rn*rn+rs*rn)+1/r[j]/rn; T[j]=aw/ap*T[i-1][j]+as/ap*T[j-1]+ae/ap*T[i+1][j]+an/ap*T[j+1]+a0[n-1][j]/ap; T0[n][j]=T[j]; } } } for(i=1;i xe=x[i+1]-x; rn=r[j+1]-r[j]; rs=r[j]-r[j-1]; ap=2/xw/xe+2/rs/rn+1/r[j]/rn; aw=2/(xw*xw+xe*xw); ae=2/(xe*xe+xe*xw); as=2/(rs*rs+rs*rn); an=2/(rn*rn+rs*rn)+1/r[j]/rn; b[j]=a0[n-1][j]-ap*T[j]+aw*T[i-1][j]+as*T[j-1]+ae*T[i+1][j]+an*T[j+1]; } } Nk=(Nk+1)/2; r[0]=x[0]=0; r[1]=x[1]=pow(2,(n-2))*detar+0.5*detar; for(i=2;i r[Nk]=x[Nk]=1; for(i=1;i a0[n][j]=a[j]; } } for(i=0;i<=Nk;i++) {for(j=0;j<=Nk;j++) T[j]=0; } } Nk=Nk*2-1; for(n=5;n>0;n--) { Nk=Nk*2-1; if (Nk==pow(2,k)+1) for(i=0;i<=Nk;i++) r=x=detar*i; else {r[0]=x[0]=0; r[1]=x[1]=pow(2,(n-2))*detar+0.5*detar; for(i=2;i r[Nk]=x[Nk]=1; } for(i=1;i<(Nk/2);i++) {for(j=1;j<(Nk/2);j++) {T0[n][2*i-1][2*j]=T0[n+1][j]+T0[n][2*i-1][2*j]; T0[n][2*i-1][2*j-1]=T0[n+1][j]+T0[n][2*i-1][2*j-1]; T0[n][2*i][2*j]=T0[n+1][j]+T0[n][2*i][2*j]; T0[n][2*i][2*j-1]=T0[n+1][j]+T0[n][2*i][2*j-1]; } } for(m=0;m<3;m++) {for(i=1;i xe=x[i+1]-x; rn=r[j+1]-r[j]; rs=r[j]-r[j-1]; ap=2/xw/xe+2/rs/rn+1/r[j]/rn; aw=2/(xw*xw+xe*xw); ae=2/(xe*xe+xe*xw); as=2/(rs*rs+rs*rn); an=2/(rn*rn+rs*rn)+1/r[j]/rn; T0[n][j]=aw/ap*T0[n][i-1][j]+as/ap*T0[n][j-1]+ae/ap*T0[n][i+1][j]+an/ap*T0[n][j+1]+a0[n][j]/ap; } } } for(i=0;i<=Nk;i++) {for(j=0;j<=Nk;j++) T[j]=T0[1][j]; } } for(i=0;i<=Nk;i++) {for(j=0;j<=Nk;j++) {e2=fabs(t[j]-T[j]); if (e2>e1) e1=e2; } } }while (e1>0.001); for(i=0;i<=Nk;i++) {for(j=0;j<=Nk;j++) (printf("T[%d][%d]=%lf ",i,j,T[j])); } } |
4楼2009-12-23 19:43:43
