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

s = Sphere(c=[2.0, 0.5], R=[0.5, 1.0])
s.search_eigenfrequencies(2, np.arange(10.5, 15.0, 0.5)*np.pi)

r = np.linspace(0, 1, 1001)
pl.plot(r, s.eigenfunction(2, 0)(r), label='eigenfunction')
pl.plot(r, s.wavefunction(2, s[2,0]+0.5*s.Delta_omega)(r),
        label='not an eigenfunction')
pl.legend()
pl.xlabel('r')
pl.ylabel('f')