Resolvimos numéricamente el conjunto de ecuaciones diferenciales parciales no lineales del material y, usando estas soluciones, incluimos la no uniformidad de la rejilla y de la magnitud y de fase de la modulación de la luz a lo largo del espesor de la muestra para calcular auto consistentemente el intercambio de energía del mezclado de dos ondas. Optimizamos la ganancia considerando condiciones fuertemente no-lineales, la variación del tamaño de la rejilla, actividad óptica, birrefringencia, absorción, ángulos de polarización, campos aplicados, y dos orientaciones cristalinas: el vector de la rejilla paralelo y perpendicular a la cara [001]. Bajo estas condiciones existe una relación compleja entre todos estos parámetros, y la predicción de las condiciones para el valor óptimo de la ganancia no es simple. Reportamos el espesor óptimo de la muestra para diferentes situaciones. Obtenemos una ganancia máxima de 5.2.We solved numerically the set of non-linear differential material rate differential equations, and using these solutions, we include the Non-uniformity of the grating and of the magnitude and phase of light modulation along a sample thickness to calculate self-consistently the energy exchange in two-wave mixing. We optimize the gain, considering strong nonlinear conditions, variation of fringe period, optical activity, birefringence, absorption, polarization angle, applied fields and two crystal orientations: the grating vector parallel and perpendicular to the face [001] Under these conditions there is a complex relationship among all these parameters, and the prediction of the conditions for the optimum value of the gain is not simple We report the optimal sample thickness for different situations We obtained a maximum gain of 5.2.