def power_spectrum(k,a):
    """Simple function for power laws

    Args:
        k (float array): wave number
        a (float): power law

    Returns:
        float array: k^-a (note minus sign is assumed!)
    """
    return k**-a