I'm writing a program which generates the emission spectrum of an element with atomic number Z. To do this, I have used the equation:

$\frac{1}{\lambda}={R}_{\mathrm{\infty}}{Z}^{2}(\frac{1}{{n}_{1}^{2}}-\frac{1}{{n}_{2}^{2}})$

Where ${n}_{1}$ is the number of the original energy level and ${n}_{2}$ is the number of the final energy level and where ${n}_{1}>{n}_{2}$

However, this does not calculate the intensity of the lines. Is there anyway to model the intensity using ${n}_{1}$, ${n}_{2}$ and/or $Z$, assuming that the brightest line has an intensity of 1?