Kolmogorov Smirnov distribution function
p = pks(x)
vector or matrix
Prob( KS < x )
Adapted from Press, Teukolsky, Vetterling and Flannery, Numerical Recipes in C p620.