| 查看: 396 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
-
现代”学阀”该如何界定
已经有14人回复
-
求合成方法
已经有8人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有4人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
我的奶奶
已经有3人回复