| 查看: 374 | 回复: 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 ] |
回复此楼» 猜你喜欢
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
-
300求调剂，材料科学英一数二
已经有8人回复
-
招收生物学/细胞生物学调剂
已经有5人回复
-
070305高分子化学与物理 304分求调剂
已经有7人回复
-
289求调剂
已经有13人回复
-
一志愿哈尔滨工业大学材料与化工方向336分
已经有9人回复
-
081200-11408-276学硕求调剂
已经有6人回复
-
调剂求院校招收
已经有5人回复
-
调剂310
已经有8人回复
