def _bifur_ ( self , mu , sigma1 , sigma2 ) :
    """Generate the bifurcated gaussian
    >>> value = bifur ( 0 , -1 , +2 )    
    """
    if sigma1 * sigma2 > 0.0 :
        raise ValueError( 'Lower and upper errors must have opposite signs' )