def METHOD_NAME(E, E0, sigma, alpha, k, beta):
    """
    return synthetic data of the form
    d = alpha * e ** (-(E - e0)**2 / (2 * sigma ** 2) + beta * sin(k * E)
    Parameters
    ----------
    E : ndarray
        The energies to compute values at
    E0 : float
       Location of the peak
    sigma : float
       Width of the peak
    alpha : float
       Height of peak
    k : float
       Frequency of oscillations
    beta : float
       Magnitude of oscillations
    """
    return alpha * np.exp(-((E - E0) ** 2) / (2 * sigma**2)) + beta * (1 + np.sin(k * E))