for i=1:5
v=simcall(10,.05,.2,63,10,2);
 disp(['standard error ' num2str(sqrt(var(v)/length(v)))])
disp(['strike any key to continue'])
 pause
end
v=simcall(10,.05,.2,63,10,100);

disp(['standard error ' num2str(sqrt(var(v)/length(v)))])

[call,put]=blsprice(10,10,.05,63/250,.2,0)