Хорошо.
Значит, Вы все знаете?
Самый простой путь, кмк, такой: пишете кубик,
который на вход принимает обычные рыночные переменные (время, цена БА, что еще может быть нужно)
и имеет параметры, которые собственно управляют формой улыбки.

На выходе этот кубик дает ЦЕНЫ опционов.
После этого стандартным кубиком мы можем превратить ЦЕНЫ опционов в их волатильность.

Пример кода с кубиком улыбки приведен чуть Выше в этой ветке.
Вам нужно только его модифицировать под формулы Хестона.


ПС Можем продолжить разговор в личку, если так будет удобней.