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

