[...] = SBProd (...)
multiplies a symmetric banded matrix with a matrix
X=SBProd(A,B)
A is mxt where t-1 is number of non-zero superdiagonals
B is mxn
X is mxn
if A would be ! 11000 ! then A= ! 11 !
! 14300 ! ! 43 !
! 03520 ! ! 52 !
! 00285 ! ! 85 !
! 00059 ! ! 90 !
B is full matrix Ax=B