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