1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
|
let pprint_histo { Gsl_histo.n = n ;
Gsl_histo.range = r ;
Gsl_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 = Gsl_histo.make n in
Gsl_histo.set_ranges_uniform h ~xmin ~xmax ;
begin try while true do
Scanf.scanf "%g"
(fun x -> Gsl_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))
|
