Смотрю цены опционов что то высокoваты..., пересчитал по другому, действительно завышены процентов на 5% и даже больше...
Нашел ошибку, когда перевожу распределение из лог пространства нарушается дискретизация (шаги становятся неравномерны, и ошибка увеличивается) — внешне выглядит вроде одно и то же (слева), но если увеличить ошибку видна проблема (справа).

И, я маловато точек при численной интеграции использовал, что еще увеличило ошибку... распределение вообще лучше не переводить в экспоненциальную форму а использовать «d.pdf(x.log()) / x»
И, это уже 2я ошибка, я использовал PMF вместо PDF и там ошибки еще сильнее были, из за того что дискретные «массы» PMF еще неравномерней «перекашивались» после экспоненты.
В симуляции начинаются сложности при высоких размерностях, но, в финансах, очень многое низкой размерности.