| 查看: 377 | 回复: 0 | |||
| 当前主题已经存档。 | |||
jwf633木虫 (正式写手)
|
[交流]
请大家帮忙看看，我程序中的龙格-库塔方法对吗？
|
||
|
公式见附件： 程序主要代码如下: for(k=0;k extract1(i,j,xtemp); //pick out X[k][j] at site(i,j) to xtemp[k]; //caculate the k1[S] for all species at site(i,j) f(sigma[k],rou[k],kc,xtemp,k1); //caculate X[k][j]+k1[k]*tao/2 increase(tao,k1,xtemp); //caculate K2 for species k at site(i,j) f(sigma[k],rou[k],kc,xtemp,k2); //caculate X[k][j]+k2[k][j]*tao/2 increase(tao,k2,xtemp); //caculate k3 for species k at site(i,j) f(sigma[k],rou[k],kc,xtemp,k3); //caculate X[k][j]+k3[k][j]*tao increase(2*tao,k3,xtemp); //caculate K4[k][j] for species k at site(i,j) f(sigma[k],rou[k],kc,xtemp,k4); XX[k][j]=X[k][j]+tao*((k1[k]+2*k2[k]+2*k3[k]+k4[k])/6+PXX[k][j]); } void f(float sg,float r, float k[S][S],float x[S],float kk[S]) { int i,j; float temp1,temp2; for(i=0;i temp1=0; temp2=0; kk=(-1)*sg*x; for(j=0;j temp1=1-temp1; for(j=0;j temp2+=k[j]*x[j]; temp2+=r; kk+=temp1*temp2*x; } } [ Last edited by csfn on 2008-12-29 at 20:08 ] |
回复此楼» 猜你喜欢
国自然评分
已经有5人回复
-
国自然上会要求
已经有11人回复
-
投稿求助，期刊
已经有6人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有15人回复
-
提交了我也来说说感想
已经有10人回复
-
评审有感
已经有26人回复
-
今年审到国自然15份，谈谈感受
已经有22人回复
-
Sci. Bull. 悲剧经验
已经有8人回复
-
河北省自然科学基金
已经有4人回复
-
投MDPI的三个期刊都被桌拒
已经有4人回复
5