En este trabajo desarrollamos un nuevo algoritmo para el cálculo eficiente de la función de error de Voigt/complex. En particular, en este enfoque proponemos un esquema de dos dominios donde el número de puntos de la cuadrícula de interpolación depende del parámetro de entrada. El análisis de error que realizamos muestra que la implementación en MATLAB cumple con los requisitos para aplicaciones de transferencia radiativa que involucran la base de datos espectroscópica molecular HITRAN. La prueba de tiempo de ejecución muestra que esta implementación en MATLAB proporciona cálculos rápidos, especialmente en rangos más pequeños del parámetro.