A recently developed, simple and easily programmable method based on the resonance technique for the analysis of nonlinear guided modes for fiber geometry is extended to include more complicated configurations. As an example of the proposed methodology, the following configuration is treated: a circular optical fiber consisting of a nonlinear core bounded by a nonlinear cladding medium. The calculation of the nonlinear modes in this case has some peculiarities and is presented in detail. The nonlinearity considered is of the more general type, including the saturable mode, which is of current significant importance. The numerical results obtained by this method correspond to the weak-guidance approximation and include both field distributions and power calculations. Some useful conclusions are deduced.