Написать код для отображения гипоциклоиды

import math
import matplotlib.pyplot as plt
import numpy as np
coordinate = int(input('ВВедите 0, если график необходим в декартовой системе, 1 - в полярной: '))
if coordinate == 0:
  def filling():
    while True:
      try:
        R = float(input('Введите значение радиуса R:'))
        break
      except:
        print('Введите число')
    return R
  def xy_astroida(R):
    x, y = [], []
    t = np.arange(0 * np.pi, 2 * np.pi, 0.01)
    x.append(R * np.cos(t)**3)
    y.append(R * np.sin(t)**3)
    return x[0], y[0]
  R = filling()
  x, y = xy_astroida(R)
  plt.plot(x, y)
  plt.show()
elif coordinate == 1:
  alpha = int(input('Введите значение альфа: '))
  theta = np.arange(0.,2.,1./180.)*np.pi
  r = alpha * ((np.cos(theta) ** (6) + np.sin(theta) ** 6)) ** (1/2)
  plt.polar(theta, r)
  plt.show()
else:
  print('Несуществующая система координат')

import numpy as np
from matplotlib import pyplot as plt
R=float(input('R: '))
r=float(input('r: '))
d=float(input('d: '))
t=np.linspace(0,100, 10000)
x=(R-r)*np.cos(t)+d*np.cos((R/r-1)*t)
y=(R-r)*np.sin(t)-d*np.sin((R/r-1)*t)
plt.plot(x,y)
plt.axis('square')
#plt.savefig('hypo.png')
plt.show()

R_=np.sqrt(x*x+y*y)
phi=np.arctan2(y, x)

plt.polar(phi, R_)