Kepler運動と軌跡

Keplerの方程式
     E=e*sin(E)+M
  の解 E を用いて
    x=a*cos(E)
   y=a*sin(E)*√(1-e2)
 を描画してみます。aを長半径 e (0〜1)を離心率といい各惑星固有の定数です。
e= ,M=として
M を連続的に変化して,太陽の周りを公転させてます。