Google
 

Trailing-Edge - PDP-10 Archives - decuslib20-01 - decus/20-0025/gamxx.for
There is 1 other file named gamxx.for in the archive. Click here to see a list.
00010	C   GAMXX.SRC   OLS   32 LINES
00020	C
00030	      FUNCTION GAM(XPL1)
00040	1     X=XPL1-1.
00050	      TEMP=1.
00060	      TMP=XPL1
00070	      IF (X-1.)  11, 11, 101
00080	11    IF (X)  201, 301, 301
00090	101   N1=X
00100	      DO 121 J1=1,N1
00110	      TMP=TMP-1.
00120	121   TEMP=TEMP*TMP
00130	      GO TO 301
00140	201   N1=1.-X
00150	      DO 221 J1=1,N1
00160	      IF (TMP)  211, 241, 211
00170	211   TEMP=TEMP/TMP
00180	221   TMP=TMP+1.
00190	      GO TO 301
00200	241   GAM=0.
00210	      GO TO 5001
00220	301   TMP=TMP-1.
00230		GAM=.035868343*TMP-.193527818
00240		GAM=GAM*TMP+.482199394
00250		GAM=GAM*TMP-.756704078
00260		GAM=GAM*TMP+.918206857
00270		GAM=GAM*TMP-.897056937
00280		GAM=GAM*TMP+.988205891
00290		GAM=GAM*TMP-.577191652
00300		GAM=GAM*TMP+1.
00310		GAM=GAM*TEMP
00320	5001  END