高斯消元_fortran

高斯消元___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