  function si=sf(u2)
  global h
  si=((h-1)*u2+h)/(2*(h-1)*u2+1);
