解对称变宽带线性方程组程序

! n: 方程组个数

! np: a数组的长度

! id: 存放方程组系数矩阵的第i行对角元在a数组的位置j，ID(i)=j

! b： 存放方程组右列向量

! ir: 计算结束标志，ir=0程序正常结束，ir=/0表示系数矩阵奇异-Solution of linear equations is ymmetrical broadband program