source('utils.R')

T = 0.5
K = 35
C = 2.15
t = T - 7/52
r = 0.07
S0 = 36.12

res = BS_Call_Implied_Volatility(C, S0, t, r, K, T)
print(sprintf('Implied vol= %f', res[[1]]))