dN3/dt = PN1 - Γ32N3
dN2/dt = Γ32N3 - Γ21N2
dN1/dt = -PN1 + Γ21N2
ΣdNi/dt = 0, i=1,2,3.
Where P:pumping rate and Γ:decay rate.
Goal: To get steady state values for Ni, i=1,2,3 and this is done by setting dNi/dt = 0 , i=1,2,3.
When we get each steady state value, we usually denote it as Nibar, i = 1,2,3.
A good exercise would be to obtain the steady state population inversion:
N2bar - N1bar > 0 , and from here we will get a constraint on the pumping rate.
You will see that that the greater P is w.r.t. the decay rate, the greater the population inversion, and hence the "gain", which will give us "lasing"
Also you can get all of the steady state values in terms of NTbar, using NTbar = ΣNibar; i=1,2,3.
If you need a schematic of this situation, here it is:(under the headline:3-level laser)
http://en.wikipedia.org/wiki/Population_inversion
p.s. I assumed that we have (photon flux)<<1, since it's not included in the 3-level population rate equations.