(use srfi-4 blas atlas-lapack)

(define order ColMajor)
(define n 4)
(define nrhs 1)

 
(define A (f64vector 1.8   5.25   1.58 -1.11  ;; column-major order
 		     2.88  -2.95 -2.69 -0.66 
 		     2.05  -0.95 -2.90 -0.59 
 		     -0.89 -3.80 -1.04  0.80))
(define b (f64vector 9.52 24.35 0.77 -6.22))
 
;; A and b are not modified
(define-values (LU x piv) (dgesv order n nrhs A b))

(pp (list LU x piv))