function f = myfun(x,a);

%f = a(2) * x.^2.*sqrt(1-1./x.^2).*exp(-x/a(1));

f = a(1) * sin(2*pi*x/a(2)+a(3)) .* exp(-x/a(4)) + a(5);
