He trabajado extensamente en este tema (precios y cálculo de IV en producción) y creo que puedo ofrecer una opinión informada. En primer lugar, Mathworks - la empresa que crea Matlab - no es una firma comercial, por lo que probablemente no deberías fiarte tanto de su consejo.
Hay pocos modelos de precios de opciones en forma cerrada, y todos tienen limitaciones prácticas. El modelo de Barone-Adesi y Whaley (por favor, corrige la ortografía de los apellidos ya que estoy escribiendo de memoria) es una aproximación simple para opciones americanas pero lamentablemente no es muy preciso, y no tiene en cuenta los dividendos. Roll-Geske-Whaley trata los dividendos, pero no muy bien - hay situaciones de arbitraje que son posibles en el modelo. Ju y Zhong tienen otra aproximación pero nuevamente no es muy precisa. Finalmente, Bjerksund y Stensland parecen tener la mejor aproximación (versión de 2002, no 1993) pero aún así no resuelve el problema de los dividendos discretos.
En mi experiencia, el árbol es la mejor opción. Los árboles CRR son lentos pero Leisen y Reimer, 1995 idearon un esquema que converge mucho más rápido. Además, Mark Joshi creó su propio esquema binomial que converge ligeramente más rápido. En lugar de un dividendo discreto puedes usar un dividendo proporcional discreto - así no terminas con un árbol frondoso. Alternativamente, puedes probar con un árbol trinomial y un factor de descuento adicional te dará una mejor resolución sobre el dividendo, pero no encontré una mejora significativa en producción. En mi opinión, esa es la mejor combinación en términos de velocidad y precisión. Si buscas opiniones alternativas, echa un vistazo a estos dos artículos - http://www.nccr-finrisk.uzh.ch/media/pdf/ODD.pdf sobre el problema del dividendo discreto, y http://ssrn.com/abstract=1567218 sobre la fijación de precios de opciones americanas.
Aun así, las mejoras más importantes en velocidad vendrán de tu código: una técnica que se usa extensivamente en los servidores de comilla es la precomputación. Básicamente, continuamente calculas los precios para el precio de la acción $\pm 0.1$ y la volatilidad $\pm 0.01$, para que cuando el spot o la volatilidad se muevan, utilices el valor previamente calculado en caché (o interpolado a partir de los valores más cercanos).