/*   function to compute students t   */
/*   Uses Numerical Recipes routine betai   */
/*   Based on numerical recipes equation 6.3.9 (p181)   */

/*   	Parameters	*t = Students t Abssicia */
/* 			*v = Degrees of freedom  */
/*	Value returned is prob(abs(t)>q) 	 */

double pst_(t,v)
double *t,*v;
{

	double betai();
        return 1.-betai((*v)/2., 0.5, (*v)/((*v)+(*t)*(*t)));
}