function y = myfun1(x,c,n,delta)
%function y = myfun(x,c,n,delta)
    %1./(x.^3-2*x-5)
    %y = sin(1*pi*x/5).*(exp(-1.*((x - 5/2).^2))./1.5);
    
    z = -1*(x.^2 + (c/2)^2 - c*x);
    z = z/delta;
    y = 2*sin(n*c).*exp(z)/c;
    