Butterfly Curve

 import numpy as np 

import matplotlib.pyplot as plt 

from math import *

  

  

 

plt.axes(projection = 'polar') 

  

rads = np.arange(0, (2 * np.pi), 0.01) 

  

for rad in rads: 

r=np.exp(sin(rad))-2*cos(4*rad)+sin((2*rad-np.pi)/24)**5

plt.polar(rad, r, 'g.')

plt.show()