功能: 求解对称3带宽带状矩阵方程

格式: SparseSolveSym3(A, B)

A : 一个列数为2的矩阵变量，且行数不少于3的矩阵变量
B : 一个与A行数一致的矩阵变量

说明:
1. 本函数主要求解如下线性方程组里的X
$$
AX = B
$$
这里A具有如下对称且带宽为3的形式
$$
A=\begin{bmatrix}
b_1& c_1
\\
c_1 & b_2 & c_2
\\
&c_2&b_3&c_3
\\
&&c_3&b_4&c_4
\\
...
\\
&&&&c_{n-2}&b_{n-1}&c_{n-1}
\\
&&&&&c_{n-1}&b_n
\end{bmatrix}
$$
由于A对称，因此函数里的A约定按如下方式存储
$$
A=\begin{bmatrix}
b_1& c_1
\\
b_2 & c_2
\\
b_3&c_3
\\
...
\\
b_{n-1} & c_{n-1}
\\
b_n&0
\end{bmatrix}
$$

2. B参数如果有多列，每一列单独进行求解

3. 当矩阵A分解失败时,本函数内部会采用拟牛顿法进行求解。

例子:
现在有如下A,b矩阵,解方程
$$
\begin{cases}
A&=\begin{bmatrix}
19&14
\\
14&13&8
\\
&8&9&6
\\
&&6&25&7
\\
&&&7&16&9
\\
&&&&9&21
\end{bmatrix}
\\
\\
b&=\begin{bmatrix}1.2&2.3&4&-5&6&-7.8\end{bmatrix}^T
\end{cases}
$$
b = [ 1.20000000000000
      2.30000000000000
      4.00000000000000
     -5.00000000000000
      6.00000000000000
     -7.80000000000000 ];
a = [ 19.0000000000000    14.0000000000000
      13.0000000000000    8.00000000000000
      9.00000000000000    6.00000000000000
      25.0000000000000    7.00000000000000
      16.0000000000000    9.00000000000000
      21.0000000000000    0.00000000000000 ];
x = SparseSolveSym3(a, b)//回车得到如下结果
x =
[-0.69053934994159
  1.02287483206359
 -0.16622773970554
 -0.44782483319313
  1.02756960972310
 -0.81181554702418 ]
