<gh_stars>1-10
def Ms(T):
    """ Returns the saturation magnetization as a function 
    of temperature
    """
    pass


def Kc1(T):
    """ Returns the 1st order cubic anisotropy term as a 
    function of temperature
    """
    pass


def Kc2(T):
    """ Returns the 1st order cubic anisotropy term as a 
    function of temperature
    """
    pass