  function sa=sf(u2)
  global      h     cm   beta    
  sa=2*cm*(h-beta)*u2/(cm^2+u2^2)+beta;
