>>>>> "C" == Conrad Schneiker <schneiker / jump.net> writes:

C>     http://sourceware.cygnus.com/gsl/ref/gsl-ref_toc.html

pigeon% cat b.rb
#!/usr/bin/ruby
require 'random'
a = Gsl::Rng.new()
4.times { p a.rand() }
p "================"
b = a.state
4.times { p a.rand() }
p "================"
a.state = b
4.times { p a.rand() }
p a.info
p a.gaussian
p a.gaussian(12)
p a.gaussian_tail(2, 12)
p a.gaussian_tail(12)
p a.bivariate_gaussian
p a.bivariate_gaussian(1)
p a.bivariate_gaussian(1, 1)
p a.bivariate_gaussian(1, 2, -1)
p a.exponential(1)
p a.laplace(1)
p "levy"
p a.levy(1, 1)
p a.levy(1, 4)
p "shuffle"
p a.shuffle([1,2,3,4,5,6,7,8,9,10])
p a.shuffle([1,2,3,4,5,6,7,8,9,10])
p a.choose(4, [1,2,3,4,5,6,7,8,9,10])
pigeon%

pigeon% b.rb
0.8173300598
0.9990608995
0.5103543724
0.1315329098
"================"
0.03541634837
0.9924695343
0.6257087034
0.06259194575
"================"
0.03541634837
0.9924695343
0.6257087034
0.06259194575
#<Struct::Rnginfo name="mt19937", min=0, max=4294967295, seed=0>
-1.037281406
-23.50162436
13.24839816
12.09667727
[0.5573033256, 0.2893942072]
[-2.10828341, 2.016359466]
[-0.6506783886, 0.6428933895]
[0.08091231178, -0.1618246236]
0.295286719
-0.6626309225
"levy"
0.7876817326
3.147055135
"shuffle"
[6, 1, 10, 2, 9, 3, 4, 7, 5, 8]
[7, 8, 5, 6, 4, 9, 2, 3, 10, 1]
[1, 3, 6, 8]
pigeon% 


 well, I've never finished the interface :-(, but if someone want the
 source for random I can send it


Guy Decoux