import numpy as np
from dysmalpy.models import LightTruncateSersic
import matplotlib.pyplot as plt

plt.figure()
plt.subplot(111, xscale='log', yscale='log')
ls1 = LightTruncateSersic(r_eff=5, n=1, r_inner=1, r_outer=20, L_tot=1.e11, tracer='halpha')
r=np.arange(0, 100, .01)

for n in range(1, 10):
     ls1.n = n
     plt.plot(r, ls1(r), color=str(float(n) / 15))

plt.axis([0.8, 27, 1e5, 1e10])
plt.xlabel('log Radius [kpc]')
plt.ylabel('log Intensity Surface Density [log Lsun/kpc^2]')
plt.text(1.1, 7.e8, 'n=1')
plt.text(1.1, 3.e9, 'n=10')
plt.show()