function gpic = gpic(x)
  gpic = 1. + 1.0*(x/(1. + x));