Presentamos un modelo numérico práctico para calcular la transmitancia y reflectancia de recubrimientos antirreflejos multicapa teniendo en cuenta las no linealidades ópticas de tercer orden. De esta manera, se investiga el impacto de diferentes tipos de discretización del perfil del índice de refracción complejo en el rendimiento del sistema predicho. Además, se discuten aspectos de la paralelización de los cálculos. Se muestra que la inclusión de la no linealidad es esencial cuando intensidades de láser grandes inciden en el recubrimiento. El método desarrollado se aplica al diseño de diferentes recubrimientos antirreflectantes que coinciden con varios tipos de objetivos.