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