c++解线性方程组的几种方法cpp

k=0,temp=0;k<n;k++)Р temp+=a[n][k]*a[k][j];?Р a[n][j]-=temp;Р }Р for(i=n+1;i<lenth;i++) //求第n+1层的下三角矩阵部分即LР {?for(k=0,temp=0;k<n;k++)Р temp+=a[i][k]*a[k][n];Р a[i][n]=(a[i][n]-temp)/a[n][n];Р }Р?}Р}Рint Doolittle_check(double temp_a[][Number],double temp_b[Number]) //若行列式不为零,将系数矩阵调整为顺序主子式大于零Р{Р?int i,j,k,maxi;double lik,temp;Р?for(k=0;k<lenth-1;k++)Р?{?Р j=k;Р for(maxi=i=k;i<lenth;i++)Р if(temp_a[i][j]>temp_a[maxi][j])?maxi=i;Р if(maxi!=k)?Р {?exchange_hang(k,maxi);Р for(j=0;j<lenth;j++)Р {?temp=temp_a[k][j];Р temp_a[k][j]=temp_a[maxi][j];Р temp_a[maxi][j]=temp;Р }Р }Р for(i=k+1;i<lenth;i++)Р {?Р lik=temp_a[i][k]/temp_a[k][k];Р for(j=k;j<lenth;j++)Р temp_a[i][j]=temp_a[i][j]-temp_a[k][j]*lik;Р temp_b[i]=temp_b[i]-temp_b[k]*lik;Р }Р?}Р?if(temp_a[lenth-1][lenth-1]==0) return 0; Р?return 1;Р}