t1=cputime;
	
U=rand(1,200000);  % crude Monte Carlo 
	
F=fn(U); % value of 200000 discounted values of the call options
	
m=mean(F) % average discounted value of call = .4743
	
   se=sqrt(var(F)/length(F)) %\ the standard error of the mean
   t1=[t1 cputime];

 %%%%%%%%%%%%  *****************************************
   
	
U=rand(1,100000);  %  use of antithetic random numbers
	
F=.5*(fn(U)+fn(1-U));
mean(F)	
 se=sqrt(var(F)/length(F)) %\ the standard error of the mean
t1=[t1 cputime];
display(['times ' mat2str(diff(t1))])

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


	
