En este trabajo se dedujo un nuevo conjunto de ecuaciones de diferencias finitas basadas en la serie de Taylor para el calculo de propiedades opticas no lineales de sistemas moleculares. Las expresiones se obtuvieron a partir de las ecuaciones generalizadas de diferencias finitas con grado y exactitud arbitrarios, corrigiendo el orden de magnitud del error de truncamiento en las componentes no-axiales de las hiperpolarizabilidades. Las mismas se validaron calculando las propiedades opticas no lineales en 2 grupos de moleculas cuyos valores de primera y segunda hiperpolarizabilidad se han medido experimentalmente. Las moleculas fueron la 4-nitro-anilina, 4-ciano-fenol. Los valores de energia de cada sistema en funcion del campo electrico se calcularon por metodos de teoria del funcional de la densidad. Los resultados se compararon con los obtenidos con las ecuaciones de Kurtz y Kamada. El mejor desempeno se obtuvo con los funcionales HSEH1PBE, MN12SX y N12SX. En conclusion, las ecuaciones generalizadas representan una alternativa viable para el calculo de β y γ con precision y exactitud ajustables.