From c888876054d713803ceaa54c4d5326942f553e53 Mon Sep 17 00:00:00 2001
From: "Harris M. Snyder" <harris.snyder@gmail.com>
Date: Sun, 22 May 2022 22:02:06 -0400
Subject: [PATCH 1/2] fix github issue 619

---
 src/stdlib_quadrature_gauss.f90     |  2 --
 src/tests/quadrature/test_gauss.f90 | 22 +++++++++++++++++++++-
 2 files changed, 21 insertions(+), 3 deletions(-)

diff --git a/src/stdlib_quadrature_gauss.f90 b/src/stdlib_quadrature_gauss.f90
index 0a346db48..fd1afa1e6 100644
--- a/src/stdlib_quadrature_gauss.f90
+++ b/src/stdlib_quadrature_gauss.f90
@@ -56,8 +56,6 @@ pure module subroutine gauss_legendre_fp64 (x, w, interval)
         if (present(interval)) then
             associate ( a => interval(1) , b => interval(2) )
             x = 0.5_dp*(b-a)*x+0.5_dp*(b+a)
-            x(1)       = interval(1)
-            x(size(x)) = interval(2)
             w = 0.5_dp*(b-a)*w
             end associate
         end if
diff --git a/src/tests/quadrature/test_gauss.f90 b/src/tests/quadrature/test_gauss.f90
index 8fce773e6..d5c843132 100644
--- a/src/tests/quadrature/test_gauss.f90
+++ b/src/tests/quadrature/test_gauss.f90
@@ -21,7 +21,8 @@ subroutine collect_gauss(testsuite)
             new_unittest("gauss-lobatto-analytic", test_gauss_lobatto_analytic), &
             new_unittest("gauss-lobatto-5", test_gauss_lobatto_5), &
             new_unittest("gauss-lobatto-32", test_gauss_lobatto_32), &
-            new_unittest("gauss-lobatto-64", test_gauss_lobatto_64) &
+            new_unittest("gauss-lobatto-64", test_gauss_lobatto_64), &
+            new_unittest("gauss-github-issue-619", test_fix_github_issue619) &
             ]
     end subroutine
 
@@ -48,6 +49,25 @@ subroutine test_gauss_analytic(error)
 
     end subroutine
 
+    subroutine test_fix_github_issue619(error)
+        !> Error handling
+        type(error_type), allocatable, intent(out) :: error
+        integer :: i
+
+        ! test the values of nodes and weights
+        i = 5
+        block 
+            real(dp), dimension(i) :: x1,w1,x2,w2
+            call gauss_legendre(x1,w1)
+            call gauss_legendre(x2,w2,interval=[-1._dp, 1._dp])
+
+            call check(error, all(abs(x1-x2) < 2*epsilon(x1(1))))
+            if (allocated(error)) return
+            call check(error, all(abs(w1-w2) < 2*epsilon(w1(1))))
+        end block
+
+    end subroutine
+
     subroutine test_gauss_5(error)
         !> Error handling
         type(error_type), allocatable, intent(out) :: error

From 1bf4e33e5138a012e9530dcc8de32972fde1af2c Mon Sep 17 00:00:00 2001
From: "Harris M. Snyder" <harris.snyder@gmail.com>
Date: Tue, 24 May 2022 00:09:56 -0400
Subject: [PATCH 2/2] add github issue link to test doc string

---
 src/tests/quadrature/test_gauss.f90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/tests/quadrature/test_gauss.f90 b/src/tests/quadrature/test_gauss.f90
index d5c843132..dee5c9fd2 100644
--- a/src/tests/quadrature/test_gauss.f90
+++ b/src/tests/quadrature/test_gauss.f90
@@ -50,7 +50,7 @@ subroutine test_gauss_analytic(error)
     end subroutine
 
     subroutine test_fix_github_issue619(error)
-        !> Error handling
+        !> See github issue https://github.com/fortran-lang/stdlib/issues/619
         type(error_type), allocatable, intent(out) :: error
         integer :: i