gsl-ocaml / examples / histo_ex.ml

open Gsl

let pprint_histo { Histo.n = n ; 
		   Histo.range = r ; 
		   Histo.bin = b } =
  for i=0 to pred n do
    Printf.printf "%g %g %g\n"
      r.(i) r.(succ i) b.(i)
  done


let main xmin xmax n =
  let h = Histo.make n in
  Histo.set_ranges_uniform h ~xmin ~xmax ;
  
  begin try while true do
    Scanf.scanf "%g" 
      (fun x -> Histo.accumulate h x)
  done
  with End_of_file -> () 
  end ;

  pprint_histo h
    

let _ = 
  if Array.length Sys.argv <> 4
  then (
    Printf.printf "Usage: gsl-histogram xmin xmax n\n" ;
    Printf.printf "Computes a histogram of the data on \
                   stdin using n bins from xmin to xmax\n" ;
    exit 1 ) ;
  main 
    (float_of_string Sys.argv.(1))
    (float_of_string Sys.argv.(2))
    (int_of_string Sys.argv.(3))
Tip: Filter by directory path e.g. /media app.js to search for public/media/app.js.
Tip: Use camelCasing e.g. ProjME to search for ProjectModifiedEvent.java.
Tip: Filter by extension type e.g. /repo .js to search for all .js files in the /repo directory.
Tip: Separate your search with spaces e.g. /ssh pom.xml to search for src/ssh/pom.xml.
Tip: Use ↑ and ↓ arrow keys to navigate and return to view the file.
Tip: You can also navigate files with Ctrl+j (next) and Ctrl+k (previous) and view the file with Ctrl+o.
Tip: You can also navigate files with Alt+j (next) and Alt+k (previous) and view the file with Alt+o.