def METHOD_NAME(self):
    '''
    Calculate the phase slippage resulting from the frequency offset for \
    each RF system
    '''
    delta_phi = (2 * np.pi * self.rf_station.harmonic[:, :self.end_turn]
                 * (self.rf_station.harmonic[:, :self.end_turn]
                     * self.new_frequency
                     - self.design_frequency)
                 / self.design_frequency)
    self.phase_slippage = np.cumsum(delta_phi, axis=1)