2630 DIFF(I)=(YACT(I)-YEST(I))»*2 



2649 IB SSR=SSR+DIFF(I) 



2650 NRITEI6,15) 



2660 15 F0RI1ftT(//,l)(," XVALUE YVALUE VEST DIFF »,/,) 



2670 DO 25 1=1, N 



2680 «RITE(6,20) X(I),YACT(I),YEST(I),SQRT(DIFF(I)) 



2690 20 F0RMAT(11(,F11.4,F11.4,F11.4,F11.4,/,) 



2700 25 CONTINUE 



2710 WRITE(6,40) CORR.SSR 



2720 40 FORHAK/.IX, "NON-LINEAR CORRELATION IS",51(,F10.B,/, 



2730 i ll(,"SUI1 SQR RESIDUALS IS",10X,F10.5,//) 



2740 RETURN 



2750 END 



2760 



2770C SUBROUTINE TO EVALUATE THE SAMMA FUNCTION 



2780C PROGRAM ADJUSTS ALPHA TO BE BETWEEN 1.0 AND 2.0 



2790C AND THEN MULTIPLIES BY 6F TO COMPENSATE 



2800 SUBROUTINE 6AMMA(ALPHA,AREA) 



2B10 DOUBLE PRECISION C(25), SUM 



2820 eF=1.0 



2830 IF(ALPHA) 1,2,3 



2840 



2850 2 PRINT, 'TROUBLE IN GAMMA' 



2860 AREA=1.0 



2870 GO TO 200 



2880 



2890C FOR GAMMA OF A NEGATIVE NUMBER 



2980 3 «=INT (ALPHA) 



2910 EPSI=ALPHA-FLOAT(M) 



2920 IF( M .EQ. 0) GF=GF/ALPHA 



2930 IF( M .EQ. 0) ALPHA=ALPHA+1.0 



2940 IF( M .EQ. 0) 60 TO 100 



2950 IF( M .EQ. 1) BF=1.0 



2960 IF( M .EQ. 1) GO TO 100 



2970 DO 10 I=2,M 



2980 10 GF=6Fi (FLOAT (I-1)+EPSI) 



2990 ALPHA=1.0+EPSI 



3000 80 TO 100 



3010 



3020C FOR GAMMA OF A POSITIVC NUMBER 



3030 1 M=INT (ALPHA) 



3040 EPSI=ALPHA-FLOAT(M) 



3050 IF( M .EQ. 0) 6F=1.0/(EPSI*(EPSI+1.0)) 



3060 IF( M .EQ. 0) ALPHA=EPSI+2.0 



3070 IF( M .EQ. 0) GO TO 100 



3080 DO 20 I=l,2-M 



3090 J=M+(I-1) 



3100 20 6F=GF/(EPSI+FL0AT(J)) 



3110 ALPHA=EPSI+2.0 



3120 



3130C COEFFICIENTS FOR SERIES EXPANSION OF THE GAMMA INTEGRAL 



3140C SEE HANDBOOK OF MATHEMATICAL FUNCTIONS BY ABRAMOHITZ AND SE6UN 



3150 100 C(l)=l. 



C11 



