end_prices = r0296122_simulateEndPrice(price, tau, sigma, days, N);
value = max(0, (mean(end_prices) - strike))* (1 + tau)^(-days/365);