sigma=0.3
f=bscall(sigma)-15 ;
h=1e-5; 

while abs(f)>1e-10

  f2 = bscall(sigma+h)-15; 
  sigma  = sigma - h*f/(f2-f) 
  f = bscall(sigma)-15 ;
  
end