A theoretical model that allows numerical simulations of optical bistability phenomena in CW-pumped erbium-doped fiber lasers is presented. Results the match correctly with previous experimental data. In order to simulate this phenomenon, population inversion and powers longitudinal distributions along active have been considered. If long enough, it has found transition between non-lasing lasin...