出现end of file during read！新手求助，这是书上的例题

sxtsmile

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	MODULE GAUSS

02

CONTAINS

03	SUBROUTINE SOLVE(A,B,X,N)

04	IMPLICIT REAL*8(A-Z)

05	INTEGER::I,K,N

06	REAL*8::A(N,N),B(N),X(N)

07	REAL*8::AUP(N,N),BUP(N)

08	REAL*8::AB(N,N+1)

09	AB(1:N,1:N)=A

10	AB(:,N+1)=B

11	DO K=1,N-1

12	DO I=K+1,N

13	TEMP=AB(I,K)/AB(K,K)

14	AB(I,:)=AB(I,:)-TEMP*AB(K,:)

15

END DO

16

END DO

17	AUP(:,:)=AB(1:N,1:N)

18	BUP(:)=AB(:,N+1)

19	CALL UPTRI(AUP,BUP,X,N)

20	END SUBROUTINE SOLVE

21	!!!!!!!!!!!!!!!!!!!!!!!

22	subroutine uptri(a,b,x,n)

23	IMPLICIT REAL*8(A-Z)

24	INTEGER::I,J,N

25	REAL*8::A(N,N),B(N),X(N)

26	X(N)=B(N)/A(N,N)

27	DO I=N-1,1,-1

28

X(I)=B(I)

29	DO J=I+1,N

30	X(I)=X(I)-A(I,J)*X(J)

31

END DO

32	X(I)=X(I)/A(I,I)

33

END DO

34	END SUBROUTINE UPTRI

35	END MODULE GAUSS

36	PROGRAM MAIN

37

USE GAUSS

38	IMPLICIT REAL*8(A-Z)

39	INTEGER,PARAMETER::N=4

40	INTEGER::I,J

41	REAL*8::A(N,N),B(N),X(N)

42	OPEN(UNIT=11,FILE='FIN.TXT')

43

44	OPEN(UNIT=12,FILE='FOUT.TXT')

45

46	READ(11,*)

47	READ(11,*)((A(I,J),J=1,N),I=1,N)

48	READ(11,*)B

49	CALL SOLVE(A,B,X,N)

50	WRITE(12,101)X

51	101 FORMAT(T5,'高斯消去法计算结果',/,T4,'X=',4(/F12.8))

52	END PROGRAM MAIN

li913

http://fcode.cn/guide-36-1.html
fortran常见错误，Q2003

sxtsmile

谢谢，我看看。