  function gi=gf(u1)
  gi=u1;