亚当高斯消元法解线性方程组

最新推荐文章于 2022-04-22 22:31:44 发布

lafanda217

最新推荐文章于 2022-04-22 22:31:44 发布

阅读量3.4k

点赞数

分类专栏： Fortran 文章标签： file pivot subroutine io integer action

本文链接：https://blog.csdn.net/lafanda217/article/details/6448778

版权

该程序使用亚当高斯消元法从文件中读取线性方程组并求解。用户输入文件名，程序检查文件有效性，然后读取方程组系数和常数项。如果遇到无效数据，程序会给出提示。通过dsimul子程序进行矩阵操作，包括交换行、消除和归一化，以找到解决方案。

摘要由CSDN通过智能技术生成

PROGRAM hw
IMPLICIT NONE
INTEGER,PARAMETER::DBL=SELECTED_REAL_KIND(p=13)
REAL(KIND=DBL),ALLOCATABLE,DIMENSION(:,:)::a
REAL(KIND=DBL),ALLOCATABLE,DIMENSION(:)::b
REAL(KIND=DBL),ALLOCATABLE,DIMENSION(:)::soln

INTEGER::IO_state
CHARACTER(LEN=128)::IO_msg

INTEGER::error_flag
INTEGER::i,j,n
CHARACTER(LEN=128)::file_name

    ! FILE valid?
    WRITE(*,*) 'Enter the file name'
    DO
        READ(*,*) file_name
        OPEN(UNIT=1,FILE=file_name,IOSTAT=IO_state,IOMSG=IO_msg,STATUS='OLD',ACTION='READ')
        IF(IO_state==0) EXIT
        WRITE(*,*) IO_msg
        WRITE(*,*) 'Re-enter the file name:'
    END DO

    READ(1,*) n
    ALLOCATE(a(n,n),b(n),soln(n))
    DO i=1,n
        READ(1,*,IOSTAT=IO_state)(a(i,j),j=1,n), b(i)
        data_valid:IF(IO_state/