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)
graphics(
    line_parametric_2d(fx, fy, ("t", 0, 12 * pi),
        use_cm=False, n=2000),
    title="Butterfly Curve",
)  # doctest: +SKIP
