哪位高手给个反演的SIRT的算法~

#include #include #define M 5 #define N 5 void SIRT(float A[M][N],float D[M],float X[N]) { float R[M]={0.0},E[N]={0.0},S[M]={0.0},C[N],a=0,B=0; int i,j,k,n=0; 初始化矩阵 /*参量说明：A[M][N]线形方程组的系数矩阵；D[M]线形方程组等号右边的数值矩阵； X[N]解向量；R[M]残差矩阵；E[N]系数矩阵每列非零元素的个数矩阵；S[M]系数矩阵的每行元素的平方和矩阵； C[N]为第K+1次迭代解与第K次迭代解之差,用于解的精度控制.*/ for(i=0;i