Skip to content

Instantly share code, notes, and snippets.

@Shirataki2
Created January 10, 2020 12:50
Show Gist options
  • Save Shirataki2/c25969355bcad908d4efb939508014bc to your computer and use it in GitHub Desktop.
Save Shirataki2/c25969355bcad908d4efb939508014bc to your computer and use it in GitHub Desktop.
import numpy as np
import matplotlib.pyplot as plt
import warnings
warnings.simplefilter('ignore')
t=np.linspace(0,112*np.pi,6000)
sin=np.sin
sqrt=np.sqrt
pi=np.pi
theta=lambda x:np.where(x>=0,1,0)
sgn=lambda x:x/abs(x)
plt.axes().set_aspect('equal')
x=((529/24*sin(1*t+11/7)+611/204*sin(2*t+49/31)+91/29*sin(3*t+49/31)+11/19*sin(4*t+49/30)+152/65*sin(5*t+46/29)+3/29*sin(6*t+76/43)+19/44*sin(7*t+89/19)+8/17*sin(8*t+155/97)+120/239*sin(9*t+35/22)+19/67*sin(10*t+62/39)+41/81*sin(11*t+62/39)+5/42*sin(12*t+57/34)+21459/56)*theta(111*pi-t)*theta(1*t-107*pi)+(-1/18*sin(38/25-12*t)-10/27*sin(61/39-6*t)+1561/33*sin(1*t+107/68)+237/26*sin(2*t+41/26)+1/4*sin(3*t+23/15)+4/31*sin(4*t+46/29)+52/37*sin(5*t+30/19)+4/19*sin(7*t+69/43)+5/27*sin(8*t+23/14)+3/32*sin(9*t+48/31)+1/148*sin(10*t+7/20)+5/23*sin(11*t+11/7)-1457/35)*theta(107*pi-t)*theta(1*t-103*pi)+(-9/40*sin(37/24-12*t)-11/43*sin(57/37-11*t)-24/47*sin(14/9-10*t)-9/7*sin(53/34-6*t)+1217/32*sin(1*t+11/7)+762/103*sin(2*t+41/26)+292/33*sin(3*t+52/33)+34/19*sin(4*t+49/31)+5/13*sin(5*t+25/16)+13/23*sin(7*t+61/39)+14/27*sin(8*t+11/7)+3/7*sin(9*t+67/43)+5/22*sin(13*t+68/43)-3569/8)*theta(103*pi-t)*theta(1*t-99*pi)+(-160/33*sin(31/21-34*t)-117/20*sin(58/37-31*t)-31/34*sin(23/19-29*t)-179/57*sin(112/75-22*t)-51/28*sin(45/32-18*t)-291/91*sin(32/21-14*t)-75/29*sin(58/37-6*t)-104/29*sin(47/30-3*t)-113/20*sin(91/58-2*t)-57/34*sin(45/29-t)+81/17*sin(4*t+74/47)+10/7*sin(5*t+46/29)+14/37*sin(7*t+12/7)+97/18*sin(8*t+80/17)+164/57*sin(9*t+35/22)+66/19*sin(10*t+58/37)+53/6*sin(11*t+41/26)+32/37*sin(12*t+143/31)+143/21*sin(13*t+52/33)+53/30*sin(15*t+58/35)+373/30*sin(16*t+127/27)+391/43*sin(17*t+14/9)+75/41*sin(19*t+56/33)+231/19*sin(20*t+146/31)+239/38*sin(21*t+36/23)+155/156*sin(23*t+71/16)+33/19*sin(24*t+158/35)+120/29*sin(25*t+20/13)+308/27*sin(26*t+14/9)+95/29*sin(27*t+215/46)+45/28*sin(28*t+128/71)+916/131*sin(30*t+39/25)+79/25*sin(32*t+719/154)+6/19*sin(33*t+108/23)+69/20*sin(35*t+43/26)+173/35*sin(36*t+169/36)+97/21*sin(37*t+31/20)+172/19*sin(38*t+63/40)+72/25*sin(39*t+108/23)+158/27*sin(40*t+51/31)+51/47*sin(41*t+33/19)+49/36*sin(42*t+71/41)+178/23*sin(43*t+41/26)+165/43*sin(44*t+71/46)+17527/56)*theta(99*pi-t)*theta(1*t-95*pi)+(-1/8*sin(36/23-22*t)+1201/13*sin(1*t+11/7)+145/34*sin(2*t+36/23)+350/29*sin(3*t+11/7)+83/26*sin(4*t+74/47)+137/18*sin(5*t+11/7)+19/9*sin(6*t+74/47)+23/9*sin(7*t+11/7)+33/43*sin(8*t+27/17)+11/12*sin(9*t+91/58)+49/37*sin(10*t+74/47)+55/63*sin(11*t+41/26)+4/25*sin(12*t+57/37)+43/27*sin(13*t+19/12)+1/15*sin(14*t+113/24)+1/5*sin(15*t+67/43)+2/5*sin(16*t+59/37)+1/98*sin(17*t+30/17)+1/82*sin(18*t+31/7)+15/22*sin(19*t+113/72)+43/85*sin(20*t+63/40)+17/52*sin(21*t+125/78)+5/22*sin(23*t+11/7)+4/37*sin(24*t+36/23)+15/59*sin(25*t+59/38)+4/17*sin(26*t+25/16)+3/13*sin(27*t+20/13)+11/31*sin(28*t+58/37)-13441/56)*theta(95*pi-t)*theta(1*t-91*pi)+(-149/224*sin(42/29-26*t)-4/13*sin(37/26-25*t)-3/20*sin(67/46-24*t)-2/23*sin(4/25-23*t)-18/19*sin(53/35-22*t)-1/22*sin(16/11-20*t)-52/29*sin(299/199-16*t)-31/21*sin(19/13-15*t)-170/37*sin(32/21-14*t)-18/7*sin(40/27-13*t)-94/35*sin(146/97-12*t)-67/29*sin(157/105-11*t)-107/14*sin(65/42-10*t)-34/21*sin(10/7-9*t)-11/36*sin(43/36-7*t)-337/20*sin(31/20-6*t)-166/35*sin(37/24-5*t)-396/41*sin(31/20-4*t)-2259/14*sin(47/30-2*t)+2567/16*sin(1*t+11/7)+370/33*sin(3*t+77/50)+11/30*sin(8*t+60/37)+5/23*sin(17*t+56/41)+4/7*sin(18*t+75/16)+19/23*sin(19*t+45/28)+1/14*sin(21*t+245/53)+5174/19)*theta(91*pi-t)*theta(1*t-87*pi)+(-26/29*sin(65/42-3*t)+97/38*sin(1*t+39/25)+149/34*sin(2*t+61/13)+35/16*sin(4*t+209/45)+293/26*sin(5*t+49/32)+407/42*sin(6*t+14/3)+269/35*sin(7*t+86/57)+133/53*sin(8*t+50/31)+201/26*sin(9*t+43/29)+805/62*sin(10*t+167/36)+205/31*sin(11*t+54/37)+321/38*sin(12*t+83/18)+41/4*sin(13*t+59/40)+11/5*sin(14*t+37/24)+96/55*sin(15*t+60/43)+231/46*sin(16*t+46/31)+175/44*sin(17*t+129/28)+499/73*sin(18*t+25/17)+133/20*sin(19*t+142/31)+53/20*sin(20*t+23/5)+613/28*sin(21*t+197/43)+23/32*sin(22*t+212/141)+47/30*sin(23*t+79/53)+205/57*sin(24*t+41/28)+40/11*sin(25*t+121/26)+43/16*sin(26*t+230/51)+422/29*sin(27*t+25/18)+559/61*sin(28*t+10/7)+31/23*sin(29*t+28/23)+117/34*sin(30*t+31/22)+332/41*sin(31*t+73/16)+133/30*sin(32*t+116/83)+49/26*sin(33*t+159/35)+41/14*sin(34*t+383/87)+74/23*sin(35*t+48/37)+75/23*sin(36*t+71/16)+121/42*sin(37*t+141/31)+50/51*sin(38*t+152/35)+76/35*sin(39*t+518/115)+72/23*sin(40*t+56/39)+95/71*sin(41*t+77/17)+182/31*sin(42*t+151/34)+768/35)*theta(87*pi-t)*theta(1*t-83*pi)+(-4117/73*sin(11/7-2*t)+2006/41*sin(1*t+11/7)+1979/51*sin(3*t+245/52)+994/3)*theta(83*pi-t)*theta(1*t-79*pi)+(-206/35*sin(39/25-5*t)-29/30*sin(41/30-4*t)-1103/27*sin(25/16-3*t)+1802/29*sin(1*t+41/26)+2604/37*sin(2*t+41/26)+62/35*sin(6*t+17/11)+919/38)*theta(79*pi-t)*theta(1*t-75*pi)+(-4319/216*sin(47/30-3*t)+2632/87*sin(1*t+49/31)+4858/55*sin(2*t+30/19)+37/21*sin(4*t+67/46)+16/21*sin(5*t+59/13)+51/13*sin(6*t+69/44)+49/26*sin(7*t+136/29)+52/17*sin(8*t+68/43)+12/19*sin(9*t+121/26)+141/65*sin(10*t+107/67)+7/50*sin(11*t+221/48)+4/9*sin(12*t+85/53)+3886/27)*theta(75*pi-t)*theta(1*t-71*pi)+(-43/38*sin(39/25-10*t)-6/29*sin(143/95-9*t)-64/51*sin(17/11-8*t)-85/24*sin(53/34-6*t)-37/53*sin(142/95-5*t)-1380/29*sin(36/23-2*t)-212/25*sin(39/25-t)+41/28*sin(3*t+47/31)+207/62*sin(4*t+174/37)+3/23*sin(7*t+41/26)+3/25*sin(11*t+37/26)+16/23*sin(12*t+146/31)+7201/21)*theta(71*pi-t)*theta(1*t-67*pi)+(509/51*sin(1*t+146/31)+3728/49*sin(2*t+74/47)+301/62*sin(3*t+68/43)+19/6*sin(4*t+52/33)+7/10*sin(5*t+108/23)+41/9*sin(6*t+30/19)+49/31*sin(7*t+49/31)+22/23*sin(8*t+46/29)+9/25*sin(9*t+41/25)+26/19*sin(10*t+27/17)+10/17*sin(11*t+53/33)+13/34*sin(12*t+71/44)+200/21)*theta(67*pi-t)*theta(1*t-63*pi)+(-41/17*sin(17/11-8*t)-199/25*sin(39/25-7*t)-57/40*sin(54/35-6*t)-487/23*sin(36/23-3*t)+1967/66*sin(1*t+96/61)+19/35*sin(2*t+50/33)+2131/87*sin(4*t+30/19)+41/12*sin(5*t+80/17)+19/14*sin(9*t+221/47)+13/22*sin(10*t+14/9)+15/11*sin(11*t+53/33)+53/23*sin(12*t+35/22)-5558/31)*theta(63*pi-t)*theta(1*t-59*pi)+(-85/37*sin(17/11-54*t)-4/27*sin(37/25-51*t)-22/31*sin(133/89-48*t)-44/27*sin(34/23-37*t)-133/61*sin(31/21-26*t)-2/5*sin(45/29-23*t)-3827/26*sin(69/44-4*t)+179*sin(1*t+113/24)+79162/263*sin(2*t+212/45)+3350/19*sin(3*t+41/26)+3820/57*sin(5*t+49/31)+1830/41*sin(6*t+245/52)+307/31*sin(7*t+245/52)+317/37*sin(8*t+51/31)+2066/81*sin(9*t+113/24)+86/35*sin(10*t+102/23)+502/37*sin(11*t+197/42)+121/91*sin(12*t+45/28)+61/20*sin(13*t+191/41)+99/19*sin(14*t+67/43)+3/8*sin(15*t+43/34)+25/11*sin(16*t+64/41)+71/18*sin(17*t+58/37)+616/137*sin(18*t+383/82)+101/21*sin(19*t+139/30)+291/53*sin(20*t+139/30)+43/13*sin(21*t+117/25)+5/6*sin(22*t+45/23)+533/178*sin(24*t+207/44)+96/23*sin(25*t+60/37)+31/47*sin(27*t+63/34)+7/50*sin(28*t+779/167)+13/27*sin(29*t+41/33)+182/81*sin(30*t+59/37)+15/46*sin(31*t+177/76)+311/44*sin(32*t+145/31)+71/34*sin(33*t+143/32)+31/28*sin(34*t+188/41)+32/21*sin(35*t+46/29)+17/14*sin(36*t+95/21)+83/40*sin(38*t+49/31)+26/31*sin(39*t+39/22)+169/38*sin(40*t+174/37)+46/35*sin(41*t+88/49)+18/31*sin(42*t+83/19)+209/104*sin(43*t+593/127)+52/33*sin(44*t+157/35)+39/35*sin(45*t+141/31)+41/28*sin(46*t+181/39)+44/29*sin(47*t+51/31)+25/36*sin(49*t+40/27)+5/22*sin(50*t+100/23)+19/58*sin(52*t+308/67)+18/19*sin(53*t+52/29)+3469/27)*theta(59*pi-t)*theta(1*t-55*pi)+(-27/35*sin(49/33-38*t)-13/15*sin(71/47-36*t)-41/31*sin(22/15-33*t)-26/17*sin(238/159-29*t)-93/35*sin(32/21-28*t)-87/40*sin(86/57-24*t)-73/85*sin(43/29-22*t)-44/17*sin(20/13-19*t)-179/52*sin(31/20-18*t)-1766/321*sin(23/15-17*t)-6/31*sin(35/33-13*t)-269/29*sin(53/34-10*t)-613/30*sin(67/43-8*t)-578/67*sin(61/39-6*t)-1059/37*sin(25/16-4*t)-1315/39*sin(58/37-2*t)-236/7*sin(80/51-t)+431/76*sin(3*t+52/33)+467/21*sin(5*t+49/31)+41/45*sin(7*t+102/65)+17/3*sin(9*t+30/19)+97/17*sin(11*t+49/31)+13/19*sin(12*t+271/58)+41/7*sin(14*t+45/28)+252/65*sin(15*t+19/12)+15/26*sin(16*t+257/55)+31/21*sin(20*t+77/48)+227/81*sin(21*t+47/29)+1/26*sin(23*t+135/32)+48/19*sin(25*t+31/19)+5/18*sin(26*t+43/26)+2/13*sin(27*t+47/10)+26/17*sin(30*t+147/88)+10/51*sin(31*t+33/25)+51/41*sin(32*t+68/41)+26/47*sin(34*t+569/122)+13/43*sin(35*t+36/23)+23/29*sin(37*t+44/27)+35/39*sin(39*t+77/46)+989/18)*theta(55*pi-t)*theta(1*t-51*pi)+(-68/19*sin(26/17-6*t)-67/19*sin(32/21-4*t)-1177/24*sin(80/51-t)+5591/51*sin(2*t+91/58)+229/28*sin(3*t+31/20)+927/59*sin(5*t+36/23)+148/19*sin(7*t+20/13)+430/29*sin(8*t+87/56)+586/35*sin(9*t+25/16)+31/7*sin(10*t+52/33)+769/62*sin(11*t+136/29)+1630/59*sin(12*t+127/27)+231/47*sin(13*t+44/27)+17/19*sin(14*t+47/27)+209/21*sin(15*t+146/31)+107/42*sin(16*t+80/17)+159/28*sin(17*t+11/7)+184/47*sin(18*t+37/24)+18/53*sin(19*t+35/27)-28410/31)*theta(51*pi-t)*theta(1*t-47*pi)+(-124/41*sin(70/47-16*t)-85/19*sin(17/11-8*t)-133/40*sin(65/42-6*t)-431/26*sin(47/30-4*t)-803/12*sin(80/51-2*t)+1477/27*sin(1*t+212/45)+446/35*sin(3*t+41/26)+273/20*sin(5*t+49/31)+287/30*sin(7*t+46/29)+1615/269*sin(9*t+36/23)+467/34*sin(10*t+80/17)+172/69*sin(11*t+155/33)+684/25*sin(12*t+11/7)+323/39*sin(13*t+25/16)+199/29*sin(14*t+249/53)+221/21*sin(15*t+49/31)+121/23*sin(17*t+55/34)+7/9*sin(18*t+49/33)+59/32*sin(19*t+25/16)+62/103*sin(20*t+44/25)+9157/22)*theta(47*pi-t)*theta(1*t-43*pi)+(298/35*sin(1*t+11/7)+816/43*sin(2*t+11/7)+73/15*sin(3*t+245/52)+21/29*sin(4*t+52/33)+453/13)*theta(43*pi-t)*theta(1*t-39*pi)+(-103/69*sin(91/58-5*t)-13/17*sin(41/27-4*t)-299/33*sin(36/23-3*t)+118/13*sin(1*t+52/33)+510/73*sin(2*t+107/68)+1751/27)*theta(39*pi-t)*theta(1*t-35*pi)+(-229/24*sin(80/51-3*t)+533/38*sin(1*t+11/7)+85/22*sin(2*t+11/7)+1/146*sin(4*t+30/29)-22247/108)*theta(35*pi-t)*theta(1*t-31*pi)+(-847/53*sin(47/30-3*t)+194/27*sin(1*t+11/7)+5/6*sin(2*t+14/9)-5001/37)*theta(31*pi-t)*theta(1*t-27*pi)+(-82/31*sin(212/141-16*t)-839/140*sin(32/21-14*t)-138/19*sin(63/41-13*t)-98/29*sin(56/37-12*t)-272/23*sin(71/46-11*t)-126/25*sin(17/11-9*t)-103/20*sin(47/30-6*t)+2768/17*sin(1*t+11/7)+40/27*sin(2*t+28/19)+478/21*sin(3*t+245/52)+521/32*sin(4*t+113/24)+254/25*sin(5*t+46/29)+290/31*sin(7*t+30/19)+148/45*sin(8*t+47/30)+148/23*sin(10*t+69/44)+7/34*sin(15*t+17/29)+124/71*sin(17*t+49/31)+26/23*sin(18*t+45/28)+46/35*sin(19*t+28/17)+940/21)*theta(27*pi-t)*theta(1*t-23*pi)+(-94/33*sin(44/29-27*t)-93/32*sin(35/23-24*t)-669/55*sin(47/30-23*t)-38/3*sin(47/30-13*t)-777/26*sin(31/20-4*t)+7013/20*sin(1*t+80/51)+1121/20*sin(2*t+36/23)+8619/44*sin(3*t+113/24)+3055/32*sin(5*t+80/17)+1212/17*sin(6*t+25/16)+630/31*sin(7*t+39/25)+1813/24*sin(8*t+197/42)+1681/35*sin(9*t+48/31)+435/28*sin(10*t+68/43)+1001/37*sin(11*t+221/47)+103/38*sin(12*t+117/25)+746/41*sin(14*t+17/11)+529/29*sin(15*t+80/17)+139/21*sin(16*t+116/25)+47/33*sin(17*t+13/11)+131/30*sin(18*t+97/21)+189/32*sin(19*t+117/25)+583/19*sin(20*t+39/25)+45/7*sin(21*t+373/80)+115/18*sin(22*t+107/68)+767/43*sin(25*t+155/33)+511/93*sin(26*t+23/15)+29/16*sin(28*t+265/59)+106/39*sin(29*t+67/46)+67/21*sin(30*t+31/19)+94/37*sin(31*t+61/13)+32/33*sin(32*t+63/37)+90/29*sin(33*t+17/12)+266/33*sin(34*t+51/11)+691/148*sin(35*t+3/2)+225/44*sin(36*t+53/35)+125/21*sin(37*t+134/29)+51/29*sin(38*t+25/18)+50/37*sin(39*t+17/11)+59/33*sin(40*t+155/33)+132/61*sin(41*t+78/17)+72/25*sin(42*t+43/29)+46/29*sin(43*t+134/29)+112/39*sin(44*t+173/37)+44/39*sin(45*t+45/32)-9397/55)*theta(23*pi-t)*theta(1*t-19*pi)+(-152/51*sin(178/119-47*t)-83/17*sin(20/13-37*t)-107/21*sin(28/19-35*t)-159/38*sin(42/29-34*t)-145/38*sin(38/25-30*t)-124/23*sin(47/32-29*t)-113/14*sin(51/35-25*t)-309/25*sin(163/109-24*t)-633/32*sin(45/29-21*t)-421/23*sin(65/42-11*t)-85/21*sin(24/17-10*t)-5200/107*sin(59/38-3*t)-6214/31*sin(36/23-2*t)+3868/19*sin(1*t+11/7)+2002/37*sin(4*t+46/29)+1720/19*sin(5*t+43/27)+578/5*sin(6*t+59/37)+247/4*sin(7*t+45/28)+2990/67*sin(8*t+51/32)+11/3*sin(9*t+75/16)+151/14*sin(12*t+36/23)+52/37*sin(13*t+58/13)+306/19*sin(14*t+50/31)+18/37*sin(15*t+79/21)+247/31*sin(16*t+65/41)+28/29*sin(17*t+117/59)+119/25*sin(18*t+41/27)+233/24*sin(19*t+136/29)+653/27*sin(20*t+46/29)+47/20*sin(22*t+68/47)+76/29*sin(23*t+63/41)+63/37*sin(26*t+59/32)+27/2*sin(27*t+18/11)+3841/960*sin(28*t+7/4)+177/32*sin(31*t+59/36)+58/11*sin(32*t+55/32)+94/43*sin(33*t+37/22)+127/42*sin(36*t+64/41)+121/31*sin(38*t+43/26)+192/47*sin(39*t+29/17)+47/35*sin(40*t+136/29)+9/14*sin(41*t+73/43)+191/49*sin(42*t+53/31)+95/61*sin(43*t+79/45)+69/55*sin(44*t+49/26)+56/23*sin(45*t+49/29)+17/15*sin(46*t+43/27)+13/42*sin(48*t+22/21)+14/33*sin(49*t+92/21)-8528/25)*theta(19*pi-t)*theta(1*t-15*pi)+(-101/50*sin(18/13-24*t)-25/27*sin(41/32-23*t)-403/48*sin(23/16-19*t)-267/44*sin(55/38-18*t)-722/73*sin(16/11-15*t)-452/35*sin(63/44-14*t)-133/26*sin(60/43-13*t)-1052/59*sin(25/17-10*t)-1128/31*sin(3/2-9*t)-713/48*sin(80/53-8*t)-779/44*sin(26/17-6*t)-391/10*sin(26/17-5*t)-89/5*sin(29/19-4*t)-1109/10*sin(25/16-2*t)+1573/24*sin(1*t+61/39)+223/34*sin(3*t+229/153)+809/35*sin(7*t+35/22)+164/23*sin(11*t+113/72)+17/10*sin(12*t+23/5)+171/86*sin(16*t+43/24)+48/5*sin(17*t+51/31)+62/39*sin(20*t+39/25)+37/12*sin(21*t+75/43)+25/29*sin(22*t+47/26)+7/13*sin(25*t+43/31)+9/8*sin(26*t+202/43)-438/17)*theta(15*pi-t)*theta(1*t-11*pi)+(-13/11*sin(317/211-23*t)-151/46*sin(17/11-20*t)-70/39*sin(32/21-17*t)-12/25*sin(43/31-16*t)-311/32*sin(45/29-11*t)-127/21*sin(31/20-9*t)-1813/47*sin(25/16-6*t)-1729/30*sin(61/39-4*t)+100/31*sin(1*t+67/43)+4206/35*sin(2*t+74/47)+1352/25*sin(3*t+63/40)+232/15*sin(5*t+30/19)+223/11*sin(7*t+49/31)+500/19*sin(8*t+68/43)+7/2*sin(10*t+102/65)+141/37*sin(12*t+113/24)+406/81*sin(13*t+67/42)+274/61*sin(14*t+43/27)+31/21*sin(15*t+46/29)+67/33*sin(18*t+80/17)+77/20*sin(19*t+45/28)+56/33*sin(21*t+43/27)+32/19*sin(22*t+11/7)+236/79*sin(24*t+29/18)+21/26*sin(25*t+8/5)-5339/12)*theta(11*pi-t)*theta(1*t-7*pi)+(-433/63*sin(71/46-20*t)-509/40*sin(39/25-15*t)-109/15*sin(17/11-13*t)-443/29*sin(58/37-6*t)-2283/47*sin(58/37-5*t)-1091/11*sin(36/23-2*t)+8067/109*sin(1*t+69/44)+18461/142*sin(3*t+63/40)+3241/46*sin(4*t+41/26)+1471/64*sin(7*t+65/41)+19/9*sin(8*t+91/57)+754/33*sin(9*t+49/31)+562/25*sin(10*t+85/54)+643/29*sin(11*t+52/33)+237/26*sin(12*t+30/19)+180/17*sin(14*t+43/27)+145/38*sin(16*t+53/32)+86/25*sin(17*t+207/44)+11/36*sin(18*t+9/2)+70/31*sin(19*t+113/24)+429/86*sin(21*t+101/63)+7/4*sin(22*t+48/31)+51/8*sin(23*t+91/57)+176/37*sin(24*t+43/27)+8021/38)*theta(7*pi-t)*theta(1*t-3*pi)+(-92/29*sin(47/30-21*t)-131/21*sin(59/38-19*t)-177/23*sin(53/34-17*t)-1580/69*sin(87/56-9*t)-1348/23*sin(53/34-8*t)-524/15*sin(69/44-5*t)-4208/35*sin(47/30-2*t)+767/21*sin(1*t+107/68)+1503/38*sin(3*t+11/7)+571/33*sin(4*t+212/45)+947/26*sin(6*t+46/29)+101/33*sin(7*t+523/112)+7/2*sin(10*t+75/16)+331/25*sin(11*t+8/5)+35/32*sin(12*t+329/72)+195/34*sin(13*t+39/25)+89/16*sin(14*t+81/52)+43/23*sin(15*t+21/13)+169/19*sin(16*t+27/17)+12/19*sin(18*t+26/17)+47/20*sin(20*t+68/43)+27/29*sin(22*t+46/27)-12713/48)*theta(3*pi-t)*theta(1*t+pi))*theta(sqrt(sgn(sin(1*t/2))))
y=((-83/13*sin(61/39-5*t)-253/59*sin(36/23-3*t)+938/75*sin(1*t+11/7)+17/16*sin(2*t+31/20)+1/23*sin(4*t+23/28)+47/20*sin(6*t+11/7)+22/15*sin(7*t+46/29)+1/24*sin(8*t+278/61)+1/4*sin(9*t+44/27)+17/18*sin(10*t+77/48)+12/23*sin(11*t+55/34)+1/9*sin(12*t+152/33)+2053/20)*theta(111*pi-t)*theta(1*t-107*pi)+(-1383/79*sin(11/7-t)+114/5*sin(2*t+107/68)+15/19*sin(3*t+14/3)+44/13*sin(4*t+107/68)+3/5*sin(5*t+46/29)+70/33*sin(6*t+65/41)+1/32*sin(7*t+23/12)+31/29*sin(8*t+52/33)+1/11*sin(9*t+169/37)+14/23*sin(10*t+30/19)+1/53*sin(11*t+109/73)+3/7*sin(12*t+35/22)-20442/13)*theta(107*pi-t)*theta(1*t-103*pi)+(-2799/55*sin(113/72-t)+409/28*sin(2*t+11/7)+22/29*sin(3*t+65/41)+27/7*sin(4*t+11/7)+21/20*sin(5*t+179/38)+91/53*sin(6*t+49/31)+21/26*sin(7*t+41/26)+543/272*sin(8*t+96/61)+4/13*sin(9*t+51/32)+40/67*sin(10*t+69/44)+1/105*sin(11*t+41/24)+30/37*sin(12*t+65/41)+31/72*sin(13*t+49/31)-57768/37)*theta(103*pi-t)*theta(1*t-99*pi)+(-29/27*sin(54/37-40*t)-23/30*sin(21/20-39*t)-503/42*sin(25/16-37*t)-49/25*sin(71/53-35*t)-115/22*sin(92/61-34*t)-225/38*sin(25/16-32*t)-188/23*sin(29/19-23*t)-456/37*sin(36/23-21*t)-205/32*sin(47/30-10*t)-32/5*sin(64/41-8*t)-107/32*sin(53/34-6*t)-63/43*sin(57/37-2*t)+223/28*sin(1*t+25/16)+273/38*sin(3*t+63/40)+616/137*sin(4*t+108/23)+64/11*sin(5*t+35/22)+487/55*sin(7*t+47/30)+386/73*sin(9*t+52/33)+44/15*sin(11*t+67/42)+74/19*sin(12*t+207/44)+179/52*sin(13*t+43/27)+110/47*sin(14*t+49/32)+394/51*sin(15*t+65/42)+74/21*sin(16*t+145/31)+247/35*sin(17*t+173/37)+125/27*sin(18*t+43/28)+9/38*sin(19*t+11/67)+49/19*sin(20*t+21/13)+90/13*sin(22*t+25/16)+71/57*sin(24*t+119/27)+133/11*sin(25*t+136/29)+81/8*sin(26*t+45/29)+63/31*sin(27*t+18/13)+54/31*sin(28*t+79/42)+97/34*sin(29*t+144/31)+379/52*sin(30*t+65/42)+2099/350*sin(31*t+25/16)+13/17*sin(33*t+33/16)+49/30*sin(36*t+163/36)+359/45*sin(38*t+95/61)+227/81*sin(41*t+23/14)+163/28*sin(42*t+135/29)+14/17*sin(43*t+100/47)+568/53*sin(44*t+20/13)+5111/21)*theta(99*pi-t)*theta(1*t-95*pi)+(-5/9*sin(63/41-26*t)-4/31*sin(48/31-25*t)-8/41*sin(77/53-23*t)-1/4*sin(20/13-15*t)-133/23*sin(47/30-12*t)-605/56*sin(64/41-11*t)-97/31*sin(48/31-8*t)-423/31*sin(36/23-7*t)-763/41*sin(58/37-5*t)-1177/43*sin(11/7-t)+31/28*sin(2*t+67/43)+287/32*sin(3*t+113/24)+915/52*sin(4*t+11/7)+3231/202*sin(6*t+52/33)+647/41*sin(9*t+49/31)+123/52*sin(10*t+74/47)+73/19*sin(13*t+212/45)+28/19*sin(14*t+19/12)+89/21*sin(16*t+19/12)+193/58*sin(17*t+41/26)+83/25*sin(18*t+11/7)+87/41*sin(19*t+19/12)+3/19*sin(20*t+17/12)+44/31*sin(21*t+43/27)+22/9*sin(22*t+46/29)+21/13*sin(24*t+19/12)+17/44*sin(27*t+51/32)+18/17*sin(28*t+27/17)+52/7)*theta(95*pi-t)*theta(1*t-91*pi)+(-124/29*sin(23/16-26*t)-131/29*sin(37/25-25*t)-76/35*sin(69/44-15*t)-73/11*sin(26/17-14*t)-106/31*sin(77/51-13*t)-29/19*sin(48/31-11*t)-50/11*sin(29/19-10*t)-308/31*sin(95/61-9*t)-1/23*sin(25/24-6*t)-1017/46*sin(36/23-5*t)-337/18*sin(61/39-3*t)-257/15*sin(25/16-2*t)-522/31*sin(47/30-t)+484/19*sin(4*t+85/54)+120/19*sin(7*t+8/5)+483/40*sin(8*t+65/41)+79/14*sin(12*t+36/23)+427/67*sin(16*t+34/21)+207/47*sin(17*t+43/26)+39/19*sin(18*t+62/37)+22/35*sin(19*t+88/19)+2/33*sin(20*t+64/33)+23/30*sin(21*t+59/33)+167/38*sin(22*t+48/29)+81/22*sin(23*t+227/136)+25/24*sin(24*t+71/44)-2798/49)*theta(91*pi-t)*theta(1*t-87*pi)+(-9/46*sin(5/12-35*t)+433/63*sin(1*t+81/52)+43/16*sin(2*t+218/47)+40/9*sin(3*t+58/39)+158/27*sin(4*t+172/37)+431/37*sin(5*t+35/23)+91/27*sin(6*t+377/82)+100/27*sin(7*t+45/31)+197/34*sin(8*t+64/41)+39/14*sin(9*t+47/32)+92/17*sin(10*t+43/28)+233/43*sin(11*t+135/29)+194/33*sin(12*t+55/37)+409/26*sin(13*t+125/27)+346/23*sin(14*t+49/33)+253/27*sin(15*t+97/21)+211/17*sin(16*t+22/15)+358/35*sin(17*t+134/29)+469/46*sin(18*t+25/17)+92/39*sin(19*t+24/17)+147/13*sin(20*t+47/32)+259/72*sin(21*t+299/64)+381/29*sin(22*t+214/47)+152/45*sin(23*t+173/37)+103/30*sin(24*t+361/82)+15/28*sin(25*t+27/32)+335/32*sin(26*t+86/19)+321/37*sin(27*t+55/12)+118/15*sin(28*t+10/7)+77/29*sin(29*t+39/29)+57/8*sin(30*t+47/33)+29/13*sin(31*t+192/41)+207/59*sin(32*t+112/25)+513/79*sin(33*t+29/22)+89/23*sin(34*t+71/16)+145/36*sin(36*t+38/27)+122/47*sin(37*t+41/9)+1/4*sin(38*t+93/23)+139/50*sin(39*t+279/61)+35/52*sin(40*t+95/23)+116/25*sin(41*t+25/19)+151/76*sin(42*t+135/31)+14846/61)*theta(87*pi-t)*theta(1*t-83*pi)+(-911/27*sin(80/51-3*t)+793/16*sin(1*t+113/72)+1600/23*sin(2*t+11/7)+15991/69)*theta(83*pi-t)*theta(1*t-79*pi)+(-974/19*sin(80/51-t)+3857/46*sin(2*t+30/19)+426/13*sin(3*t+27/17)+82/39*sin(4*t+5/3)+68/9*sin(5*t+61/38)+27/37*sin(6*t+57/32)+8477/37)*theta(79*pi-t)*theta(1*t-75*pi)+(-3/25*sin(23/15-11*t)-47/43*sin(32/21-8*t)-851/21*sin(36/23-t)+781/26*sin(2*t+49/31)+775/22*sin(3*t+46/29)+23/4*sin(4*t+27/17)+10/13*sin(5*t+61/39)+607/91*sin(6*t+123/77)+47/38*sin(7*t+53/33)+110/47*sin(9*t+66/41)+41/38*sin(10*t+70/43)+14/11*sin(12*t+21/13)-81/16)*theta(75*pi-t)*theta(1*t-71*pi)+(-17/37*sin(43/28-12*t)-52/67*sin(43/28-11*t)-13/10*sin(45/29-8*t)-31/18*sin(37/24-7*t)-306/31*sin(39/25-4*t)-43/28*sin(48/31-3*t)+91/10*sin(1*t+74/47)+454/43*sin(2*t+30/19)+69/19*sin(5*t+49/31)+91/73*sin(6*t+61/38)+39/37*sin(9*t+8/5)+13/41*sin(10*t+133/83)+59329/144)*theta(71*pi-t)*theta(1*t-67*pi)+(-30/17*sin(61/39-12*t)-15/29*sin(43/28-9*t)-239/50*sin(25/16-8*t)-595/31*sin(36/23-4*t)+49/99*sin(1*t+47/10)+523/49*sin(2*t+63/40)+93/46*sin(3*t+36/23)+31/23*sin(5*t+58/37)+13/11*sin(6*t+27/17)+33/17*sin(7*t+102/65)+2/9*sin(10*t+117/73)+35/27*sin(11*t+11/7)+3599/9)*theta(67*pi-t)*theta(1*t-63*pi)+(-437/131*sin(25/16-10*t)-269/43*sin(59/38-8*t)-151/49*sin(45/29-4*t)-451/38*sin(39/25-3*t)-509/30*sin(11/7-t)+173/3*sin(2*t+52/33)+1214/47*sin(5*t+30/19)+397/30*sin(6*t+65/41)+71/9*sin(7*t+49/31)+8/29*sin(9*t+69/47)+39/43*sin(11*t+49/31)+7/23*sin(12*t+93/56)+4273/38)*theta(63*pi-t)*theta(1*t-59*pi)+(-21/29*sin(44/45-36*t)-433/21*sin(53/34-30*t)-587/36*sin(47/31-8*t)-3175/26*sin(25/16-6*t)-2589/32*sin(11/7-t)+34/31*sin(2*t+125/38)+3176/19*sin(3*t+179/38)+223/52*sin(4*t+77/46)+1655/16*sin(5*t+11/7)+9074/75*sin(7*t+63/40)+158/11*sin(9*t+63/40)+874/27*sin(10*t+155/33)+265/28*sin(11*t+93/20)+756/25*sin(12*t+289/62)+474/31*sin(13*t+65/14)+744/41*sin(14*t+229/49)+74/13*sin(15*t+133/29)+1576/59*sin(16*t+113/24)+377/41*sin(17*t+46/29)+206/15*sin(18*t+52/33)+311/29*sin(19*t+26/17)+53/88*sin(20*t+13/28)+491/70*sin(21*t+10/7)+464/33*sin(22*t+49/33)+872/39*sin(23*t+45/29)+20/27*sin(24*t+49/23)+765/52*sin(25*t+129/28)+1493/39*sin(26*t+135/29)+817/33*sin(27*t+261/56)+803/42*sin(28*t+19/12)+725/41*sin(29*t+102/65)+31/10*sin(31*t+302/67)+37/18*sin(32*t+67/33)+1312/99*sin(33*t+51/11)+230/31*sin(34*t+68/15)+13/25*sin(35*t+43/15)+479/49*sin(37*t+69/44)+181/37*sin(38*t+75/16)+520/57*sin(39*t+131/28)+302/55*sin(40*t+26/15)+39/47*sin(41*t+539/135)+1232/93*sin(42*t+191/41)+33/35*sin(43*t+40/11)+22/17*sin(44*t+35/8)+377/60*sin(45*t+152/33)+69/22*sin(46*t+179/40)+67/42*sin(47*t+30/19)+7/31*sin(48*t+77/29)+33/31*sin(49*t+93/20)+41/31*sin(50*t+174/37)+21/29*sin(51*t+317/88)+868/137*sin(52*t+129/28)+71/19*sin(53*t+159/35)+13/20*sin(54*t+70/27)+11225/33)*theta(59*pi-t)*theta(1*t-55*pi)+(-25/18*sin(26/17-39*t)-5/26*sin(40/29-38*t)-4/23*sin(61/43-36*t)-54/163*sin(33/23-31*t)-1/66*sin(19/24-27*t)-16/11*sin(83/55-26*t)-100/27*sin(55/36-25*t)-97/49*sin(181/121-24*t)-54/29*sin(340/227-20*t)-13/17*sin(43/28-19*t)-680/47*sin(31/20-13*t)-329/66*sin(14/9-11*t)-3254/35*sin(36/23-3*t)+106/17*sin(1*t+46/29)+1372/13*sin(2*t+63/40)+645/52*sin(4*t+41/26)+2806/37*sin(5*t+49/31)+1066/73*sin(6*t+52/33)+109/9*sin(7*t+41/26)+1043/48*sin(8*t+19/12)+252/31*sin(9*t+101/63)+67/23*sin(10*t+41/26)+103/33*sin(12*t+21/13)+111/32*sin(14*t+59/37)+72/23*sin(15*t+43/27)+217/34*sin(16*t+51/32)+35/24*sin(17*t+49/30)+45/17*sin(18*t+19/12)+1589/318*sin(21*t+47/29)+391/78*sin(22*t+37/23)+7/30*sin(23*t+108/65)+59/12*sin(28*t+21/13)+9/20*sin(29*t+60/37)+1/12*sin(30*t+47/27)+23/32*sin(32*t+95/61)+5/36*sin(33*t+71/39)+49/50*sin(34*t+37/23)+43/36*sin(35*t+178/107)+4/45*sin(37*t+205/44)-7753/33)*theta(55*pi-t)*theta(1*t-51*pi)+(-69/28*sin(73/47-14*t)-103/90*sin(116/83-13*t)-427/37*sin(39/25-9*t)-125/12*sin(47/30-7*t)-56/41*sin(374/249-5*t)+313/19*sin(1*t+146/31)+1543/23*sin(2*t+179/38)+289/15*sin(3*t+11/7)+52/7*sin(4*t+41/26)+495/26*sin(6*t+91/58)+491/29*sin(8*t+80/51)+853/46*sin(10*t+61/39)+371/74*sin(11*t+83/55)+1137/62*sin(12*t+36/23)+202/27*sin(15*t+47/30)+42/13*sin(16*t+57/37)+26/53*sin(17*t+8/5)+56/55*sin(18*t+133/89)+9/32*sin(19*t+164/37)-46465/88)*theta(51*pi-t)*theta(1*t-47*pi)+(-38/15*sin(59/38-18*t)-208/29*sin(59/38-16*t)-73/39*sin(64/41-14*t)-13/59*sin(37/30-8*t)-297/32*sin(47/30-6*t)+2519/28*sin(1*t+245/52)+83/15*sin(2*t+179/38)+603/31*sin(3*t+69/44)+128/27*sin(4*t+52/33)+41/21*sin(5*t+61/13)+53/44*sin(7*t+53/33)+215/29*sin(9*t+63/40)+37/26*sin(10*t+37/23)+419/37*sin(11*t+41/26)+489/31*sin(12*t+74/47)+177/23*sin(13*t+85/54)+43/18*sin(15*t+35/22)+14/17*sin(17*t+37/22)+16/25*sin(19*t+34/21)+35/34*sin(20*t+37/23)-23041/45)*theta(47*pi-t)*theta(1*t-43*pi)+(389/32*sin(1*t+212/45)+674/55*sin(2*t+11/7)+279/29*sin(3*t+11/7)+41/30*sin(4*t+245/52)-29349/38)*theta(43*pi-t)*theta(1*t-39*pi)+(33/5*sin(1*t+245/52)+379/22*sin(2*t+52/33)+115/32*sin(3*t+27/17)+1/7*sin(4*t+29/19)+8/25*sin(5*t+35/22)-21578/35)*theta(39*pi-t)*theta(1*t-35*pi)+(-47/39*sin(39/25-4*t)+14/43*sin(1*t+69/44)+444/19*sin(2*t+11/7)+19/21*sin(3*t+11/7)-11681/15)*theta(35*pi-t)*theta(1*t-31*pi)+(5/14*sin(1*t+117/25)+715/38*sin(2*t+85/54)+8/37*sin(3*t+145/31)-10395/17)*theta(31*pi-t)*theta(1*t-27*pi)+(-233/51*sin(113/72-9*t)-563/18*sin(36/23-t)+2883/16*sin(2*t+107/68)+321/20*sin(3*t+46/29)+154/17*sin(4*t+14/3)+53/18*sin(5*t+114/25)+593/30*sin(6*t+46/29)+46/17*sin(7*t+47/29)+94/15*sin(8*t+19/12)+193/19*sin(10*t+85/54)+45/44*sin(11*t+257/55)+23/9*sin(12*t+53/33)+49/19*sin(13*t+89/19)+79/21*sin(14*t+101/63)+65/54*sin(15*t+50/29)+75/113*sin(16*t+41/25)+26/51*sin(17*t+19/11)+70/47*sin(18*t+8/5)+25/14*sin(19*t+57/35)-4289/14)*theta(27*pi-t)*theta(1*t-23*pi)+(-24/31*sin(46/35-34*t)-163/17*sin(47/30-26*t)-123/44*sin(57/37-21*t)-101/12*sin(14/9-19*t)-8/21*sin(23/15-17*t)-779/25*sin(36/23-4*t)+689/38*sin(1*t+145/31)+2649/10*sin(2*t+80/51)+113/6*sin(3*t+69/44)+1707/13*sin(5*t+113/24)+281/16*sin(6*t+25/16)+1817/34*sin(7*t+45/29)+453/52*sin(8*t+65/14)+55/13*sin(9*t+41/27)+1927/29*sin(10*t+25/16)+1073/31*sin(11*t+69/44)+2221/33*sin(12*t+36/23)+287/34*sin(13*t+53/33)+9639/482*sin(14*t+61/39)+99/13*sin(15*t+86/57)+41/40*sin(16*t+40/23)+222/25*sin(18*t+63/41)+434/27*sin(20*t+108/23)+17/32*sin(22*t+149/42)+791/34*sin(23*t+65/42)+19/4*sin(24*t+143/31)+63/41*sin(25*t+45/31)+185/27*sin(27*t+131/28)+46/33*sin(28*t+76/59)+149/34*sin(29*t+46/29)+208/43*sin(30*t+353/235)+76/33*sin(31*t+32/19)+87/17*sin(32*t+40/27)+73/49*sin(33*t+203/45)+77/47*sin(35*t+84/19)+283/29*sin(36*t+61/40)+88/19*sin(37*t+53/33)+71/28*sin(38*t+20/13)+35/18*sin(39*t+10/7)+54/13*sin(40*t+181/39)+24/13*sin(41*t+34/23)+39/25*sin(42*t+117/70)+23/26*sin(43*t+32/25)+212/37*sin(44*t+125/27)+21/20*sin(45*t+37/29)-3585/8)*theta(23*pi-t)*theta(1*t-19*pi)+(-71/85*sin(14/9-48*t)-17/29*sin(10/7-44*t)-150/37*sin(19/13-43*t)-51/8*sin(43/28-39*t)-296/31*sin(76/51-37*t)-147/32*sin(33/23-36*t)-8/33*sin(12/31-33*t)-127/20*sin(53/37-32*t)-573/38*sin(47/32-31*t)-25/36*sin(35/26-30*t)-774/31*sin(26/17-26*t)-52/17*sin(13/10-23*t)-811/26*sin(41/27-22*t)-3257/71*sin(45/29-19*t)-921/47*sin(65/42-9*t)-522/31*sin(26/17-8*t)-861/22*sin(17/11-7*t)-897/35*sin(66/43-6*t)-677/15*sin(53/34-5*t)-33/40*sin(37/34-4*t)+433/18*sin(1*t+51/32)+8978/57*sin(2*t+52/33)+839/18*sin(3*t+146/31)+190/23*sin(10*t+81/52)+17/12*sin(11*t+59/13)+401/20*sin(12*t+50/31)+461/31*sin(13*t+59/36)+266/11*sin(14*t+34/21)+222/11*sin(15*t+31/19)+183/46*sin(16*t+58/37)+241/22*sin(17*t+103/22)+865/32*sin(18*t+25/16)+699/32*sin(20*t+49/31)+306/41*sin(21*t+49/30)+255/31*sin(24*t+63/37)+227/38*sin(25*t+31/21)+60/31*sin(27*t+24/19)+392/45*sin(28*t+77/45)+390/31*sin(29*t+107/64)+28/25*sin(34*t+53/27)+513/110*sin(35*t+51/31)+18/7*sin(38*t+34/25)+127/25*sin(40*t+67/41)+31/9*sin(41*t+76/41)+61/28*sin(42*t+53/33)+52/31*sin(45*t+33/19)+5/16*sin(46*t+134/73)+33/83*sin(47*t+39/29)+19/27*sin(49*t+55/36)-23029/24)*theta(19*pi-t)*theta(1*t-15*pi)+(-7/4*sin(42/29-23*t)-7/10*sin(57/49-22*t)-188/39*sin(55/37-17*t)-59/25*sin(31/23-16*t)-52/21*sin(39/28-15*t)-543/52*sin(49/32-13*t)-243/20*sin(19/13-11*t)-1019/58*sin(142/95-10*t)-1465/157*sin(67/43-9*t)-1276/39*sin(23/15-7*t)-1517/40*sin(49/32-6*t)+5219/33*sin(1*t+35/22)+15493/36*sin(2*t+49/31)+1583/21*sin(3*t+103/22)+656/7*sin(4*t+46/29)+524/105*sin(5*t+83/51)+2064/59*sin(8*t+65/41)+39/46*sin(12*t+28/65)+54/23*sin(14*t+44/29)+sin(18*t+40/27)+9/22*sin(19*t+137/31)+269/55*sin(20*t+47/28)+33/83*sin(21*t+146/43)+103/24*sin(24*t+68/39)+31/13*sin(25*t+100/53)+1/18*sin(26*t+4/29)-67025/39)*theta(15*pi-t)*theta(1*t-11*pi)+(-103/23*sin(71/46-25*t)-127/19*sin(54/35-20*t)-79/28*sin(72/47-19*t)-79/27*sin(38/25-18*t)-158/35*sin(43/28-15*t)-236/37*sin(71/46-14*t)-229/33*sin(66/43-12*t)-1322/33*sin(14/9-9*t)-263/7*sin(61/39-7*t)-447/25*sin(58/37-5*t)-8509/39*sin(47/30-3*t)-11816/41*sin(58/37-2*t)+13897/38*sin(1*t+85/54)+399/11*sin(4*t+179/38)+737/36*sin(6*t+30/19)+1248/79*sin(8*t+46/29)+205/12*sin(10*t+59/37)+273/38*sin(11*t+212/45)+269/55*sin(13*t+146/31)+143/18*sin(16*t+83/52)+87/23*sin(17*t+193/41)+95/46*sin(21*t+61/37)+158/105*sin(22*t+37/24)+81/43*sin(23*t+108/23)+25/41*sin(24*t+74/45)-64516/39)*theta(11*pi-t)*theta(1*t-7*pi)+(-71/21*sin(54/35-23*t)-55/21*sin(38/25-18*t)-251/42*sin(36/23-17*t)-71/22*sin(39/25-16*t)-120/23*sin(48/31-13*t)-67/23*sin(3/2-11*t)-354/41*sin(67/43-9*t)-497/19*sin(39/25-7*t)-5852/97*sin(58/37-3*t)-593/29*sin(69/44-2*t)+3680/49*sin(1*t+11/7)+155/8*sin(4*t+30/19)+181/32*sin(5*t+46/29)+229/21*sin(6*t+179/38)+1101/55*sin(8*t+49/31)+398/19*sin(10*t+46/29)+121/23*sin(12*t+39/25)+783/94*sin(14*t+123/77)+74/47*sin(15*t+49/29)+266/43*sin(19*t+93/58)+164/33*sin(20*t+19/12)+56/25*sin(21*t+141/88)+57/25*sin(22*t+14/9)+1/23*sin(24*t+26/11)-39578/17)*theta(7*pi-t)*theta(1*t-3*pi)+(-155/43*sin(45/29-21*t)-139/36*sin(102/65-19*t)-112/17*sin(71/46-16*t)-267/22*sin(45/29-12*t)-1556/63*sin(14/9-11*t)-872/39*sin(69/44-9*t)-143/12*sin(25/16-5*t)-908/27*sin(47/30-3*t)+661/7*sin(1*t+96/61)+327/17*sin(2*t+74/47)+1106/37*sin(4*t+96/61)+58/15*sin(6*t+52/33)+41/35*sin(7*t+139/30)+487/73*sin(8*t+91/57)+5/26*sin(10*t+12/7)+143/25*sin(13*t+127/27)+255/29*sin(14*t+29/18)+444/71*sin(15*t+45/28)+27/25*sin(17*t+3/2)+113/42*sin(18*t+30/19)+275/42*sin(20*t+67/42)+79/38*sin(22*t+48/31)-48856/21)*theta(3*pi-t)*theta(1*t+pi))*theta(sqrt(sgn(sin(1*t/2))))
x[x==0]=np.nan
y[y==0]=np.nan
plt.plot(x,y)
plt.savefig("/images/hoge.png")
@eggplants
Copy link

eggplants commented Jan 10, 2020

Wolframの公式変換です

echo "${wolfram_fomula}" > f
sed -r "s_π_pi_g;s_θ_theta_g;s_\)\(_)*(_g;s_(x|y)\(t\)_\1_g;s_([0-9]|\)) (theta|sqrt|sin|sgn|pi|t)_\1*\2_g;s_ __g" f > ff

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment