function [theta] = ouo_beta_M(gamma, beta, M)
% angoli in radianti
%
theta = atan(2.*cot(beta).*((M.*sin(beta)).^2-1)./(M.^2.*(gamma+cos(2*beta))+2));
end