U=unifrnd(0,1,1,100000);  
v=callopt2(U,10,10,.05,.2,.25);
disp(['average discounted payoff ' num2str(mean(v))])
Sf=sqrt(var(v));
disp(['Standard error of estimate ' num2str(Sf/sqrt(length(U)) )])
[CALL,PUT] =BLSPRICE(10,10,.05,.25,.2,0) 