I'm quite surprised that two option pricers, ostensibly using identical closed end solving mechanisms, sometimes give different results.
You have at least two different packages offering option pricers based on the Black-Scholes formula in R:
I've generated some extremely wild option parameters for the the pricers:
(option type = put)
and got the following results:
One may say that the difference of 134.69 is not much for the option valued some 72.000 (134.69 / 71863.77 = 1.8742e-3). It's probably true.
But what's even more amazing is that both these results are different from what you can get from the "pure" Black-Scholes formula:
So, which value is correct?