question: public static double logcdf ( double val , double loc , double scale , double shape ) { val = ( val - loc ) / scale ; return Math . log1p ( Math . exp ( - val ) ) * - shape ; }
target:  log Cumulative density function. TODO: untested.