function quadresult

p=[.5,1,2];
w=[1,10,100];


for i=1:3
    for j=1:3
        alpha=p(i)
        omega=w(j)
        f=@(x)(abs(sin(w(j)*pi*x))).^p(i);
        quad(f,0,1)
    end
end

alpha=.5
omega=2
f=@(x)sqrt(abs(sin(2*pi*x)));
quad(f,0,1)