#### Iterated printing of cubehelix.gp
#### (why doesn't gnuplot provide multicolumn table plot?)
x255=1.0*i
x=x255/imax
r0=f0red(x);g0=f0grn(x);b0=f0blu(x)
r=fred(x);  g=fgrn(x);  b=fblu(x)
gray_ideal=fxgry(x)
gray_appar=fgry(x)
defect=gray_appar-gray_ideal
defred=r0-r; defgrn=g0-g; defblu=b0-b
if(defect<mindefect) mindefect=defect;xatmin=x
if(defect>maxdefect) maxdefect=defect;xatmax=x
if(defred<mindefred) mindefred=defred;xatmin=x
if(defred>maxdefred) maxdefred=defred;xatmax=x
if(defgrn<mindefgrn) mindefgrn=defgrn;xatmin=x
if(defgrn>maxdefgrn) maxdefgrn=defgrn;xatmax=x
if(defblu<mindefblu) mindefblu=defblu;xatmin=x
if(defblu>maxdefblu) maxdefblu=defblu;xatmax=x
if(abs(defect)>absdefect) absdefect=abs(defect);xatabs=x
if (r0<0) nrlo=nrlo+1
if (g0<0) nglo=nglo+1
if (b0<0) nblo=nblo+1
if (r0>1) nrhi=nrhi+1
if (g0>1) nghi=nghi+1
if (b0>1) nbhi=nbhi+1
isfalling=0;graydrop=0.
if (gray_appar<=gray_last) ndegrade=ndegrade+1;isfalling=1;\
   graydrop=gray_appar-gray_beforefalling
if (graydrop<worstgraydrop) worstgraydrop=graydrop
print sprintf("%5.f  %9.6f %9.6f %9.6f",x255,r,g,b)
#print sprintf("%9.6f  %9.6f %9.6f %9.6f",x,r,g,b)
if (isfalling==0) gray_beforefalling=gray_appar
gray_last=gray_appar
i=i+1
if (i<=imax) reread