Fortran Coder

查看: 6833|回复: 2
打印 上一主题 下一主题

[文件读写] 出现end of file during read!新手求助,这是书上的例题

[复制链接]

3

帖子

2

主题

0

精华

新人

F 币
27 元
贡献
12 点
跳转到指定楼层
楼主
发表于 2015-2-2 20:08:36 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
[Fortran] 纯文本查看 复制代码
 MODULE GAUSS
        CONTAINS
        SUBROUTINE SOLVE(A,B,X,N)
        IMPLICIT REAL*8(A-Z)
        INTEGER::I,K,N
        REAL*8::A(N,N),B(N),X(N)
        REAL*8::AUP(N,N),BUP(N)
        REAL*8::AB(N,N+1)
        AB(1:N,1:N)=A
        AB(:,N+1)=B
        DO K=1,N-1
          DO I=K+1,N
            TEMP=AB(I,K)/AB(K,K)
                AB(I,:)=AB(I,:)-TEMP*AB(K,:)
      END DO
    END DO
        AUP(:,:)=AB(1:N,1:N)
        BUP(:)=AB(:,N+1)
        CALL UPTRI(AUP,BUP,X,N)
        END SUBROUTINE SOLVE
        !!!!!!!!!!!!!!!!!!!!!!!
        subroutine uptri(a,b,x,n)
        IMPLICIT REAL*8(A-Z)
        INTEGER::I,J,N
        REAL*8::A(N,N),B(N),X(N)
        X(N)=B(N)/A(N,N)
        DO I=N-1,1,-1
          X(I)=B(I)
          DO J=I+1,N
          X(I)=X(I)-A(I,J)*X(J)
          END DO
          X(I)=X(I)/A(I,I)
    END DO
        END SUBROUTINE UPTRI
        END MODULE GAUSS
        PROGRAM MAIN
        USE GAUSS
        IMPLICIT REAL*8(A-Z)
        INTEGER,PARAMETER::N=4
        INTEGER::I,J
        REAL*8::A(N,N),B(N),X(N)
        OPEN(UNIT=11,FILE='FIN.TXT')
        
        OPEN(UNIT=12,FILE='FOUT.TXT')
        
    READ(11,*)
        READ(11,*)((A(I,J),J=1,N),I=1,N)
        READ(11,*)B
        CALL SOLVE(A,B,X,N)
        WRITE(12,101)X
        101 FORMAT(T5,'高斯消去法计算结果',/,T4,'X=',4(/F12.8))
        END PROGRAM MAIN
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

835

帖子

2

主题

0

精华

大宗师

F 币
3926 元
贡献
2334 点
沙发
发表于 2015-2-3 09:19:00 | 只看该作者
http://fcode.cn/guide-36-1.html
fortran常见错误,Q2003

3

帖子

2

主题

0

精华

新人

F 币
27 元
贡献
12 点
板凳
 楼主| 发表于 2015-2-3 22:23:15 | 只看该作者
谢谢,我看看。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-12-24 08:29

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表