A partir de un modelo de interaccion de cuatro ondas basado en las ecuaciones de Kukhtarev para la respuesta de un material fotorrefractivo a un patron de intensidad incidente, y en la ecuacion de onda, se analiza numericamente el intercambio de energia entre los haces incidentes y los dos armonicos mas cercanos +1 y +2 en un cristal BSO. El modelo permite una interpretacion agil, util en el diseno de esquemas experimentales.