import numpy as np
import matplotlib.pyplot as plt



a=np.linspace(0,2*np.pi)
x=16*(np.sin(a))**3
y=13*np.cos(a)-5*np.cos(2*a)-2*np.cos(3*a)-np.cos(4*a)

plt.plot(x,y,'r')
plt.axis([-25,25,-20,15])
plt.show ()

输出结果：