      SUBROUTINE TUTEST(DATA1,N1,DATA2,N2,T,PROB)
      DIMENSION DATA1(N1),DATA2(N2)
      CALL AVEVAR(DATA1,N1,AVE1,VAR1)
      CALL AVEVAR(DATA2,N2,AVE2,VAR2)
      T=(AVE1-AVE2)/SQRT(VAR1/N1+VAR2/N2)
      DF=(VAR1/N1+VAR2/N2)**2/((VAR1/N1)**2/(N1-1)+(VAR2/N2)**2/(N2-1))
      PROB=BETAI(0.5*DF,0.5,DF/(DF+T**2))
      RETURN
      END
