def madau_14(z, **kwargs):
    """Star formation history as a function of redshift, from
    Madau & Dickinson 2014 (http://arxiv.org/abs/1403.0007v3)

    """
    rate = (1+z)**2.7 / (1 + ((1+z)/2.9)**5.6)
    return rate