We present an algorithm for calculating steady states in the dynamic PDE model for SLs admitting gain compression, spatial hole burning, and multilevel carrier rate equations. Presented example simulations rely on 1(time)+1(space)-dimensional traveling-wave- and Lang-Kobayashi-type models.