def METHOD_NAME(Ring, RFStation, dt, dE,
                                           timestep=0, Np_histogram=None):
    '''
    Returns the oscillation amplitude in time for given particle coordinates,
    assuming single-harmonic RF system and no intensity effects. 
    Optional: RF parameters at a given timestep (default = 0) are used.
    Optional: Number of points for histogram output
    '''
    omega_rf = RFStation.omega_rf[0, timestep]
    phi_rf = RFStation.phi_rf[0, timestep]
    phi_s = RFStation.phi_s[timestep]
    eta = RFStation.eta_0[0]
    T0 = Ring.t_rev[0]
    V = RFStation.voltage[0, 0]
    beta_sq = RFStation.beta[0]**2
    E = RFStation.energy[0]
    const = eta * T0 * omega_rf / (2. * V * beta_sq * E)
    dtmax = np.fabs(np.arccos(np.cos(omega_rf * dt + phi_rf) + const * dE**2)
                    - phi_rf - phi_s) / omega_rf
    if Np_histogram is not None:
        histogram, bins = np.histogram(dtmax, Np_histogram, (0,
                                                             np.pi / omega_rf))
        histogram = np.double(histogram) / np.sum(histogram[:])
        bin_centres = 0.5 * (bins[0:-1] + bins[1:])
        return dtmax, bin_centres, histogram
    else:
        return dtmax