diff --git a/test/fc_lapack/src/example.F90 b/test/fc_lapack/src/example.F90
index 07a5df6..09f2b26 100644
--- a/test/fc_lapack/src/example.F90
+++ b/test/fc_lapack/src/example.F90
@@ -1,24 +1,40 @@
 program example
+
    implicit none
+
    integer, parameter :: n = 3
-   integer :: info
+   real(8), parameter :: small = 1.0d-12
+   integer :: ierr
    real(8) :: a(n, n)
    real(8) :: b(n)
    integer :: ipiv(n)
-   logical :: roots_ok=.false.
-  
-   data a /2,1,3,2,6,8,6,8,18/
-   data b /1,3,5/
+   logical :: roots_ok
 
-   call dgesv( n, 1, a, n, ipiv, b, n , info )
-   if (info.gt.0) stop "error in dgesv routine !"
+   a(1, 1) =  2.0d0
+   a(2, 1) =  1.0d0
+   a(3, 1) =  3.0d0
+   a(1, 2) =  2.0d0
+   a(2, 2) =  6.0d0
+   a(3, 2) =  8.0d0
+   a(1, 3) =  6.0d0
+   a(2, 3) =  8.0d0
+   a(3, 3) = 18.0d0
 
-!    roots are -0.5, 0.25, 0.25
-   roots_ok=dabs(b(1)+0.5).le.1d-12.and.dabs(b(2)-0.25).le.1d-12.and.dabs(b(3)-0.25).le.1d-12
+   b(1) = 1.0d0
+   b(2) = 3.0d0
+   b(3) = 5.0d0
+
+   call dgesv(n, 1, a, n, ipiv, b, n, ierr)
+   if (ierr /= 0) stop "error in dgesv routine!"
+
+   roots_ok = dabs(b(1) + 0.50d0) <= small .and. &
+              dabs(b(2) - 0.25d0) <= small .and. &
+              dabs(b(3) - 0.25d0) <= small
 
    if (roots_ok) then
-    print *, 'dgesv test ok'
+      print *, 'dgesv test ok'
    else
-    stop 'dgesv test failed!'
+      stop 'ERROR: dgesv test failed!'
    endif
+
 end program