| 查看: 1099 | 回复: 7 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】c++编程，请高手帮我看看这个程序应该怎么写
|
|||
|
程序中的公式为： p[I][j]=(p[j]+p[i-1][j])/2; s[I][j]=-2*miu*h*p[I][j]*p[I][j]*(1-p[I][j])/(4*dx*4*dx); fw[i-1][j]=((k[j]+k[i-1][j])/2*v[I][j]+(k[j]+k[i-2][j])/2*v[I-1][j])/2; fe[j]=((k[i+1][j]+k[j])/2*v[I+1][j]+(k[j]+k[i-1][j])/2*v[I][j])/2; aw[I-1][j]=dw[i-1][j]+fw[i-1][j]/2; ae[I+1][j]=de[j]-fe[j]/2; 在这个程序中包含i，j，I,J。我可以像下边这样写程序吗？ for(I=3;I<=m;I++) for(J=3;J<=n;J++) { for(i=2;i<=m;i++) for(j=2;j<=n;j++) { p[I][j]=(p[j]+p[i-1][j])/2; s[I][j]=-2*miu*h*p[I][j]*p[I][j]*(1-p[I][j])/(4*dx*4*dx); fw[i-1][j]=((k[j]+k[i-1][j])/2*v[I][j]+(k[j]+k[i-2][j])/2*v[I-1][j])/2; fe[j]=((k[i+1][j]+k[j])/2*v[I+1][j]+(k[j]+k[i-1][j])/2*v[I][j])/2; aw[I-1][j]=dw[i-1][j]+fw[i-1][j]/2; ae[I+1][j]=de[j]-fe[j]/2; } } [ Last edited by 飘雪139 on 2008-3-14 at 20:47 ] |
回复此楼» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有37人回复
-
好奇怪的filecode
已经有7人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有43人回复
-
filecode
已经有9人回复
-
请教兼职经验
已经有3人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有5人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有13人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有11人回复
-
读博的好处
已经有4人回复
-
filecode
已经有17人回复
6楼2008-03-15 18:09:06
2楼2008-03-15 08:53:11
3楼2008-03-15 17:46:06
|
for(I=3;I<=m;I++) for(J=3;J<=n;J++) { for(i=2;i<=m;i++) for(j=2;j<=n;j++) { p[I][j]=(p[j]+p[i-1][j])/2; s[I][j]=-2*miu*h*p[I][j]*p[I][j]*(1-p[I][j])/(4*dx*4*dx); dw[i-1][j]=miu*k[i-1][j]/dx; de[j]=miu*k[j]/dx; fw[i-1][j]=((k[j]+k[i-1][j])/2*v[I][j]+(k[j]+k[i-2][j])/2*v[I-1][j])/2； fe[j]=((k[i+1][j]+k[j])/2*v[I+1][j]+(k[j]+k[i-1][j])/2*v[I][j])/2; aw[I-1][j]=dw[i-1][j]+fw[i-1][j]/2; ae[I+1][j]=de[j]-fe[j]/2; au[I][j]=(k[j]+k[i-1][j])/2; } } |
4楼2008-03-15 17:54:57
