terrible attempt of plotting a Likitung parametric curve

likitung.py

import numpy as np
import matplotlib.pyplot as plt
import time

def H(x):
    return 0.5 * (np.sign(x) + 1)

def x(t):
    return ((-15/16 * np.sin(3/2 - 25 * t) - 5/7 * np.sin(16/11 - 23 * t) - 12/7 * np.sin(3/2 - 21 * t) - np.sin(14/9 - 19 * t) - 4/5 * np.sin(11/7 - 18 * t) - 7/5 * np.sin(3/2 - 17 * t) - 49/16 * np.sin(14/9 - 15 * t) - 7/5 * np.sin(11/7 - 14 * t) - 22/9 * np.sin(14/9 - 13 * t) - 7/4 * np.sin(14/9 - 12 * t) - 79/7 * np.sin(14/9 - 10 * t) - 47/8 * np.sin(17/11 - 8 * t) - 52/5 * np.sin(11/7 - 6 * t) - 108/5 * np.sin(11/7 - 5 * t) + 643/6 * np.sin(t + 11/7) + 431/18 * np.sin(2 * t + 11/7) + 111/5 * np.sin(3 * t + 8/5) + 179/8 * np.sin(4 * t + 8/5) + 152/17 * np.sin(7 * t + 8/5) + 5/6 * np.sin(9 * t + 7/4) + 22/9 * np.sin(11 * t + 13/8) + 1/18 * np.sin(16 * t + 12/7) + 2/5 * np.sin(20 * t + 14/3) + 18/17 * np.sin(22 * t + 5/3) + 1/8 * np.sin(24 * t + 33/16) - 64) * H(67 * np.pi - t) * H(t - 63 * np.pi) + (-1/2 * np.sin(3/2 - 24 * t) - 19/9 * np.sin(14/9 - 16 * t) - 6/11 * np.sin(10/7 - 15 * t) - 53/6 * np.sin(14/9 - 10 * t) - 40/7 * np.sin(14/9 - 9 * t) - 85/6 * np.sin(11/7 - 4 * t) - 141 * np.sin(11/7 - t) + 59/4 * np.sin(2 * t + 11/7) + 25/4 * np.sin(3 * t + 8/5) + 17/5 * np.sin(5 * t + 11/7) + 41/4 * np.sin(6 * t + 8/5) + 102/5 * np.sin(7 * t + 8/5) + 38/5 * np.sin(8 * t + 8/5) + 34/9 * np.sin(11 * t + 8/5) + 13/2 * np.sin(12 * t + 8/5) + 22/5 * np.sin(13 * t + 8/5) + 2/5 * np.sin(14 * t + 13/8) + 4/3 * np.sin(17 * t + 8/5) + 7/4 * np.sin(18 * t + 8/5) + 5/4 * np.sin(19 * t + 8/5) + 3/4 * np.sin(20 * t + 14/3) + 8/7 * np.sin(21 * t + 8/5) + 1/3 * np.sin(22 * t + 5/3) + 1/5 * np.sin(23 * t + 11/7) + 676/3) * H(63 * np.pi - t) * H(t - 59 * np.pi) + (-21/5 * np.sin(11/7 - 5 * t) - 68/5 * np.sin(11/7 - 3 * t) + 307/2 * np.sin(t + 8/5) + 45/7 * np.sin(2 * t + 13/8) + 29/8 * np.sin(4 * t + 3/2) + 31/6 * np.sin(6 * t + 8/5) + 13/9 * np.sin(7 * t + 14/3) + 61/11 * np.sin(8 * t + 8/5) + 5/6 * np.sin(9 * t + 41/9) + 13/4 * np.sin(10 * t + 13/8) + 1/7 * np.sin(11 * t + 26/7) + 18/17 * np.sin(12 * t + 8/5) - 61/6) * H(59 * np.pi - t) * H(t - 55 * np.pi) + (11/3 * np.sin(t + 35/8) + 1/7 * np.sin(2 * t + 4/5) + 1/2 * np.sin(3 * t + 23/6) + 1/9 * np.sin(4 * t + 3/4) + 1/5 * np.sin(5 * t + 7/2) + 1/11 * np.sin(6 * t + 2/3) + 1/7 * np.sin(7 * t + 10/3) + 1/12 * np.sin(8 * t + 4/7) + 1/11 * np.sin(9 * t + 13/4) + 1/14 * np.sin(10 * t + 5/11) + 1/14 * np.sin(11 * t + 16/5) + 1/16 * np.sin(12 * t + 2/5) - 480/7) * H(55 * np.pi - t) * H(t - 51 * np.pi) + (-1/31 * np.sin(3/5 - 12 * t) - 1/31 * np.sin(5/7 - 10 * t) - 1/32 * np.sin(4/5 - 8 * t) + 3 * np.sin(t + 14/9) + 1/4 * np.sin(2 * t + 13/14) + 2/5 * np.sin(3 * t + 29/14) + 1/36 * np.sin(4 * t + 3/7) + 1/6 * np.sin(5 * t + 5/2) + 1/11 * np.sin(7 * t + 14/5) + 1/18 * np.sin(9 * t + 3) + 1/27 * np.sin(11 * t + 22/7) - 298/3) * H(51 * np.pi - t) * H(t - 47 * np.pi) + (-5/7 * np.sin(3/5 - 27 * t) - 4/5 * np.sin(4/3 - 25 * t) - 1/7 * np.sin(1/3 - 23 * t) - 3/5 * np.sin(4/5 - 20 * t) - 5/7 * np.sin(1/5 - 19 * t) - 9/5 * np.sin(3/2 - 16 * t) - 26/5 * np.sin(1 - 7 * t) - 17/4 * np.sin(19/18 - 6 * t) - 49/10 * np.sin(1 - 5 * t) - 76/5 * np.sin(3/5 - 4 * t) + 443/6 * np.sin(t + 11/4) + 293/4 * np.sin(2 * t + 22/5) + 306/7 * np.sin(3 * t + 9/5) + 4 * np.sin(8 * t + 19/5) + 35/6 * np.sin(9 * t + 13/8) + 31/7 * np.sin(10 * t + 3/4) + 9/2 * np.sin(11 * t + 1/12) + 13/3 * np.sin(12 * t + 9/2) + np.sin(13 * t + 20/7) + 13/9 * np.sin(14 * t + 35/8) + 2/3 * np.sin(15 * t + 8/3) + 12/7 * np.sin(17 * t + 33/16) + 27/13 * np.sin(18 * t + 2/7) + 19/9 * np.sin(21 * t + 20/7) + 9/5 * np.sin(22 * t + 8/5) + 1/2 * np.sin(24 * t + 2/7) + 5/9 * np.sin(26 * t + 7/4) + 1423/5) * H(47 * np.pi - t) * H(t - 43 * np.pi) + (-10/9 * np.sin(1/6 - 13 * t) - 6/5 * np.sin(1/13 - 11 * t) - 5/4 * np.sin(4/7 - 9 * t) - 6/5 * np.sin(3/5 - 7 * t) - 1/2 * np.sin(7/6 - 6 * t) - 32/9 * np.sin(3/8 - 5 * t) - 11/4 * np.sin(7/6 - 3 * t) + 207/2 * np.sin(t + 25/13) + 23/7 * np.sin(2 * t + 4) + 17/5 * np.sin(4 * t + 22/5) + 17/7 * np.sin(8 * t + 16/7) + 13/4 * np.sin(10 * t + 20/7) + 16/7 * np.sin(12 * t + 17/5) + 32/31 * np.sin(14 * t + 11/3) + 7/5 * np.sin(15 * t + 1/9) + 3/7 * np.sin(16 * t + 10/3) + 16/15 * np.sin(17 * t + 1/2) + 4/7 * np.sin(18 * t + 14/5) + 3/7 * np.sin(19 * t + 4/5) - 65/3) * H(43 * np.pi - t) * H(t - 39 * np.pi) + (-141/28 * np.sin(3/2 - 3 * t) + 199/3 * np.sin(t + 33/7) + 1/3 * np.sin(2 * t + 19/7) + 9/7 * np.sin(4 * t + 18/17) + 11/4 * np.sin(5 * t + 37/8) + 4/5 * np.sin(6 * t + 5/6) + 7/6 * np.sin(7 * t + 17/4) + 2/7 * np.sin(8 * t + 12/13) + 2/3 * np.sin(9 * t + 35/8) + 1/6 * np.sin(10 * t + 6/5) + 3/8 * np.sin(11 * t + 30/7) + 1/7 * np.sin(12 * t + 5/4) - 746/5) * H(39 * np.pi - t) * H(t - 35 * np.pi) + (-1/6 * np.sin(2/3 - 21 * t) - 2/5 * np.sin(1/15 - 13 * t) - 9/10 * np.sin(1/6 - 11 * t) + 5/3 * np.sin(t + 38/11) + 11/7 * np.sin(2 * t + 41/9) + 11/6 * np.sin(3 * t + 49/11) + 17/6 * np.sin(4 * t + 39/10) + 6/7 * np.sin(5 * t + 9/4) + 74/15 * np.sin(6 * t + 2/5) + 6/7 * np.sin(7 * t + 2) + 1/4 * np.sin(8 * t + 1/6) + 11/8 * np.sin(9 * t + 1/4) + 1/6 * np.sin(10 * t + 7/4) + 2/7 * np.sin(12 * t + 13/6) + 1/4 * np.sin(14 * t + 5/7) + 1/24 * np.sin(15 * t + 27/7) + 2/5 * np.sin(16 * t + 20/7) + 1/5 * np.sin(17 * t + 7/4) + 1/8 * np.sin(18 * t + 13/8) + 1/6 * np.sin(19 * t + 53/18) + 1/8 * np.sin(20 * t + 13/8) + 1/4 * np.sin(22 * t + 9/5) + 1/10 * np.sin(23 * t + 13/3) - 648/5) * H(35 * np.pi - t) * H(t - 31 * np.pi) + (-6/5 * np.sin(3/2 - 10 * t) + 481/5 * np.sin(t + 8/5) + 142/13 * np.sin(2 * t + 23/5) + 5 * np.sin(3 * t + 55/27) + 23/6 * np.sin(4 * t + 17/4) + 22/9 * np.sin(5 * t + 7/3) + 9/7 * np.sin(6 * t + 29/7) + 14/9 * np.sin(7 * t + 9/4) + 2/3 * np.sin(8 * t + 23/5) + 1/3 * np.sin(9 * t + 14/5) + 4/9 * np.sin(11 * t + 63/16) + 4/5 * np.sin(12 * t + 14/3) + 47/3) * H(31 * np.pi - t) * H(t - 27 * np.pi) + (-2 * np.sin(5/6 - 7 * t) - 166/11 * np.sin(4/5 - 2 * t) + 518/5 * np.sin(t + 23/12) + 32/7 * np.sin(3 * t + 8/3) + 25/13 * np.sin(4 * t + 14/13) + 9/7 * np.sin(5 * t + 89/22) + 7/5 * np.sin(6 * t + 17/5) + 6/5 * np.sin(8 * t + 23/6) + 33/32 * np.sin(9 * t + 1/8) + 23/22 * np.sin(10 * t + 4) + 2/3 * np.sin(11 * t + 9/7) + 3/7 * np.sin(12 * t + 35/8) + 35) * H(27 * np.pi - t) * H(t - 23 * np.pi) + (-1/9 * np.sin(16/17 - 11 * t) - 5/7 * np.sin(2/3 - 7 * t) + 381/5 * np.sin(t + 8/3) + 21/4 * np.sin(2 * t + 4/5) + 23/6 * np.sin(3 * t + 1/5) + 23/7 * np.sin(4 * t + 25/6) + 5/3 * np.sin(5 * t + 23/8) + 7/5 * np.sin(6 * t + 1/15) + 4/7 * np.sin(8 * t + 11/3) + 2/5 * np.sin(9 * t + 12/7) + 2/5 * np.sin(10 * t + 1/4) + 1/3 * np.sin(12 * t + 3) + 77/3) * H(23 * np.pi - t) * H(t - 19 * np.pi) + (-33/8 * np.sin(1/2 - 6 * t) - 1651/11 * np.sin(1/5 - t) + 74/25 * np.sin(9 * t) + 77/4 * np.sin(2 * t + 9/4) + 253/12 * np.sin(3 * t + 7/4) + 27/8 * np.sin(4 * t + 17/4) + 59/9 * np.sin(5 * t + 21/8) + 7/2 * np.sin(7 * t + 26/7) + 17/6 * np.sin(8 * t + 17/9) + 7/3 * np.sin(10 * t + 23/5) + 3 * np.sin(11 * t + 13/5) + 11/5 * np.sin(12 * t + 3/7) - 377/3) * H(19 * np.pi - t) * H(t - 15 * np.pi) + (-4/5 * np.sin(11/12 - 12 * t) - 2/7 * np.sin(15/16 - 11 * t) - 26/27 * np.sin(5/9 - 10 * t) - 373/5 * np.sin(8/7 - t) + 170/13 * np.sin(2 * t + 13/7) + 11/3 * np.sin(3 * t + 7/5) + 25/4 * np.sin(4 * t + 17/9) + 23/7 * np.sin(5 * t + 5/4) + 19/9 * np.sin(6 * t + 16/11) + 27/13 * np.sin(7 * t + 8/7) + 13/14 * np.sin(8 * t + 4/5) + 4/5 * np.sin(9 * t + 5/6) - 401/2) * H(15 * np.pi - t) * H(t - 11 * np.pi) + (-2/3 * np.sin(4/5 - 19 * t) - 1/5 * np.sin(5/4 - 18 * t) - 4/3 * np.sin(1/12 - 13 * t) - 11/8 * np.sin(2/7 - 11 * t) - 19/6 * np.sin(6/7 - 7 * t) - 23/12 * np.sin(6/5 - 6 * t) - 23/6 * np.sin(3/8 - 4 * t) - 10/7 * np.sin(11/12 - t) + 2/3 * np.sin(2 * t + 17/6) + 11/6 * np.sin(3 * t + 11/7) + 2/3 * np.sin(5 * t + 21/5) + 33/7 * np.sin(8 * t + 33/7) + 7/3 * np.sin(9 * t + 19/7) + 32/7 * np.sin(10 * t + 14/3) + 43/21 * np.sin(12 * t + 10/7) + 5/7 * np.sin(14 * t + 19/9) + 9/7 * np.sin(15 * t + 35/8) + 3/4 * np.sin(16 * t + 11/5) + 3/7 * np.sin(17 * t + 1/4) + 3/5 * np.sin(20 * t + 11/5) + 1/4 * np.sin(21 * t + 3/4) + 1/4 * np.sin(22 * t + 16/5) + 1/3 * np.sin(23 * t + 11/6) + 2/7 * np.sin(24 * t + 10/3) + 201/5) * H(11 * np.pi - t) * H(t - 7 * np.pi) + (325/7 * np.sin(t + 3/8) + np.sin(2 * t + 11/4) + 7/5 * np.sin(3 * t + 2/3) + 3/5 * np.sin(4 * t + 1) + 3/5 * np.sin(5 * t + 18/5) + 769/5) * H(7 * np.pi - t) * H(t - 3 * np.pi) + (-5/6 * np.sin(2/5 - 2 * t) + 80/3 * np.sin(t + 1/2) + 1/3 * np.sin(3 * t + 5/4) + 3/5 * np.sin(4 * t + 6/5) + 5/9 * np.sin(5 * t + 45/11) + 1058/7) * H(3 * np.pi - t) * H(t + np.pi)) * H(np.sqrt(np.sign(np.sin(t/2))))

