From 3d10dd905f051e609deb6de825d5c4b0d4db1482 Mon Sep 17 00:00:00 2001
From: wo80 <christian.woltering@tu-dortmund.de>
Date: Wed, 22 Nov 2023 13:21:20 +0100
Subject: [PATCH] Fix examples BAND _nbdr2 returning error info = -9

Both examples access the M matrix, so we initialize M to be the identity matrix. See #443 for discussion.
---
 EXAMPLES/BAND/cnbdr2.f | 1 +
 EXAMPLES/BAND/znbdr2.f | 1 +
 2 files changed, 2 insertions(+)

diff --git a/EXAMPLES/BAND/cnbdr2.f b/EXAMPLES/BAND/cnbdr2.f
index 6e0d80bcc..30e089ea2 100644
--- a/EXAMPLES/BAND/cnbdr2.f
+++ b/EXAMPLES/BAND/cnbdr2.f
@@ -206,6 +206,7 @@ program cnbdr2
       idiag = kl+ku+1
       do 30 j = 1, n
          a(idiag,j) = (4.0E+0, 0.0E+0)  / h2
+         m(idiag,j) = (1.0E+0, 0.0E+0)
   30  continue
 c
 c     %-------------------------------------%
diff --git a/EXAMPLES/BAND/znbdr2.f b/EXAMPLES/BAND/znbdr2.f
index 61a0f6f39..b9b5377ec 100644
--- a/EXAMPLES/BAND/znbdr2.f
+++ b/EXAMPLES/BAND/znbdr2.f
@@ -206,6 +206,7 @@ program znbdr2
       idiag = kl+ku+1
       do 30 j = 1, n
          a(idiag,j) = (4.0D+0, 0.0D+0)  / h2
+         m(idiag,j) = (1.0D+0, 0.0D+0)
   30  continue
 c
 c     %-------------------------------------%