question: public static double cdf ( double x , double k , double theta , double shift ) { if ( x <= shift ) { return 0. ; } final double e = Math . exp ( ( x - shift ) * theta ) ; return e < Double . POSITIVE_INFINITY ? GammaDistribution . regularizedGammaP ( k , e ) : 1. ; }
target:  The CDF, static version.