def imf_flat(x):
    """
    Compute a flat IMF (useful for simulations, not for normal BEAST runs)

    Parameters
    ----------
    x : numpy vector
      masses

    Returns
    -------
    imf : numpy vector
      unformalized IMF
    """
    return 1.0