S0 = 100;
K = 95;
r = 0.12;
T = 1.0;
sigma = 0.2; 

R = 2;
L=100;
U=100;

M=10^6

tic();
price = BS_fourier_call(S0,K,r,T,sigma,R,L,U)
time1=toc()

tic();
trueprice = BS_call(S0,K,T,r,sigma)
time2=toc()

tic();
[MCprice, RMSE]=BS_MC_call(S0,K,T,r,sigma,M)
time3=toc()

disp([price trueprice])
disp([time1 time2])
disp([MCprice-2.58*RMSE MCprice MCprice+2.58*RMSE])
disp(time3)


