Estoy tratando de calcular la volatilidad implícita libre del modelo $\sigma_{MF}$ para un índice de rendimiento relativo utilizando el siguiente método:
$\sigma_{MF}^2=2\sum_{i} [\frac{C(T,K_{i})}{K_{i}^2} - \frac{max(0,F-K_{i})}{K_{i}^2}]\Delta K_{i}$ , donde $F=Ie^{(\frac{\sigma_{M}^2-\sigma_{S}^2+\sigma_{MF}^2}{2})T}$
La única incógnita aquí es $\sigma_{MF}$ ¿Cómo puedo implementar esto usando Matlab? Estoy confundido en cuanto a cómo puedo utilizar funciones de optimización no lineal cuando la incógnita $\sigma_{MF}$ está dentro de un bucle.