1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
#include <math.h>
#include <stddef.h>
#include <stdlib.h>
#include <string.h>
#include <stdio.h>
void compute_root_from_lambda(double, double *, double *);
void
compute_root_from_lambda(lambda, r, omega)
double lambda;
double *r;
double *omega;
{
double xi;
double tmp, tmp2;
tmp = sqrt(3 + 144*lambda);
xi = 1 - 96*lambda + 24*lambda * tmp;
*omega = atan(sqrt((144*lambda - 1.0)/xi));
tmp2 = sqrt(xi);
*r = (24*lambda - 1 - tmp2)/(24*lambda) \
* sqrt((48*lambda + 24*lambda*tmp))/tmp2;
return;
}
|
