Author: Piotr A. Zolnierczuk (zolnierczukp at ornl dot gov)
Drawing an elephant with four complex parameters
Jurgen Mayer, Khaled Khairy, and Jonathon Howard,
Am. J. Phys. 78, 648 (2010), DOI:10.1119/1.3254017
import matplotlib.pyplot as pylab
p1, p2, p3, p4 = (50 - 30j, 18 + 8j, 12 - 10j, -14 - 60j )
f = f + A[k]*np.cos(k*t) + B[k]*np.sin(k*t)
def elephant(t, p1, p2, p3, p4, p5):
Cx = np.zeros((npar,), dtype='complex')
Cy = np.zeros((npar,), dtype='complex')
Cy[1] = p4.imag + p1.imag*1j
x = np.append(fourier(t,Cx), [-p5.imag])
y = np.append(fourier(t,Cy), [p5.imag])
x, y = elephant(np.linspace(0,2*np.pi,1000), p1, p2, p3, p4, p5)
pylab.savefig('elephant.png')