Merge pull request #22 from jacobwilliams/21-bug-fix

fixed bug in dqnc79.

Author: jacobwilliams

src/quadpack_generic.F90

@@ -7827,8 +7827,8 @@ subroutine dqnc79(fun,a,b,err,ans,ierr,k)
     procedure(func) :: fun !! function subprogram defining the integrand function `f(x)`.
     real(wp),intent(in) :: a !! lower limit of integration
     real(wp),intent(in) :: b !! upper limit of integration (may be less than `A`)
-    real(wp),intent(out) :: err !! a requested error tolerance.  Normally, pick a value
-                                !! `0 < ERR < 1.0e-8`.
+    real(wp),intent(in) :: err !! a requested error tolerance.  Normally, pick a value
+                               !! `0 < ERR < 1.0e-8`.
     real(wp),intent(out) :: ans !! computed value of the integral.  Hopefully, `ANS` is
                                 !! accurate to within `ERR *` integral of `ABS(FUN(X))`.
     integer,intent(out) :: ierr !! a status code: