This paper is the second and the last part of the study. It presents the results of tests of the numerical scheme proposed in the first part. The approximate procedure for estimation of the power of interacting waves at the exit from the crystal is discussed. This procedure allows the computation time to be shortened by more than an order of magnitude by means of a decrease in the minimal needed number of steps of the recurrent process.
second harmonic generation, system of nonlinear equations, numerical scheme