in 0.6.0
  - sync with GSL 1.9
    - Nonsymmetric eigensystems in Gsl_eigen
    - Basis splines
    - misc. other things, cf. GSL 1.9 NEWS file

in 0.5.3
  - compile with ocaml 3.10
  - fix GC bugs with C functions taking a callbak parameter (min, multimin, 
    multifit, etc.)

in 0.5.2
  - fix Gsl_sf.legendre_array_size

in 0.5.1
  - fix wrong declaration of externals in gsl_wavelet (Will M. Farr)
  - rewrite the blurb about callbacks raising exceptions (see below)
  - drop the mode argument Gsl_cheb.eval ("not for casual use")
  - fix GC bug in Gsl_cheb

in 0.5.0
  - sync with GSL 1.8
    - fix Gsl_cdf
    - wrap new functions: 
       gsl_multifit_linear_est, gsl_linalg_cholesky_decomp_unit,
       gsl_ran_gamma_mt, gsl_ran_gaussian_ziggurat, 
       gsl_sf_debye_5, gsl_sf_debye_6
  - sync with GSL 1.7
    - add Gsl_randist.binomial_knuth, add Gsl_multifit._linear_svd
in 0.4.1
  - make several record types 'private'. This disallows building bogus values
    that could crash the program.
  - remove several unsafe manual memory-management functions from
    module signatures
  - added Associated Legendre Polynomials and Spherical Harmonics in
    Gsl_sf (Will M. Farr)
  - fixed the type Gsl_fft.Complex.direction (Joe Shirron)
  - fixed the type of Gsl_randist.negative_binomial and negative_binomial_pdf 
    (Martin Willensdorfer)

in 0.4.0
  - sync with GSL 1.6
    - new module Gsl_wavelet for DWT
    - new module Gsl_sort
    - in Gsl_linalg, support for LQ and P^T LQ decompositions
    - add a mode argument in Gsl_cheb.eval
    - add RK2SIMP in Gsl_odeiv.step_kind
    - a couple of other small additions
  - bugfix for functions in Gsl_sf returning a result_e10
  - better support for Cygwin (James Scott, Lexifi)

in 0.3.5
  - improve build system a bit so that it works better on cygwin 
    (thanks to Brian Wilfley)
  - fix bugs in Gsl_odeiv (thanks to Will M. Farr)

in 0.3.4
  - fix a GC bug in the error handler, simplify exception raising code
  - add Qrng.dimension and Qrng.sample

in 0.3.3
  - report an error when building on a platform with ARCH_ALIGN_DOUBLE
    defined
  - findlib support
	 
in 0.3.2
  - complex functions (contributed by Paul Pelzl)

in 0.3.1 :
  - bugfix in Gsl_interp.eval_array
  - mlgsl_ieee.c now compiles with gcc 2.9x
  - build system improvements

in 0.3.0 :
  - sync with GSL 1.4
    - new module Gsl_cdf for cumulative distributions
    - new function Gsl_randist.binomial_tpe
  - compiles with MSVC (contributed by Lexifi)
  - memory bugfix in adaptative integration routines
  - bugfix in Gsl_ieee.set_mode, added FPU status word querying
  - changed arguments order in Gsl_matrix.transpose :
      first arg is destination, second is source

in 0.2.2 :
  - sync with GSL 1.3
    - new multidimensional minimizer (Nelder Mead Simplex algorithm)
    - new random distributions : Dirichlet and multinomial
    - new function Gsl_math.fcmp for approximate floating point values
    comparisons
  - fixed some potential problems with the GC

in 0.2.1 :
  - Gsl_linalg.matmult is now Gsl_linalg.matmult
  - Gsl_matrix.mul is now Gsl_matrix.mul_elements (same for
  Gsl_matrix.div)
  - vector/matrix macros work with gcc 2.9x (old_gcc target
  in Makefile)

in 0.2 :
  - rewrote the vector/matrix modules to add single precision
  bigarrays and complex values
  - added complex functions in Gsl_linalg and Gsl_eigen
  - added Ordinary Differential Equations
  - added Simulated Annealing
  - added Statistics and Histograms

in 0.1.1 :
  - fixed install target in Makefile
  - fixed several C stub function names
  - fixed a bug in ext_quot in quot.ml