  function ga=gf(u1)
  global alpha
  ga=u1*(1+alpha^2)/(u1^2+alpha^2);