复数追赶法解三对角矩阵方程c c++源码 Complex chasing method solution three diagonal matrix equation c/c++ source code