*
* $Id: tprob.F,v 1.1.1.1 1996/02/15 17:55:02 mclareni Exp $
*
* $Log: tprob.F,v $
* Revision 1.1.1.1  1996/02/15 17:55:02  mclareni
* Kernlib
*
*
#include "sys/CERNLIB_machine.h"
#include "_kerngent/pilot.h"
      SUBROUTINE TPROB

#include "mkcde.inc"

      DIMENSION APROB(16),NU(10),CHI2(16,10),CHI2A(16,5),CHI2B(16,5)
      EQUIVALENCE (CHI2(1,1),CHI2A(1,1)), (CHI2(1,6),CHI2B(1,1))

      DATA APROB/  .995, .990, .975, .950, .900, .750, .500, .250,
     1             .100, .050, .025, .010, .005, .001, .0005, .0001/

      DATA NU / 1, 2, 5, 10, 25, 30, 40, 60, 80, 100/

      DATA CHI2A /
     A   3.92704E-5, 1.57088E-4, 9.82069E-4, 3.93214E-3, .0157908,        1A
     B   .101531, .454937, 1.32330, 2.70554, 3.84146, 5.02389, 6.63490,
     C   7.87944, 10.8280, 12.1160, 15.1370,
     A   1.00251E-2, 2.01007E-2, 5.06356E-2, .102587, .210720, .575364,   2A
     B   1.38629, 2.77259, 4.60517, 5.99147, 7.37776, 9.21034, 10.5966,
     C   13.8160, 15.2020, 18.4210,
     A    .41174, .55430, .831211, 1.145476, 1.61031, 2.67460, 4.35146,   5A
     B   6.62568, 9.23635, 11.0705, 12.8325, 15.0863, 16.7496, 20.5150,
     C   22.1050, 25.7450,
     A   2.15585, 2.55821, 3.24697, 3.94030, 4.86518, 6.73720, 9.34182,   10A
     B   12.5489, 15.9871, 18.3070, 20.4831, 23.2093, 25.1882, 29.5880,
     C   31.4200, 35.5640,
     A   10.5197, 11.5240, 13.1197, 14.6114, 16.4734, 19.9393, 24.3366,   25A
     B   29.3389, 34.3816, 37.6525, 40.6465, 44.3141, 46.9278, 52.6200,
     C   54.9470, 60.1400/

      DATA CHI2B /
     A   13.7867, 14.9535, 16.7908, 18.4926, 20.5992, 24.4776, 29.3360,   30A
     B   34.7998, 40.2560, 43.7729, 46.9792, 50.8922, 53.6720, 59.7030,
     C   62.1620, 67.6330,
     A   20.7065, 22.1643, 24.4331, 26.5093, 29.0505, 33.6603, 39.3354,   40A
     B   45.6160, 51.8050, 55.7585, 59.3417, 63.6907, 66.7659, 73.4020,
     C   76.0950, 82.0620,
     A   35.5346, 37.4848, 40.4817, 43.1879, 46.4589, 52.2938, 59.3347,   60A
     B   66.9814, 74.3970, 79.0819, 83.2976, 88.3794, 91.9517, 99.6070,
     C   102.695, 109.503,
     A   51.1720, 53.5400, 57.1532, 60.3915, 64.2778, 71.1445, 79.3343,   80A
     B   88.1303, 96.5782, 101.879, 106.629, 112.329, 116.321, 124.839,
     C   128.261, 135.783,
     A   67.3276, 70.0648, 74.2219, 77.9295, 82.3581, 90.1332, 99.3341,   100A
     B   109.141, 118.498, 124.342, 129.561, 135.807, 140.169, 149.449,
     C   153.167, 161.319/

      CALL NEWGUY ('PROB.','TPROB   ')
      ZERLEV = ZEROV(5)

         DO 50 ILOT=1,10
         DO 40 IC=1,16
   40 A(IC)=PROB (CHI2(IC,ILOT),NU(ILOT))
   50 CALL MVERIF (NU(ILOT),A,APROB,16)
      RETURN
      END