Full_Name: Morten Welinder Version: 2.2.0 OS: Linux Submission from: (NULL) (216.223.241.229) The qt function for df=1 is implemented as... q = - tan((P+1) * M_PI_2); Adding 1 kills accuracy for P near 0. A better way is to use q = cot(P * M_PI_2); or q = 1/tan(P * M_PI_2);