A method for simulating the propagation of coherent optical beams with Gaussian irradiance patterns in free space and through unaberrated lenses is described. This method is based on the skew-ray representation of Gaussian beams. The rotation in space of a collimated laser beam with skew axis generates a Gaussian beam profile. The phase of the optical field is given by the angular position of the laser-beam center.