## test DJCM histogram program
   a = rand(2000,3) + rand(2000,3) + rand(2000,3);
   [H,X,G] = djcmhisto(a,50);
   f=1;   
##
##   hold off; plot( X(:,f), H(:,f) ); hold on; plot( X(:,f), G(:,f));
##   
##  or alternatively using the gplot function (which is faster than plot)

   XH =  [ X(:,f) , H(:,f) , G(:,f) ] ;
   gset yrange [0:]
   gplot XH u 1:2 t 'Data' w histeps,  XH u 1:3 t 'Gaussian fit' w l 4 ;
