三体问题数值解法源代码

/***************************************************************** 拿来检验算法的 进行10000次迭代大约耗时60秒（CPU：PentiumIII E 864MHz,256 MB SDRAM,很老的机器了……） 估算了一下，要模拟太阳、地球、金星三个天体在一年内的运动情况（间隔一秒）， 所需耗时达数月之久……囧…… 看来pc和巨型计算机还是有不小的差距的 由于在下不是计算机专业的，程序写得很烂…… 在vc6.0下编译通过 ******************************************************************/ #include"element.h" #include"TriObject.h" int main(void) { TriObject triobject; double A[3]={1.9891e+30,1.90e+27,5.69e+26}; triobject.Initialize(6.67e-11,A); triobject.SetTime(500,1e+7); double I[2][9]={1.25e+3,-2.20e+3,5.56e+3,4.49e+11,4.00e+11,-5.0e+11,8.25e+11,-7.95e+11,8.13e+11, 1.10e+4,9.70e+3,9.99e+3,7.85e+3,6.63e+3,7.23e+3,5.56e+3,5.40e+3,6.00e+3}; Matrix Info(2,9); for(unsigned i=0;i<=1;i++) for(unsigned j=0;j<=8;j++) Info.MSetAt(i,j,I[i][j]); //triobject.Format('m','k'); triobject.Emulate(Info); return 1; }