%%HP: T(3)A(R)F(.);
\<< \-> z1
  \<< z1 1 - \-> z
    \<<
      IF z RE 0 <
      THEN \pi \->NUM z
* DUP SIN / 1 z -
GAMMA /
      ELSE z 5.5 +
z .5 + ^ z 5.5 +
NEG EXP * 2 \pi \->NUM
* \v/ * 1 76.18009173
z 1 + / +
-86.50532033 z 2 +
/ + 24.01409822 z 3
+ / + -1.231739516
z 4 + / +
.00120858003 z 5 +
/ + -.00000536382 z
6 + / + *
      END
    \>>
  \>>
\>>
