function f = h(alpha,t)
f = integral(@(x) sin(x).*x.^(-1.3691),t,alpha*t);
