| a = kwargs.pop('a',1)
|
| if not kwargs.pop('center', False):
|
| - return a*(1−np.e**2)/(1+np.e*np.cos(phi))
|
| + return a*(1-np.e**2)/(1+np.e*np.cos(phi))
|
| else:
|
| return pow( a**2*(1−np.e**2)/(1−np.e**2*cos(pi)**2), 1/2 )
|
| a = kwargs.pop('a',1)
|
| if not kwargs.pop('center', False):
|
| - return a*(1−np.e**2)/(1+np.e*np.cos(phi))
|
| + return a*(1-np.e**2)/(1+np.e*np.cos(phi))
|
| else:
|
| return pow( a**2*(1−np.e**2)/(1−np.e**2*cos(pi)**2), 1/2 )
|