2016-05-09から1日間の記事一覧

ピタゴラス3体問題2

Δt 時間後の m3、m4 の位置、速度を求める函数。 m5 の位置は、微分方程式を解かなくても m3、m4 の位置から一意に求まるので別に計算する。 Define pyt(x30,vx30,y30,vy30,x40,vx40,y40,vy40,dt,tol)= Func :Local r34_3,r35_3,r45_3 :Define r34_3=√((x3-x…