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;
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.