def load_data():
    """reads in lyman alpha data in the units we want
    returns k and p(k) in proper units
    """

def plot_data():
    """plot k vs p(k) observed
    """
    ...


def plot_pk():
    """ plot k vs p(k) theory
     k: array-like
        wave number k
     pk: array-like"""
    ...

def plot_transfer():
    """ 
    """
    ...
