The equations of motion for the propagation of finite amplitude elastic waves in crystals of tetragonal symmetry have been derived starting from the expression for the elastic strain energy. The equations have been solved for a finite amplitude sinusoidal wave propagating along the pure mode directions which are [100], [110] and [001] for the tetragonal group TI. The solutions corresponding to longitudinal wave propagation yield expressions for the amplitudes of the fundamental and generated second harmonic for these directions in terms of certain combinations of second and third order elastic constants of the medium. The results will aid the experimenter to determine these constants using ultrasonic harmonic generation technique.