from mmf.physics.seq import ScarfII
x = np.linspace(-3,3,100)
for B in [-1,0,1,2]:
   V = ScarfII(A=1, B=B)
   plt.plot(x, V(x), label='B=%i'%(B,))

plt.legend()