Abstract: We explain, motivation behind this work and briefly describe foundation of new method which we have developed for efficient solution in PC environment of the nonlinear propagation equation with the boundary conditions applied for both circular and not circular transducers (like array). Comparison between new and old method will be presented for strongly nonlinear disturbance. At the end we will demonstrate the results of the numerical calculations of the nonlinear field propagating from the array.