With the introduction of harmonic imaging, the design of new array transducers for tissue and contrast imaging became indispensable. Hence, prior knowledge of harmonic beams is essential to attain optimal harmonic performances. For that purpose, a new numerical algorithm that solves the parabolic nonlinear wave equation is developed. The algorithm is based on finite differences and performs exclusively in time domain. Pulsed harmonic fields emitted by a medical transducer were measured and computed at different mechanical indices. Simulations and measurements showed very good agreement for all the harmonic components.