macro geffic
  v/cre x(5) r .2 .5 .8 .9 .95
* v/cre y(5) r .835 .846 .868 .882 .893
  v/cre y(5) r .80258 .81382 .83705 .85047 .86253
  fun1 1000000 interp.f77 19 .775 .955
  h/copy 1000000 50 'geometric efficiency'
  v/cre geffic(19)
  get/c 50 geffic
  sigma o=array(19,19#1)
  sigma reffic=order(geffic,o)
  v/write reffic geffic.dat
return