nmax = 8
lmax = 8
coeff = scf.compute_coeffs(flattened_hernquist_density, nmax=nmax, lmax=lmax,
                           M=M, r_s=a, args=(M,a,q), skip_m=True)
(S,Serr),(T,Terr) = coeff

plt.figure(figsize=(6,4))
plt.semilogy(np.abs(S[:,0,0]), marker=None, lw=2)
plt.xlabel("$n$")
plt.ylabel("$S_{n00}$")
plt.tight_layout()