Example program for dgbtrs

This example solves the system of equations AX=B, where
A= -0.23 2.54 -3.66 0.00 -6.98 2.46 -2.73 -2.13 0.00 2.56 2.46 4.07 0.00 0.00 -4.78 -3.82   and   B= 4.42 -36.01 27.13 -31.67 -6.14 -1.16 10.50 -25.82 .  
Here A is nonsymmetric and is treated as a band matrix, which must first be factorized by dgbtrf.