import numpy as np
from matplotlib import pyplot as plt

angle = np.linspace(0, 2 * np.pi, 100)
x_val = 16 * (np.sin(angle) ** 3)
y_val = 13 * np.cos(angle) - 5 * np.cos(2 * angle) - 2 * np.cos(3 * angle) - np.cos(4 * angle)

plt.plot(x_val, y_val)
plt.show()