A method and system for beam shaping employing a non-collinear quasi phase-matched interaction in a crystal whose nonlinear coefficient was encoded by a computer generated hologram is provided herein. The same axis is used for both satisfying the phase-matching requirements and encoding the holographic information. This allows one-dimensional beam shaping using a very simple to fabricate nonlinear crystal pattern and two-dimensional beam shaping with high conversion efficiency. Both are demonstrated by converting a fundamental Gaussian beam into Hermite-Gaussian and Laguerre-Gaussian beams at the second harmonic in KTiOPO
4 and stoichiometric lithium tantalate. The suggested scheme enables broad wavelength tuning by simply tilting the crystal.