| 查看: 348 | 回复: 9 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖 | |||
lansesanji银虫 (小有名气)
|
[交流]
c语言求助
|
||
| 各位精通C语言的大虾,我有个问题,我的程序定义了将近10个数组,但是数组的维数如果定义大了就不能正常运行,维数如果小的话就没有问题,到底是怎么回事啊?帮帮忙啊,不胜感激啊…… |
» 猜你喜欢
孩子确诊有中度注意力缺陷
已经有12人回复
2025冷门绝学什么时候出结果
已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
AI论文写作工具:是科研加速器还是学术作弊器?
已经有3人回复
2026博士申请-功能高分子,水凝胶方向
已经有6人回复
论文投稿,期刊推荐
已经有4人回复
硕士和导师闹得不愉快
已经有13人回复
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
同一篇文章,用不同账号投稿对编辑决定是否送审有没有影响?
已经有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













回复此楼