%FindBetaFromMu
function [beta] = FindBetaFromMu(mu,alpha)
sig = 1;
[I rho I1 I2 I3] = GaussRomberg(@(x) abs(x).^alpha, mu, sig, 1000);
[J rho J1 J2 J3] = GaussRomberg(@(x) abs(x).^alpha.*sign(x), mu, sig, 1000);
J = min(I, J);
J = max(-I, J);
beta = J/I;
end


