功能: 矩阵乘法
 
格式: SparseMul(A, b)

A : 一个列数为3的系数矩阵
b : 一个行数与A一致的常数项矩阵

说明:
1. 本函数本质是求解如下方程组
$$
Ax=b
$$

2. 本函数主要稀疏矩阵进行乘法运算,其中要求A的第1行为非0元素行索引,第2行为非0元素列索引,第3行为对应元素值。

3. A当中存在非数字,或者索引值异常,或者存在重复索引,会返回失败。

4. A当中的行列索引,均从0开始。

例子:
A =
[ 40.0000000000000    5.00000000000000    0.00000000000000    0.00000000000000
  5.00000000000000    15.0000000000000    2.00000000000000    0.00000000000000
  0.00000000000000    2.00000000000000    10.0000000000000    0.00000000000000
  0.00000000000000    0.00000000000000    0.00000000000000    6.00000000000000 ]
b =
[ 4
  9
  2
  1 ]
//假如系统存在如上形式,求A*b
A = [0,0,40
     0,1,5
     1,0,5
     1,1,15
     1,2,2
     2,1,2
     2,2,10
     3,3,6];
b = [ 4
      9
      2
      1 ];
x = SparseMul(A, b)//回车得到如下结果
x =
[ 205.000000000000
  159.000000000000
  38.0000000000000
  6.00000000000000 ]