RESUMEN Se realizo el estudio de la propagacion de un haz Gaussiano por medio de la implementacion del algoritmo denominado metodo de la propagacion del haz, Beam Propagation Method, BPM. En el estudio se consideran frentes de fase planos en la region paraxial. Para la implementacion del BPM se desarrollo un software en MATLAB que muestra la propagacion del haz en funcion de la distancia de propagacion. Para realizar la simulacion del haz, el algoritmo soluciona la ecuacion de difraccion de Fresnel. Para los calculos numericos se utilizo una longitud de onda de 632,8 nm y una cintura inicial del haz de 1 mm. A traves de la implementacion de procesos de convolucion por transformada de Fourier el programa muestra las amplitudes de los picos Gaussianos inicial y difractado, permitiendo obtener resultados coherentes con los reportados en la literatura. Palabras clave: BPM, Difraccion, Laser, Haz, Fourier.