楕円 r=a(1-e^2)/(1+e*cos(theta)) の方程式を使用して、オブジェクトの楕円パスをプロットするコードを作成しようとしています。また、このデータを他の用途のために配列に入れたいと思います。
from numpy import *#Imports Python mathematical functions library
import matplotlib.pyplot as plt #Imports plot library
from pylab import *
a = 5
e = 0.3
theta = 0
while theta <= 2*pi:
r = (a*(1-e**2))/(1+e*cos(theta))
print("r = ",r,"theta = ",theta)
plt.polar(theta, r)
theta += pi/180
plt.show()
コードは r と theta の正しい値を吐き出しますが、プロットは空白です。極座標プロット ウィンドウが表示されますが、何もプロットされていません。
助けてください。前もって感謝します。