def y(t):
    return ((-1/3 * np.sin(7/5 - 25 * t) - 5/8 * np.sin(13/9 - 23 * t) - 1/5 * np.sin(13/9 - 17 * t) - 9/5 * np.sin(11/7 - 13 * t) - 17/6 * np.sin(14/9 - 12 * t) - 4/5 * np.sin(14/9 - 10 * t) - 17/5 * np.sin(14/9 - 9 * t) - 83/7 * np.sin(11/7 - 7 * t) - 88/5 * np.sin(11/7 - 3 * t) - 64/5 * np.sin(11/7 - 2 * t) + 339/13 * np.sin(t + 11/7) + 73/7 * np.sin(4 * t + 8/5) + 37/4 * np.sin(5 * t + 8/5) + 12/7 * np.sin(6 * t + 8/5) + 8/3 * np.sin(8 * t + 8/5) + 7/8 * np.sin(11 * t + 8/5) + 3/2 * np.sin(14 * t + 8/5) + 1/3 * np.sin(15 * t + 7/5) + 28/11 * np.sin(16 * t + 13/8) + 3/5 * np.sin(18 * t + 8/5) + 7/4 * np.sin(19 * t + 8/5) + 5/3 * np.sin(20 * t + 8/5) + 1/8 * np.sin(21 * t + 7/5) + 5/6 * np.sin(22 * t + 5/3) + 7/8 * np.sin(24 * t + 5/3) - 873/2) * H(67 * np.pi - t) * H(t - 63 * np.pi) + (-7/13 * np.sin(14/9 - 19 * t) - 5/6 * np.sin(3/2 - 17 * t) - 1/5 * np.sin(4/3 - 16 * t) - 1/4 * np.sin(10/7 - 12 * t) - 139/20 * np.sin(11/7 - 11 * t) - 39/8 * np.sin(11/7 - 7 * t) + 927/8 * np.sin(t + 33/7) + 2775/19 * np.sin(2 * t + 11/7) + 29/6 * np.sin(3 * t + 5/3) + 243/8 * np.sin(4 * t + 8/5) + 7 * np.sin(5 * t + 8/5) + 101/7 * np.sin(6 * t + 8/5) + 39/11 * np.sin(8 * t + 8/5) + 10/3 * np.sin(9 * t + 8/5) + 18/5 * np.sin(10 * t + 8/5) + 35/8 * np.sin(13 * t + 8/5) + 26/5 * np.sin(14 * t + 8/5) + 12/13 * np.sin(15 * t + 14/3) + 1/4 * np.sin(18 * t + 8/5) + 23/22 * np.sin(20 * t + 8/5) + 2/3 * np.sin(21 * t + 12/7) + 7/5 * np.sin(22 * t + 8/5) + 1/11 * np.sin(23 * t + 35/8) + 1/8 * np.sin(24 * t + 10/7) - 346) * H(63 * np.pi - t) * H(t - 59 * np.pi) + (-6/5 * np.sin(3/2 - 12 * t) - 1/4 * np.sin(6/5 - 11 * t) - 11/7 * np.sin(3/2 - 10 * t) - 7/3 * np.sin(11/7 - 8 * t) - 13/3 * np.sin(14/9 - 6 * t) - 11/5 * np.sin(3/2 - 5 * t) - 88/7 * np.sin(14/9 - 4 * t) - 3 * np.sin(11/8 - 3 * t) - 1847/22 * np.sin(11/7 - 2 * t) + 13/3 * np.sin(t + 13/9) + 2/5 * np.sin(7 * t + 7/5) + 4/7 * np.sin(9 * t + 3/2) + 1214/15) * H(59 * np.pi - t) * H(t - 55 * np.pi) + (37/7 * np.sin(t + 32/7) + 1/18 * np.sin(2 * t + 14/5) + 3/5 * np.sin(3 * t + 17/4) + 1/23 * np.sin(4 * t + 9/4) + 1/4 * np.sin(5 * t + 4) + 1/25 * np.sin(6 * t + 13/7) + 1/8 * np.sin(7 * t + 26/7) + 1/26 * np.sin(8 * t + 8/5) + 1/12 * np.sin(9 * t + 7/2) + 1/27 * np.sin(10 * t + 4/3) + 1/19 * np.sin(11 * t + 10/3) + 1/29 * np.sin(12 * t + 8/7) + 33) * H(55 * np.pi - t) * H(t - 51 * np.pi) + (-1/11 * np.sin(1/20 - 11 * t) - 1/7 * np.sin(1/16 - 9 * t) - 1/5 * np.sin(1/6 - 7 * t) - 1/3 * np.sin(2/5 - 5 * t) - 5/6 * np.sin(5/7 - 3 * t) - 34/5 * np.sin(6/5 - t) + 1/4 * np.sin(2 * t + 14/9) + 1/6 * np.sin(4 * t + 2) + 1/6 * np.sin(6 * t + 9/4) + 1/7 * np.sin(8 * t + 12/5) + 1/8 * np.sin(10 * t + 13/5) + 1/9 * np.sin(12 * t + 11/4) + 147/4) * H(51 * np.pi - t) * H(t - 47 * np.pi) + (-2/3 * np.sin(1/31 - 21 * t) - 26/9 * np.sin(3/5 - 13 * t) - 9/4 * np.sin(1/3 - 12 * t) - 27/7 * np.sin(5/4 - 11 * t) - 17/4 * np.sin(10/11 - 10 * t) - 124/7 * np.sin(5/4 - 4 * t) - 19 * np.sin(1/10 - 3 * t) - 241/3 * np.sin(1/33 - 2 * t) - 542/3 * np.sin(10/7 - t) + 129/10 * np.sin(5 * t + 8/3) + 37/9 * np.sin(6 * t + 7/4) + 30/7 * np.sin(7 * t + 19/5) + 11/6 * np.sin(8 * t + 18/5) + 3 * np.sin(9 * t + 14/3) + 11/10 * np.sin(14 * t + 1/25) + 10/11 * np.sin(15 * t + 2/3) + 8/5 * np.sin(16 * t + 4/5) + 6/7 * np.sin(17 * t + 1/5) + 1/2 * np.sin(18 * t + 2) + 8/7 * np.sin(19 * t + 1) + 4/5 * np.sin(20 * t + 2/5) + 1/2 * np.sin(22 * t + 28/11) + 7/5 * np.sin(23 * t + 8/5) + 12/11 * np.sin(24 * t + 6/5) + 2/7 * np.sin(25 * t + 3/4) + 2/3 * np.sin(26 * t + 11/4) + 5/8 * np.sin(27 * t + 11/5) - 1199/8) * H(47 * np.pi - t) * H(t - 43 * np.pi) + (-2/7 * np.sin(1/9 - 18 * t) - 12/5 * np.sin(2/3 - 8 * t) - 25/2 * np.sin(4/3 - 3 * t) + 188/7 * np.sin(t + 31/8) + 17/2 * np.sin(2 * t + 23/6) + 15/4 * np.sin(4 * t + 15/4) + 19/8 * np.sin(5 * t + 61/15) + 45/8 * np.sin(6 * t + 53/13) + 13/5 * np.sin(7 * t + 16/5) + 5/3 * np.sin(9 * t + 23/6) + 5/6 * np.sin(10 * t + 13/8) + np.sin(11 * t + 7/3) + 1/3 * np.sin(12 * t + 11/7) + 11/8 * np.sin(13 * t + 18/7) + 1/2 * np.sin(14 * t + 5/4) + 8/9 * np.sin(15 * t + 31/9) + 1/5 * np.sin(16 * t + 4/5) + 1/7 * np.sin(17 * t + 11/3) + 2/7 * np.sin(19 * t + 16/7) - 37/7) * H(43 * np.pi - t) * H(t - 39 * np.pi) + (290/3 * np.sin(t + 14/3) + 1/9 * np.sin(2 * t + 1/18) + 37/3 * np.sin(3 * t + 41/9) + 3/4 * np.sin(4 * t + 25/6) + 19/5 * np.sin(5 * t + 31/7) + 1/2 * np.sin(6 * t + 15/4) + 15/7 * np.sin(7 * t + 9/2) + 1/6 * np.sin(8 * t + 11/3) + 9/7 * np.sin(9 * t + 13/3) + 1/13 * np.sin(10 * t + 24/7) + 10/11 * np.sin(11 * t + 30/7) + 1/21 * np.sin(12 * t + 17/6) - 1185/7) * H(39 * np.pi - t) * H(t - 35 * np.pi) + (-1/4 * np.sin(5/8 - 23 * t) - 1/3 * np.sin(1/6 - 21 * t) - 1/5 * np.sin(3/5 - 17 * t) - 39/8 * np.sin(11/8 - 4 * t) - 1/10 * np.sin(3/2 - 2 * t) + 13/4 * np.sin(t + 17/4) + 15/4 * np.sin(3 * t + 22/5) + 15/8 * np.sin(5 * t + 11/3) + 55/7 * np.sin(6 * t + 7/5) + 7/6 * np.sin(7 * t + 9/4) + 2/5 * np.sin(8 * t + 13/9) + 8/5 * np.sin(9 * t + 4/3) + 4/5 * np.sin(10 * t + 10/3) + 13/12 * np.sin(11 * t + 13/12) + 2/5 * np.sin(12 * t + 11/4) + 7/8 * np.sin(13 * t + 17/18) + 2/5 * np.sin(14 * t + 11/5) + 1/4 * np.sin(15 * t + 1/3) + 3/4 * np.sin(16 * t + 129/32) + 1/22 * np.sin(18 * t + 15/4) + 1/5 * np.sin(19 * t + 2/3) + 1/4 * np.sin(20 * t + 3) + 1/3 * np.sin(22 * t + 3) + 495/7) * H(35 * np.pi - t) * H(t - 31 * np.pi) + (-3/7 * np.sin(1/8 - 8 * t) + 255/7 * np.sin(t + 16/7) + 157/6 * np.sin(2 * t + 7/4) + 24/5 * np.sin(3 * t + 14/5) + 32/7 * np.sin(4 * t + 15/8) + 23/7 * np.sin(5 * t + 16/5) + 8/7 * np.sin(6 * t + 7/5) + 27/11 * np.sin(7 * t + 19/5) + 11/7 * np.sin(9 * t + 18/5) + 3/8 * np.sin(10 * t + 2/5) + 5/7 * np.sin(11 * t + 24/7) + 3/5 * np.sin(12 * t + 1/14) - 4295/12) * H(31 * np.pi - t) * H(t - 27 * np.pi) + (-2 * np.sin(5/8 - 7 * t) - 38/7 * np.sin(3/2 - 5 * t) + 1231/16 * np.sin(t + 5/2) + 129/5 * np.sin(2 * t + 7/3) + 56/5 * np.sin(3 * t + 81/20) + 14/5 * np.sin(4 * t + 15/4) + 2 * np.sin(6 * t + 3/4) + 3/2 * np.sin(8 * t + 10/7) + 9/10 * np.sin(9 * t + 4/7) + 8/9 * np.sin(10 * t + 27/14) + 1/3 * np.sin(11 * t + 50/17) + 1/5 * np.sin(12 * t + 11/4) - 1208/5) * H(27 * np.pi - t) * H(t - 23 * np.pi) + (-6/11 * np.sin(4/7 - 11 * t) - 5/11 * np.sin(3/7 - 10 * t) - 3/2 * np.sin(1/4 - 7 * t) - 1/4 * np.sin(3/4 - 6 * t) - 46/5 * np.sin(1/3 - 3 * t) + 287/4 * np.sin(t + 7/2) + 17/3 * np.sin(2 * t + 18/5) + 11/3 * np.sin(4 * t + 9/5) + 13/5 * np.sin(5 * t + 10/3) + 23/22 * np.sin(8 * t + 29/14) + 7/8 * np.sin(9 * t + 10/3) + 2/5 * np.sin(12 * t + 9/4) - 1205/8) * H(23 * np.pi - t) * H(t - 19 * np.pi) + (-3/4 * np.sin(1/7 - 12 * t) - 13/5 * np.sin(5/8 - 9 * t) - 11/5 * np.sin(6/7 - 6 * t) - 54/5 * np.sin(4/9 - 3 * t) - 85/3 * np.sin(1/5 - 2 * t) - 1402/7 * np.sin(17/18 - t) + 134/15 * np.sin(4 * t + 23/12) + 17/4 * np.sin(5 * t + 8/9) + 22/5 * np.sin(7 * t + 26/7) + 23/7 * np.sin(8 * t + 4/3) + 7/4 * np.sin(10 * t + 26/7) + 3/2 * np.sin(11 * t + 7/5) - 1034/5) * H(19 * np.pi - t) * H(t - 15 * np.pi) + (-6/5 * np.sin(6/5 - 9 * t) - 4/3 * np.sin(9/7 - 8 * t) - 11/5 * np.sin(3/8 - 7 * t) - 28/5 * np.sin(2/5 - 3 * t) - 31/7 * np.sin(1/7 - 2 * t) + 219/4 * np.sin(t + 25/6) + 9/7 * np.sin(4 * t + 3/7) + 13/5 * np.sin(5 * t + 1/4) + 20/19 * np.sin(6 * t + 5/6) + 5/4 * np.sin(10 * t + 30/7) + 4/7 * np.sin(11 * t + 22/5) + 12/13 * np.sin(12 * t + 30/7) - 683/8) * H(15 * np.pi - t) * H(t - 11 * np.pi) + (-4/7 * np.sin(13/9 - 20 * t) - 13/7 * np.sin(3/5 - 15 * t) - 5/4 * np.sin(11/8 - 5 * t) - np.sin(11/7 - 2 * t) + 11/5 * np.sin(t + 7/2) + 1/2 * np.sin(3 * t + 33/8) + 10/3 * np.sin(4 * t + 5/4) + 6/5 * np.sin(6 * t + 3/5) + 13/3 * np.sin(7 * t + 2/3) + 46/9 * np.sin(8 * t + 1/9) + 23/6 * np.sin(9 * t + 33/7) + 50/9 * np.sin(10 * t + 1/6) + 6/5 * np.sin(11 * t + 8/7) + 18/5 * np.sin(12 * t + 29/10) + 4/5 * np.sin(13 * t + 14/9) + 1/2 * np.sin(14 * t + 11/4) + np.sin(16 * t + 27/7) + 5/6 * np.sin(17 * t + 31/15) + 1/8 * np.sin(18 * t + 5/3) + 1/2 * np.sin(19 * t + 1/9) + 1/3 * np.sin(21 * t + 7/4) + 1/3 * np.sin(22 * t + 11/5) + 2/5 * np.sin(23 * t + 2/3) + 5/11 * np.sin(24 * t + 105/26) + 842/11) * H(11 * np.pi - t) * H(t - 7 * np.pi) + (-399/8 * np.sin(7/6 - t) + 3/4 * np.sin(2 * t + 11/3) + 3/4 * np.sin(3 * t + 1) + 3/5 * np.sin(4 * t + 18/5) + 2/5 * np.sin(5 * t + 22/7) - 2610/7) * H(7 * np.pi - t) * H(t - 3 * np.pi) + (-2/5 * np.sin(5/7 - 4 * t) - 5/4 * np.sin(1 - 2 * t) - 167/6 * np.sin(11/10 - t) + 16/15 * np.sin(3 * t + 39/20) + 3/5 * np.sin(5 * t + 5/2) - 1493/4) * H(3 * np.pi - t) * H(t + np.pi)) * H(np.sqrt(np.sign(np.sin(t/2))))


xdata = []
ydata = []

plt.show()
 
axes = plt.gca()
axes.set_xlim(-500, 500)
axes.set_ylim(-500, 500)
line, = axes.plot(xdata, ydata, 'r-')
t_values = np.arange(0.0, np.pi * 68, 0.1)
for t in t_values:
    xdata.append(x(t))
    ydata.append(y(t))
    line.set_xdata(xdata)
    line.set_ydata(ydata)
    plt.draw()
    plt.pause(1e-17)
    time.sleep(0.1)

plt.show()

Details: https://www.wolframalpha.com/input/?i=likitung+curve

REAL LIKITUNG:

MY LIKITUNG: