def schmidt_law(Ak, theta):
    '\n    schmidt_law(Ak,(beta,kappa))\n    beta is the power law index (same as alpha)\n    '
    if (len(theta) == 2):
        (beta, kappa) = theta
        return (kappa * (Ak ** beta))
    elif (len(theta) == 3):
        (beta, kappa, Ak0) = theta
        sfr = ((Heaviside((Ak - Ak0)) * kappa) * (Ak ** beta))
        sfr[(Ak < Ak0)] = 0
        return sfr
schmidt_law(Ak,(beta,kappa))
beta is the power law index (same as alpha)
def schmidt_law(Ak, theta):
    '\n    schmidt_law(Ak,(beta,kappa))\n    beta is the power law index (same as alpha)\n    '
    if (len(theta) == 2):
        (beta, kappa) = theta
        return (kappa * (Ak ** beta))
    elif (len(theta) == 3):
        (beta, kappa, Ak0) = theta
        sfr = ((Heaviside((Ak - Ak0)) * kappa) * (Ak ** beta))
        sfr[(Ak < Ak0)] = 0
        return sfr