S0 = 100;
K = 130;
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();
//[price, CI99]=BS_MC_call(S0,K,T,r,sigma,M)
//time3=toc()
