double dx=0.001;
  for(double x=-1; x<+1; x+=dx) {
    double ll=0;

    for(int i=0; i<G.size(); i++) {
      if ((dx - abs(G[i] - x))>0)
    ll++;
    }
    cout << x << " " << ll/(double)G.size() << endl;
  }