anchor.diff:      call least_squares(X, Y, TX, B)
       
       print*, "--------------------------"
  -    print*, "      X             Y      "
  -    print*, "--------------------------"
  -    do i=1,m
  -        print*, X(i,2), Y(i)
  -    enddo
  -    print*, "--------------------------"
  -    print*, " NR based Coefficients"
  -    print*, B
  +    print*, "  Testing Regression "
  +    if (verbose) then
  +        print*, "--------------------------"
  +        print*, "      X             Y      "
  +        print*, "--------------------------"
  +        do i=1,m
  +            print*, X(i,2), Y(i)
  +        enddo
  +        print*, "--------------------------"
  +        print*, " NR based Coefficients"
  +        print*, B
  +    endif
   
       Y = Y*2 + 5
       call least_squares(X, Y, TX, B)
  -    print*, "--------------------------"
  -    print*, "--------------------------"
  -    print*, " Y = Y * 2 + 5"
  -    print*, "--------------------------"
  -    print*, " NR based Coefficients new Y values"
  -    print*, "--------------------------"
  -    print*, B
  +    if (verbose) then
  +        print*, "--------------------------"
  +        print*, "--------------------------"
  +        print*, " Y = Y * 2 + 5"
  +        print*, "--------------------------"
  +        print*, " NR based Coefficients new Y values"
  +        print*, "--------------------------"
  +        print*, B
  +    endif
       
       call lapack_least_squares(X, Y, B_lapack)
  -    print*, "--------------------------"
  -    print*, " LAPACK based regression "
  -    print*, "--------------------------"
  -    print*, B_lapack
  -    
  -    print*, "--------------------------"
  +    if (verbose) then
  +        print*, "--------------------------"
  +        print*, " LAPACK based regression "
  +        print*, "--------------------------"
  +        print*, B_lapack
  +        print*, "--------------------------"
  +    endif
       pass_fail = .True.
       do i=1,n
           if (abs(B_lapack(i)-B(i)) > MAX_ERROR) then
  
positive.code:     endif
      print*, "--------------------------"
      
  end program test_regression
  
negative.code:     integer, parameter :: n=2, m=10
      real :: random
      integer :: i
      logical :: pass_fail, verbose
  
      allocate(X(m, n))
      allocate(TX(n,m))
      allocate(Y(m))
      allocate(B(n))
      allocate(B_lapack(n))
      
      verbose = read_verbosity()
      
      do i=1,m
          X(i,1) = 1
  