function y = func(x,p,w) if w~=2 y=(abs(sin(w*pi*x)))^p; else y=sqrt(sin(2*pi*x)) - (2^(1/2)*pi^(1/2)*x^(1/2)-1/3*2^(1/2)*pi^(5/2)*x^(5/2)); end