double probabilityDensity(double x, double mean, double variance ) { return ( 1 / math.sqrt( 2 * math.pi * variance ) ) * math.exp( - ( math.pow( ( x - mean ), 2 ) ) / ( 2 * variance ) ); }