isochrones/starmodel_old.py
def q_prior(q, m=1, gamma=0.3, qmin=0.1):
    """Default prior on mass ratio q ~ q^gamma
    """
    if q < qmin or q > 1:
        return 0
    C = 1/(1/(gamma+1)*(1 - qmin**(gamma+1)))
    return C*q**gamma
Default prior on mass ratio q ~ q^gamma
    
Default prior on mass ratio q ~ q^gamma
Default prior on mass ratio q ~ q^gamma
def q_prior(q, m=1, gamma=0.3, qmin=0.1):
    if q < qmin or q > 1:
        return 0
    C = 1/(1/(gamma+1)*(1 - qmin**(gamma+1)))
    return C*q**gamma
def q_prior ( q , m = 1 , gamma = 0.3 , qmin = 0.1 ) : if q < qmin or q > 1 : return 0 C = 1 / ( 1 / ( gamma + 1 ) * ( 1 - qmin ** ( gamma + 1 ) ) ) return C * q ** gamma
Default prior on mass ratio q ~ q^gamma
Default prior on mass ratio q ~ q^gamma
"""Default prior on mass ratio q ~ q^gamma
    """
q
m
gamma
qmin