function [x, y] = elephant(t, p1, p2, p3, p4, p5) % Draws the Fermi - Non Neumann elephant % using Fourier series and epicycloids npar = 6; Cx = zeros(npar, 1) + zeros(npar, 1)*1i; Cy = zeros(npar, 1) + zeros(npar, 1)*1i; Cx(1) = real(p1)*1i; Cx(2) = real(p2)*1i; Cx(3) = real(p3); Cx(5) = real(p4); Cy(1) = imag(p4) + imag(p1)*1i; Cy(2) = imag(p2)*1i; Cy(3) = imag(p3)*1i; x = [ev_fourier(t,Cx), -imag(p5)]; y = [ev_fourier(t,Cy), imag(p5)]; end