7
void compute_root_from_lambda(double, double *, double *);
11
compute_root_from_lambda(lambda, r, omega)
19
tmp = sqrt(3 + 144*lambda);
20
xi = 1 - 96*lambda + 24*lambda * tmp;
21
*omega = atan(sqrt((144*lambda - 1.0)/xi));
23
*r = (24*lambda - 1 - tmp2)/(24*lambda) \
24
* sqrt((48*lambda + 24*lambda*tmp))/tmp2;