行星轨迹制作_行星运动轨迹的程序实现

//三体问题的导数曲线

class ThreeBody : publicDifferentiationFunction

{public:

ThreeBody()

{

m_star1_x= -10.0f;

m_star1_y= 0.0f;

m_star1_z= 0.0f;

m_star1_weight= 1.0f;

m_star2_x= 10.0f;

m_star2_y= 0.0f;

m_star2_z= 0.0f;

m_star2_weight= 1.0f;

m_planet_x= 5.0f;

m_planet_y= 5.0f;

m_planet_z= 0.1f;

m_planet_weight= 0.1f;

m_planet_speed_x= 0.0f;

m_planet_speed_y= 2.0f;

m_planet_speed_z= 0.0f;

m_g= 50.0f;

m_ek= 0.5f*m_planet_weight*(m_planet_speed_x*m_planet_speed_x + m_planet_speed_y*m_planet_speed_y + m_planet_speed_z*m_planet_speed_z); //1/2*m*v*v

float d1x = m_star1_x -m_planet_x;float d1y = m_star1_y -m_planet_y;float d1z = m_star1_z -m_planet_z;float sqd1 = d1x*d1x