{
  if ($1 == "ATOM" && $3 == "CA") {
    nres = nres + 1
    x[nres] = substr($0,31,8)
    y[nres] = substr($0,39,8)
    z[nres] = substr($0,47,8)
            
#    print "nres",nres," x",x[nres],y[nres],z[nres]
  }
}END{
#  print "nres",nres
  sigma = 2.50
  for (i=1;i<=nres;i++) {
    for (j=1;j<i;j++) {
       dij = sqrt ((x[i]-x[j])^2+(y[i]-y[j])^2+(z[i]-z[j])^2)
#       print i,j,dij
       ibin = int(dij)
       h[ibin] = h[ibin] + 1
       for (k=0;k<=100;k++) {
          dd = k+0.5
#  print k,dd,dij,dd-dij,(dd-dij)^2/(2*sigma^2)
          if ((dd-dij)^2/(2*sigma^2)<100) hh[k] = hh[k] + exp(-(dd-dij)^2/(2*sigma^2)) 
       } 
     }
   }
  norm = 0;
  nnorm = 0;
  for (k=0;k<=100;k++) {
    norm = norm + h[k];
    nnorm = nnorm + hh[k];
  }
  for (k=0;k<=100;k++) printf("%10.5f%15.5e%15.5e\n",k+0.5,h[k]/norm,hh[k]/nnorm)
}