import numpy as np
fx = lambda t: np.sin(t) * (np.exp(np.cos(t)) - 2 * np.cos(4 * t) - np.sin(t / 12)**5)
fy = lambda t: np.cos(t) * (np.exp(np.cos(t)) - 2 * np.cos(4 * t) - np.sin(t / 12)**5)
p = plot_parametric(fx, fy, ("t", 0, 12 * pi),
    title="Butterfly Curve", use_cm=False, n=2000)