def METHOD_NAME(phimax):
    '''
    Find the tune w.r.t. the central synchrotron frequency corresponding to a
    given amplitude of synchrotron oscillations in phase 
    '''
    return 0.5 * np.pi / ellipk(x(phimax))