import numpy as np
import matplotlib.pyplot as pl
from allium import Sphere

s = Sphere(c=[1.05, 1.0], R=[0.1, 1.0])
s.search_eigenfrequencies([0,1,2,3], np.arange(0.5, 20.0)*np.pi)
for l in range(4):
    w = np.array(s[l])/s.Delta_omega
    pl.plot(np.mod(w-0.01, 1.0)+0.01, w, 'o', label='ℓ=%i' % l)

pl.xlabel('ω/Δω mod 1')
pl.ylabel('ω/Δω')
pl.xlim(-0.05, 1.05)
pl.legend()