高斯消元___fortran版本
DIMENSION A(10,10),B(10),X(10)
INTEGER N,K,I,J,MID,TEMP
WRITE(*,*)"输入矩阵的阶数:"
READ(*,*)N
WRITE(*,*)"请输入每行的系数:"
READ(*,*)((A(I,J),I=1,N),J=1,N)
WRITE(*,*)"请输入右侧的值"
READ(*,*)(B(I),I=1,N)
DO K=1,N
IF(A(K,K)==0)then
WRITE(*,*)"无解"
ENDIF
DO I=K+1,N
MID=A(I,k)/A(k,K)
DO J=K,N
A(I,J)=A(I,J)-MID*A(K,J)
ENDDO
B(I)=B(I)-MID*B(K)
ENDDO
ENDDO
DO I=1,N
DO J=1,N
WRITE(*,*)A(I,J)
ENDDO
WRITE(*,*)"\"
ENDDO
DO I=N,1,-1
TEMP=B(I)
DO J=N,I+1,-1
TEMP=TEMP-A(I,J)
ENDDO
X(I)=TEMP/A(I,I)
ENDDO
DO I=1,N
WRITE(*,*)X(I)
ENDDO
END