A simple mean-field model of a nonlinear stellar α-ω dynamo is considered, in which dynamo action is supposed to occur in a spherical shell, and where the main nonlinearity retained is the influence of the Lorentz force on the zonal flow field. The equations are simplified by truncating in the radial direction, while full latitudinal dependence is retained. The resulting nonlinear p.d.e.'s in latitude and time are solved numerically, and it is found that while regular dynamo wave type solutions are stable when the dynamo number D is sufficiently close to its critical value, there is a wide variety of stable solutions at larger values of D. Furthermore, two different types of dynamo can coexist at the same parameter values. Implications for fields in late-type stars are discussed.