code 35060;
real procedure RECIP GAMMA(X, ODD, EVEN);
value X; real X, ODD, EVEN;
begin integer I;
real ALFA, BETA, X2;
array B[1:12];
B[ 1]:= -.28387 65422 76024; B[ 2]:= -.07685 28408 44786;
B[ 3]:= +.00170 63050 71096; B[ 4]:= +.00127 19271 36655;
B[ 5]:= +.00007 63095 97586; B[ 6]:= -.00000 49717 36704;
B[ 7]:= -.00000 08659 20800; B[ 8]:= -.00000 00331 26120;
B[ 9]:= +.00000 00017 45136; B[10]:= +.00000 00002 42310;
B[11]:= +.00000 00000 09161; B[12]:= -.00000 00000 00170;
X2:= X * X * 8;
ALFA:= -.00000 00000 00001; BETA:= 0;
for I:= 12 step - 2 until 2 do
begin BETA:= -(ALFA * 2 + BETA); ALFA:= - BETA * X2 - ALFA + B[I]
end;
EVEN:= (BETA / 2 + ALFA) * X2 - ALFA + .92187 02936 50453;
ALFA:= -.00000 00000 00034; BETA:= 0;
for I:= 11 step - 2 until 1 do
begin BETA:= -(ALFA * 2 + BETA); ALFA:= - BETA * X2 - ALFA + B[I]
end;
ODD:= (ALFA + BETA) * 2;
RECIP GAMMA:= ODD * X + EVEN
end RECIP GAMMA;
eop