Go to Qiita Advent Calendar 2024 Top
LoginSignup
4
5

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?

FortranAdvent Calendar 2019

Day 18この記事誰得？ 私しか得しないニッチな技術で記事投稿！
@kaizen_nagoya(Dr. Kiyoshi Ogawa)

久々のFORTRANで出したエラー群（未解決有）docker(112) coding(145)

Last updated at Posted at 2021-03-27

41Hb9INzhYL.SL500.jpg
演習 大学院入試問題[数学]II
姫野 俊一
https://bookmeter.com/books/233806

のFORTRANプログラムを動かそうとして出たエラー群。

この記事は、fortranに限定して記載し、
fortranに再び馴染むために、打ち間違いなどを記録して、どういうコンパイルエラーがでるかを
確認するために記録しようとしていました。
論理は別記事にする予定で、演算子などあまり深く考えずに入力していました。

使った資料は第一版で、文字形(font)が読みにくい上に、印刷がかすれていて、文字が判別しにくく、コンパイルエラーが出ないと誤植か、誤入力かの判定ができない状態。
30年前の印刷の書籍で、文字形も1とIが判定しずらく、KとRも見間違え。
-と=は単純な打ち間違い（シフトの押すタイミングがずれて）。
まだ、もうすこしプログラムがあり、順次追加します。
統計のプログラム全部打ち終えたら、統計の記事を別途書く予定です。
１日、１プログラムしか入力できそうになく、時間がかかりますが、よろしくお願いします。

docker

使ったのは、

$ docker run -v /tmp:/tmp -it kaizenjapan/fortran /bin/bash

使う前に更新。

# apt update; apt -y upgrade

pushしてある。
/home/fortran
にソースをコピーしておいた。入っているはず。

最初にdockerでFORTRANについて利用した記事は下記。

プログラムちょい替え(11)Fortran 多次元配列のアクセス順序による計算時間の違い: docker(96)
https://qiita.com/kaizen_nagoya/items/3d657649d74fdd753bad

p487.f

bash
# gfortran p487.f 
p487.f:12:11:

   C(N,q)=N
           1
Warning: Legacy Extension: REAL array index at (1)
/tmp/ccq2t893.o: In function `insatu_':
p487.f:(.text+0xb2): undefined reference to `r_'
collect2: error: ld returned 1 exit status

ソースコードは下記。

p487.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806
C	PASKAL
	CALL PASKAL
	CALL INSATU
	STOP
	END
	
	SUBROUTINE PASKAL
	INTEGER C, R
	COMMON /NCR/C(20,20)
	NMAX=20
	DO 1 N=1, NMAX
		C(N,q)=N
		C(N,N)=1
1	CONTINUE
	DO 2 N=2, NMAX
		DO 3 R=2, N-1
			C(N,R)=C(N-1,R-1)*C(N-1,R)
3		CONTINUE
2	CONTINUE
	RETURN
	END
	
	SUBROUTINE INSATU
	COMMON /NCR/K(20,20)
	DO 1 I=1,10
		WRITE(6,601) (R(I,J),J=1,I)
601		FORMAT(' 1',I3,I4,I5,10I7,I5,I5,I5,I4,I3,I2)
1	CONTINUE
	RETURN
	END

単純な打ち誤り。

# gfortran p487.f 
/tmp/cc5XGojg.o: In function `insatu_':
p487.f:(.text+0xb2): undefined reference to `r_'
collect2: error: ld returned 1 exit status

1017を10I7に変えても同じエラー。

10I7をI7に変えても同じエラー。

K(20,20)でなくてR(20,20)だった。

直してコンパイルして実行。

./a.out
At line 26 of file p487.f (unit = 6, file = 'stdout')
Fortran runtime error: Expected INTEGER for item 1 in formatted transfer, got REAL
(' 1',I3,I4,I5,I7,I5,I5,I5,I4,I3,I2)
      ^

Error termination. Backtrace:
#0  0x7f636cbae32a
#1  0x7f636cbaeed5
#2  0x7f636cbaf69d
#3  0x7f636cd18cb1
#4  0x7f636cd25522
#5  0x7f636cd28222
#6  0x7f636cd286dc
#7  0x55911f8e2935
#8  0x55911f8e2ae4
#9  0x55911f8e2b29
#10  0x7f636c7c3bf6
#11  0x55911f8e2769
#12  0xffffffffffffffff

N=2, ではなく、N=3,だった。

直しても同じエラー。

I3じゃなくてR3にしてみる。

# gfortran p487.f 
p487.f:27:22:

 601  FORMAT(' 1', R3, I4, I5, I7, I5, I5, I5, I4, I3, I2)
                      1
Error: Unexpected element '3' in format string at (1)
p487.f:26:72:

   WRITE(6,601) (R(I,J),J=1,I)
                                                                        1
Error: FORMAT label 601 at (1) not defined

エラーが取れない。
ひとまず投稿してから考える。

複数のR()をK()に変えて、コンパイル。実行。

bash
# gfortran p487.f
# ./a.out
 1  1
 1  2   1
 1  3   2    1
 1  4   6    2      1
 1  5  24   12      2    1
 1  6 120  288     24    2    1
 1  7 72034560   6912   48    2    1
 1  85040*****************   96    2   1
 1  9**************************  192   2  1
 1 10*********      0    0    0***** 384  2 1

あと一息。@cure_honey さんの指摘に基づき*を+に変更。

# gfortran p487.F
root@0546c6fb4c77:/tmp# ./a.out
 1  1
 1  2   1
 1  3   3    1
 1  4   6    4      1
 1  5  10   10      5      1
 1  6  15   20     15      6      1
 1  7  21   35     35     21      7      1
 1  8  28   56     70     56     28      8      1
 1  9  36   84    126    126     84     36      9      1
 1 10  45  120    210    252    210    120     45     10      1

現状のソースは

p487.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806  
C	PASKAL
	CALL PASKAL
	CALL INSATU
	STOP
	END
	
	SUBROUTINE PASKAL
	INTEGER C, R
	COMMON /NCR/C(20,20)
	NMAX=20
	DO 1 N=1, NMAX
		C(N,1)=N
		C(N,N)=1
1	CONTINUE
	DO 2 N=3, NMAX
		DO 3 R=2, N-1
			C(N,R)=C(N-1,R-1)+C(N-1,R)
3		CONTINUE
2	CONTINUE
	RETURN
	END
	
	SUBROUTINE INSATU
	COMMON /NCR/K(20,20)
	DO 1 I=1,10
		WRITE(6,601) (K(I,J),J=1,I)
601		FORMAT(' 1', I3, I4, I5, 10I7, I5, I5, I5, I4, I3, I2)
1	CONTINUE
	RETURN
	END

p487-2.f

# gfortran p487-2.f 
p487-2.f:2:19:

 50 READA(5,100) N
                   1
Error: Expected comma in I/O list at (1)
p487-2.f:3:5:

 100 FORMAT(I3)
     1
Error: Non-numeric character in statement label at (1)
p487-2.f:3:5:

 100 FORMAT(I3)
     1
Error: Unclassifiable statement at (1)
p487-2.f:5:5:

 200 FORMAT(' ',5X,'N',7X,'KAIJOU',5X,'STIRLING')
     1
Error: Non-numeric character in statement label at (1)
p487-2.f:5:5:

 200 FORMAT(' ',5X,'N',7X,'KAIJOU',5X,'STIRLING')
     1
Error: Unclassifiable statement at (1)
p487-2.f:13:5:

 201 FORMAT(' ',I6,I13,F13,0)
     1
Error: Non-numeric character in statement label at (1)
p487-2.f:13:5:

 201 FORMAT(' ',I6,I13,F13,0)
     1
Error: Unclassifiable statement at (1)
p487-2.f:4:72:

  WRITE(6,200)
                                                                        1
Error: FORMAT label 200 at (1) not defined
p487-2.f:12:72:

  WRITE(6,201)N,KAI,STIRL
                                                                        1
Error: FORMAT label 201 at (1) not defined
p487-2.f:14:72:

  GOTO 50
                                                                        1
Error: Label 50 referenced at (1) is never defined

ソースは

p487.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806    
C	FACTIORIAL
50	READA(5,100) N
100 FORMAT(I3)
	WRITE(6,200)
200 FORMAT(' ',5X,'N',7X,'KAIJOU',5X,'STIRLING')
	IF(N.LE.0) STOP
	KAI=1
	DO 10 I=1,N
10	KAI=KAI*I
	AN=A
	STIRL=SQRT(2.*3.14159*AN)*(AN/2.71828)**N
	WRITE(6,201)N,KAI,STIRL
201 FORMAT(' ',I6,I13,F13,0)
	GOTO 50
	STOP
	END

50 READA(5,100) N
と本には,がない。
50 READA(5,100) , N
にしないとgfortranはコンパイルが通らないのかも。

201 FORMAT(' ',I6,I13,F13,0)
は,を.に訂正（誤植）。
201 FORMAT(' ',I6,I13,F13.0)

AN=A
のAをNに訂正（誤植）。
AN=N

再度コンパイル

# gfortran p487-2.f 
p487-2.f:4:11:

 50  READA(5,100) , N
           1
Error: FORMAT tag at (1) must be of type default-kind CHARACTER or of INTEGER
#

この時点のソース

p487-2.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806  
C		FACTIORIAL
50		READA(5,100) , N
100 	FORMAT(I3)
		WRITE(6,200)
200 	FORMAT(' ',5X,'N',7X,'KAIJOU',5X,'STIRLING')
		IF(N.LE.0) STOP
		KAI=1
		DO 10 I=1,N
10			KAI=KAI*I

		STIRL=SQRT(2.*3.14159*AN)*(AN/2.71828)**N
		WRITE(6,201)N,KAI,STIRL
201 	FORMAT(' ',I6,I13,F13.0)
		GOTO 50
		STOP
		END

50 READA(5,100),N
READAじゃなくてREADだった。
50 READ(5,100)N

bash
# gfortran p487-2.f 
# ./a.out
10
      N       KAIJOU     STIRLING
     10      3628800          NaN

現在のソース

p487-2.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806  
C		FACTIORIAL
 50		READ(5,100)N
 100 	FORMAT(I3)
		WRITE(6,200)
 200 	FORMAT(' ',5X,'N',7X,'KAIJOU',5X,'STIRLING')
		IF(N.LE.0) STOP
		KAI=1
		DO 10 I=1,N
 10			KAI=KAI*I
		STIRL=SQRT(2.*3.14159*AN)*(AN/2.71828)**N
		WRITE(6,201)N,KAI,STIRL
 201 	FORMAT(' ',I6,I13,F13.0)
		GOTO 50
		STOP
		END

#p491.f

# gfortran p491.f 
p487-2.f:4:11:

 50  READA(5,100) , N
           1
Error: FORMAT tag at (1) must be of type default-kind CHARACTER or of INTEGER
root@0546c6fb4c77:/tmp# gfortran p491.f 
p491.f:7:19:

  90  FORMAT(' ',　’ーMR COURSE-'）
                   1
Error: Unexpected element '\xe3' in format string at (1)
p491.f:9:6:

  100 FORMT(' ','I=DOBOK', 1X,'2=KIKAI',1X,'3=DENKI',1X,
      1
Error: Unclassifiable statement at (1)
p491.f:10:9:

     +'4=OBUTU',1X,'5=GENSI',1X,'6=JOUHO')
         1
Error: Unclassifiable statement at (1)
p491.f:15:7:

  20  FIRNAT(I5,F10.2)
       1
Error: Unclassifiable statement at (1)
p491.f:17:7:

   DO 2 I-1,N
       1
Error: Unclassifiable statement at (1)
p491.f:24:8:

    S-S+HH(I)
        1
Error: Unclassifiable statement at (1)
p491.f:28:7:

   DO 12 I-1,N
       1
Error: Unclassifiable statement at (1)
p491.f:29:24:

    SS=(HH(I)-HEIKIN(N)**2+SS
                        1
Error: Expected a right parenthesis in expression at (1)
p491.f:6:72:

   WRITE(6,90)
                                                                        1
Error: FORMAT label 90 at (1) not defined
p491.f:8:72:

   WRITE(6,100)
                                                                        1
Error: FORMAT label 100 at (1) not defined
p491.f:14:72:

    READ(5,20)LL(I),HH(I)
                                                                        1
Error: FORMAT label 20 at (1) not defined
p491.f:23:12:

   DO 11 U=1,N
            1
Warning: Deleted feature: Loop variable at (1) must be integer

最初の入力は

p491.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806  
C		URIAGE
C		HEIKIN, HYOUJUN HENSA, SAIDAI, SAISHOU
		INTEGER LL(200)
		REAL HH(200), S, SS, HEIKIN, BUNSAN, SIGMA
		INTEGER N, I
		WRITE(6,90)
 90		FORMAT(' ',　’ーMR COURSE-'）
		WRITE(6,100)
 100	FORMT(' ','I=DOBOK', 1X,'2=KIKAI',1X,'3=DENKI',1X,
	   +'4=OBUTU',1X,'5=GENSI',1X,'6=JOUHO')
	   READ(5,10)N
 10		FORMAT(I5)
		DO 1 I=1,N
			READ(5,20)LL(I),HH(I)
 20		FIRNAT(I5,F10.2)
 1		CONTINUE
		DO 2 I-1,N
			WRITE(6,30)LL(I),HH(I)
30			FORMAT('  SENKOU',I4,4X,'WARIAI'F6.2)
 2		CONTINUE
C		HEIKIN,HYOUJUNHENSA
		S=0.0
		DO 11 U=1,N
			S-S+HH(I)
 11		CONTINUE
 		HEIKIN=S/REAL(N)
		 SS=0.0
		DO 12 I-1,N
			SS=(HH(I)-HEIKIN(N)**2+SS
 12		CONTINUE
 		BUNSAN=SS/REAL(N)
		 SIGMA=SQRT(BUNSAN)
		WRITE(6,40)HEIKIN
 40		FORMAT(' HEININ0',7X,F4.2)
 		WRITE(6,50)SIGMA
 50		FORMAT(' HYOUJUN HENSA=',F4.2)
C		saidai, saishou
		KMIN=1
		KMAX=1
		DO 22 I=2,N
			IF(HH(I).LT.HH(KMIN))KMIN=1
			IF(HH(I).GT.HH(KMAX))KMAX=1
 22		CONTINUE
		WRITE(6,60)HH(KMIN),LL(KMIN)
 60		FORMAT(' MIN=',4X,F10.2,3X,'SENKOU=',I5)
 		WRITE(6,80)HH(KMAX),LL(KMAX)
 80		FORMAT(' MAX=',4X,F10.2,3X,'SENKOU',I5)
 		STOP
		END

継続行の+印の前の空白文字の数

 +'4=OBUTU',1X,'5=GENSI',1X,'6=JOUHO')

+の前の空白が正確に５文字必要で、タブで代替すると今のところうまく行っていない。

２バイト文字混入

90 FORMAT(' ',　’ーMR COURSE-'）
２バイト文字が入っていた。
90 FORMAT(' ',　’-MR COURSE-'）

予約語(FORMAT)打ち間違い

100 FORMT(' ','I=DOBOK', 1X,'2=KIKAI',1X,'3=DENKI',1X,
は
100 FORMAT(' ','I=DOBOK', 1X,'2=KIKAI',1X,'3=DENKI',1X,

20 FIRNAT(I5,F10.2)
は
20 FORMAT(I5,F10.2)

変数間違い

DO 11 U=1,N
は、
DO 11 I=1,N

演算子間違い

	DO 2 I-1,N

は
DO 2 I=1,N

	DO 12 I-1,N

は
DO 12 I=1,N

		S-S+HH(I)

は
S=S+HH(I)

その他の打ち間違い

		SS=(HH(I)-HEIKIN(N)**2+SS

は
SS=(HH(I)-HEIKIN)**2+SS

bash
# gfortran p491.f 
# ./a.out
 -MR course-
 I=DOBOK 2=KIKAI 3=DENKI 4=OBUTU 5=GENSI 6=JOUHO
6
5
4
3
2
1
0
  SENKOU   5    WARIAI  0.00
  SENKOU   4    WARIAI  0.00
  SENKOU   3    WARIAI  0.00
  SENKOU   2    WARIAI  0.00
  SENKOU   1    WARIAI  0.00
  SENKOU   0    WARIAI  0.00
 HEININ0       0.00
 HYOUJUN HENSA=0.00
 MIN=          0.00   SENKOU=    5
 MAX=          0.00   SENKOU    5

文字の誤植を直し、コンパイル、実行。

bash
# fortran p491.f
# ./a.out
 -MR COURSE-
 I=DOBOK 2=KIKAI 3=DENKI 4=OBUTU 5=GENSI 6=JOUHO
6
1	0.25
At line 14 of file p491.F (unit = 5, file = 'stdin')
Fortran runtime error: Bad value during integer read

Error termination. Backtrace:
#0  0x7f3fb5daa32a
#1  0x7f3fb5daaed5
#2  0x7f3fb5dab69d
#3  0x7f3fb5f1f34c
#4  0x7f3fb5f2339c
#5  0x7f3fb5f2472c
#6  0x564623933bf9
#7  0x564623934164
#8  0x7f3fb59bfbf6
#9  0x564623933929
#10  0xffffffffffffffff

入力は、最初に何個値を入れるか。
次に、２つの数値の値を順に入れるはず。

p491.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806  
C		URIAGE
C		HEIKIN, HYOUJUN HENSA, SAIDAI, SAISHOU
		INTEGER LL(200)
		REAL HH(200), S, SS, HEIKIN, BUNSAN, SIGMA
		INTEGER N, I
		WRITE(6,90)
90		FORmat(' ','-MR COURSE-')
		WRITE(6,100)
100		FORMAT(' ','I=DOBOK', 1X,'2=KIKAI',1X,'3=DENKI',1X,
     +'4=OBUTU',1X,'5=GENSI',1X,'6=JOUHO')
	   READ(5,10)N
10		FORMAT(I5)
		DO 1 I=1,N
			READ(5,20)LL(I),HH(I)
20		FORMAT(I5,F10.2)
1		CONTINUE
		DO 2 I=1,N
			WRITE(6,30)LL(I),HH(I)
30			FORMAT('  SENKOU',I4,4X,'WARIAI'F6.2)
2		CONTINUE
C		HEIKIN,HYOUJUNHENSA
		S=0.0
		DO 11 I=1,N
			S=S+HH(I)
11		CONTINUE
 		HEIKIN=S/REAL(N)
		 SS=0.0
		DO 12 I=1,N
			SS=(HH(I)-HEIKIN)**2+SS
12		CONTINUE
 		BUNSAN=SS/REAL(N)
		 SIGMA=SQRT(BUNSAN)
		WRITE(6,40)HEIKIN
40		FORMAT(' HEININ0',7X,F4.2)
 		WRITE(6,50)SIGMA
50		FORMAT(' HYOUJUN HENSA=',F4.2)
C		saidai, saishou
		KMIN=1
		KMAX=1
		DO 22 I=2,N
			IF(HH(I).LT.HH(KMIN))KMIN=1
			IF(HH(I).GT.HH(KMAX))KMAX=1
 22		CONTINUE
		WRITE(6,60)HH(KMIN),LL(KMIN)
 60		FORMAT(' MIN=',4X,F10.2,3X,'SENKOU=',I5)
 		WRITE(6,80)HH(KMAX),LL(KMAX)
 80		FORMAT(' MAX=',4X,F10.2,3X,'SENKOU=',I5)
 		STOP
		END

p491-2.f

# gfortran p491-2.f
p491-2.f:4:17:

   READ(5,100）
                 1
Error: Syntax error in READ statement at (1)
p491-2.f:26:7:

   DO 10 I-1,N
       1
Error: Unclassifiable statement at (1)
p491-2.f:31:23:

    SYY=SYY+YY(I)*:2
                       1
Error: Invalid character in name at (1)
p491-2.f:37:7:

     1 /(SQRT(SXX-FN*XMEAN**2)
       1
Error: Invalid character in name at (1)
p491-2.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806  
C		correlation coefficient
C		MAIN PROGRAM
		DIMENSION X(99),Y(99)
		READ(5,100）
100		FORMAT(I5)
		READ(5,110)(X(I),Y(I),I=1,N)
110		FORMAT(2F10.3)
		WRITE(6,105)
105		FORMAT(' OCCUPIED RATE',3X,'UNEXAMINED ENTRANT RATE')
		WRITE(6,120)(X(I),Y(I),I=1,N)
120		FORMAT(F10.3,4X,F10.3)
		R=SOKAN(N,X,Y)
		WRITE(6,200)R
200 	FORMAT(' ','SOUKAN KEISUU=',F10.4)
		STOP
		END
		
C 		SUBPROGRAM
		FUNCTION SOUKAN(N,XX,YY)
		DIMENSION XX(99),YY(99)
		SX=0.0
		SY=0.0
		SXY=0.0
		SXX=0.0
		SYY=0.0
		DO 10 I-1,N
			SX=SX+XX(I)
			SY=SY+YY(I)
			SXY=SXY+XX(I)*YY(I)
			SXX=SXX+XX(I)**2
			SYY=SYY+YY(I)*:2
10		CONTINUE
		FN=N
		XMEAN=SX/FN
		YMEAN=SY/FN
		SOKAN=(SXY-FN*XMEAN*YMEAN)
    1	/(SQRT(SXX-FN*XMEAN**2)
	2	*SQRT(SYY-FN*YMEAN**2))
	 RETURN
	 END

p502.f

# gfortran p502.f
p502.f:3:7:

   DIMENSION X(100),(100)
       1
Error: Unclassifiable statement at (1)
p502.f:6:7:

    DO 1 I-1,N
       1
Error: Unclassifiable statement at (1)
p502.f:7:19:

    READ(5,20)X(I),Y(I)
                   1
Error: Syntax error in READ statement at (1)
p502.f:9:12:

   WRITE'6,22)(I,X(I),Y(I),I-1,N)
            1
Error: Syntax error in WRITE statement at (1)
p502.f:21:7:

   DO 3 I-1,N
       1
Error: Unclassifiable statement at (1)

D(5,20)X(I),Y(I)
                   1
Error: Syntax error in READ statement at (1)
p502.f:9:12:

   WRITE'6,22)(I,X(I),Y(I),I-1,N)
            1
Error: Syntax error in WRITE statement at (1)
p502.f:21:7:

最初のソースは

p502.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806  
C		LEAST SQUARE, Y=AX+B
C		NUMBER OF MR YS NUMBER OF DR
		DIMENSION X(100),(100)
		READ(5,10)N
 10		FORMAT(I5)
 		DO 1 I-1,N
			READ(5,20)X(I),Y(I)
 20		FORMAT(' ','DEPT.',3X,'NO.OF MR',3X,'NO.OF DR')
		WRITE'6,22)(I,X(I),Y(I),I-1,N)
 22		FORMAT(I5,3X,F8.0,3X,F8.0)
		SX=0.0
		SY=0.0
		DO 2 I=1,N
			SX=SX+X(I)
			SY=SY+Y(I)
 2		CONTINUE
		XM=SX/REAL(N)
		YM=SY/REAL(N)	
		S1=0.0
		S2=0.0 
		DO 3 I-1,N
			S1=S1+(X(I)-XM)*(Y(I)-YM)
			S2=S2+(X(I)-XM)**2
 3		CONTINUE
 		A=S1/S2
		B=YM-A*XM
		WRITE(6,30)A
 30		FORMAT(' ','A=',F12.5)
 		WRITE(6,40)B
 40		FORMAT(' ','B=',F12.5)
 		STOP
		END
bash
./a.out
9
34.	11.
At line 7 of file p502.f (unit = 5, file = 'stdin')
Fortran runtime error: Bad value during floating point read

Error termination. Backtrace:
#0  0x7f595bf1332a
#1  0x7f595bf13ed5
#2  0x7f595bf1469d
#3  0x7f595c088d6d
#4  0x7f595c08c0a8
#5  0x7f595c08d72c
#6  0x555aa5d17b31
#7  0x555aa5d17fb1
#8  0x7f595bb28bf6
#9  0x555aa5d17929
#10  0xffffffffffffffff

現状のソース

p502.f
C       演習 大学院入試問題[数学]II姫野 俊一 p.487
C       https://bookmeter.com/books/233806  

C		LEAST SQUARE, Y=AX+B
C		NUMBER OF MR VS NUMBER OF DR
		DIMENSION X(100),Y(100)
		READ(5,10)N
 10		FORMAT(I5)
 		DO 1 I=1,N
			READ(5,20)X(I),Y(I)
 20		FORMAT(2F10.0)
 1		CONTINUE
 		WRITE(6,21)
 21		FORMAT(' ','DEPT.',3X,'NO.OF MR',3X,'NO.OF DR')
		WRITE(6,22)(I,X(I),Y(I),I=1,N)
 22		FORMAT(I5,3X,F8.0,3X,F8.0)
		SX=0.0
		SY=0.0
		DO 2 I=1,N
			SX=SX+X(I)
			SY=SY+Y(I)
 2		CONTINUE
		XM=SX/REAL(N)
		YM=SY/REAL(N)	
		S1=0.0
		S2=0.0 
		DO 3 I=1,N
			S1=S1+(X(I)-XM)*(Y(I)-YM)
			S2=S2+(X(I)-XM)**2
 3		CONTINUE
 		A=S1/S2
		B=YM-A*XM
		WRITE(6,30)A
 30		FORMAT(' ','A=',F12.5)
 		WRITE(6,40)B
 40		FORMAT(' ','B=',F12.5)
 		STOP
		END

関連資料

' @kazuo_reve 私が効果を確認した「小川メソッド」
https://qiita.com/kazuo_reve/items/a3ea1d9171deeccc04da

' @kazuo_reve 新人の方によく展開している有益な情報
https://qiita.com/kazuo_reve/items/d1a3f0ee48e24bba38f1

' @kazuo_reve Vモデルについて勘違いしていたと思ったこと
https://qiita.com/kazuo_reve/items/46fddb094563bd9b2e1e

Engineering Festa 2024前に必読記事一覧

登壇直後版 色使い(JIS安全色) Qiita Engineer Festa 2023〜私しか得しないニッチな技術でLT〜 スライド編 0.15
https://qiita.com/kaizen_nagoya/items/f0d3070d839f4f735b2b

プログラマが知っていると良い「公序良俗」
https://qiita.com/kaizen_nagoya/items/9fe7c0dfac2fbd77a945

逆も真：社会人が最初に確かめるとよいこと。OSEK(69)、Ethernet(59)
https://qiita.com/kaizen_nagoya/items/39afe4a728a31b903ddc

統計の嘘。仮説（127）
https://qiita.com/kaizen_nagoya/items/63b48ecf258a3471c51b

自分の言葉だけで論理展開できるのが天才なら、文章の引用だけで論理展開できるのが秀才だ。仮説（136）
https://qiita.com/kaizen_nagoya/items/97cf07b9e24f860624dd

参考文献駆動執筆(references driven writing)・デンソークリエイト編
https://qiita.com/kaizen_nagoya/items/b27b3f58b8bf265a5cd1

「何を」よりも「誰を」。１０年後のために今見習いたい人たち
https://qiita.com/kaizen_nagoya/items/8045978b16eb49d572b2

Qiitaの記事に３段階または５段階で到達するための方法
https://qiita.com/kaizen_nagoya/items/6e9298296852325adc5e

出力(output)と呼ばないで。これは状態(state)です。
https://qiita.com/kaizen_nagoya/items/80b8b5913b2748867840

coding (101) 一覧を作成し始めた。omake:最近のQiitaで表示しない5つの事象
https://qiita.com/kaizen_nagoya/items/20667f09f19598aedb68

あなたは「勘違いまとめ」から、勘違いだと言っていることが勘違いだといくつ見つけられますか。人間の間違い(human error(125))の種類と対策
https://qiita.com/kaizen_nagoya/items/ae391b77fffb098b8fb4

プログラマの「プログラムが書ける」思い込みは強みだ。３つの理由。仮説（168）統計と確率(17) , OSEK(79)
https://qiita.com/kaizen_nagoya/items/bc5dd86e414de402ec29

出力(output)と呼ばないで。これは状態(state)です。
https://qiita.com/kaizen_nagoya/items/80b8b5913b2748867840

これからの情報伝達手段の在り方について考えてみよう。炎上と便乗。
https://qiita.com/kaizen_nagoya/items/71a09077ac195214f0db

ISO/IEC JTC1 SC7 Software and System Engineering
https://qiita.com/kaizen_nagoya/items/48b43f0f6976a078d907

アクセシビリティの知見を発信しよう！（再び）
https://qiita.com/kaizen_nagoya/items/03457eb9ee74105ee618

統計論及確率論輪講（再び）
https://qiita.com/kaizen_nagoya/items/590874ccfca988e85ea3

読者の心をグッと惹き寄せる7つの魔法
https://qiita.com/kaizen_nagoya/items/b1b5e89bd5c0a211d862

「@kazuo_reve 新人の方によく展開している有益な情報」確認一覧
https://qiita.com/kaizen_nagoya/items/b9380888d1e5a042646b

ソースコードで議論しよう。日本語で議論するの止めましょう（あるプログラミング技術の議論報告）
https://qiita.com/kaizen_nagoya/items/8b9811c80f3338c6c0b0

脳内コンパイラの3つの危険
https://qiita.com/kaizen_nagoya/items/7025cf2d7bd9f276e382

心理学の本を読むよりはコンパイラ書いた方がよくね。仮説（34）
https://qiita.com/kaizen_nagoya/items/fa715732cc148e48880e

NASAを超えるつもりがあれば読んでください。
https://qiita.com/kaizen_nagoya/items/e81669f9cb53109157f6

データサイエンティストの気づき!「勉強して仕事に役立てない人。大嫌い!!」『それ自分かも?』ってなった!!!
https://qiita.com/kaizen_nagoya/items/d85830d58d8dd7f71d07

「ぼくの好きな先生」「人がやらないことをやれ」プログラマになるまで。仮説（37）　
https://qiita.com/kaizen_nagoya/items/53e4bded9fe5f724b3c4

なぜ経済学徒を辞め、計算機屋になったか（経済学部入学前・入学後・卒業後対応） 転職(1)
https://qiita.com/kaizen_nagoya/items/06335a1d24c099733f64

プログラミング言語教育のXYZ。 仮説（52）
https://qiita.com/kaizen_nagoya/items/1950c5810fb5c0b07be4

【24卒向け】9ヶ月後に年収１０００万円を目指す。二つの関門と三つの道。
https://qiita.com/kaizen_nagoya/items/fb5bff147193f726ad25

「【25卒向け】Qiita Career Meetup for STUDENT」予習の勧め
https://qiita.com/kaizen_nagoya/items/00eadb8a6e738cb6336f

大学入試不合格でも筆記試験のない大学に入って卒業できる。卒業しなくても博士になれる。
https://qiita.com/kaizen_nagoya/items/74adec99f396d64b5fd5

全世界の不登校の子供たち「博士論文」を書こう。世界子供博士論文遠隔実践中心 安全(99)
https://qiita.com/kaizen_nagoya/items/912d69032c012bcc84f2

小川メソッド　覚え（書きかけ）
https://qiita.com/kaizen_nagoya/items/3593d72eca551742df68

DoCAP（ドゥーキャップ）って何ですか？
https://qiita.com/kaizen_nagoya/items/47e0e6509ab792c43327

views 20,000越え自己記事一覧
https://qiita.com/kaizen_nagoya/items/58e8bd6450957cdecd81

Views１万越え、もうすぐ１万記事一覧 最近いいねをいただいた213記事
https://qiita.com/kaizen_nagoya/items/d2b805717a92459ce853

自己記事一覧

Qiitaで逆リンクを表示しなくなったような気がする。時々、スマフォで表示するとあらわっることがあり、完全に削除したのではなさそう。

４月以降、せっせとリンクリストを作り、統計を取って確率を説明しようとしている。
2025年２月末を目標にしている。

物理記事　上位100
https://qiita.com/kaizen_nagoya/items/66e90fe31fbe3facc6ff

量子(0) 計算機, 量子力学
https://qiita.com/kaizen_nagoya/items/1cd954cb0eed92879fd4

数学関連記事１００
https://qiita.com/kaizen_nagoya/items/d8dadb49a6397e854c6d

統計(0)一覧
https://qiita.com/kaizen_nagoya/items/80d3b221807e53e88aba

図(0) state, sequence and timing. UML and お絵描き
https://qiita.com/kaizen_nagoya/items/60440a882146aeee9e8f

品質一覧
https://qiita.com/kaizen_nagoya/items/2b99b8e9db6d94b2e971

言語・文学記事　１００
https://qiita.com/kaizen_nagoya/items/42d58d5ef7fb53c407d6

医工連携関連記事一覧
https://qiita.com/kaizen_nagoya/items/6ab51c12ba51bc260a82

自動車　記事　１００
https://qiita.com/kaizen_nagoya/items/f7f0b9ab36569ad409c5

通信記事１００
https://qiita.com/kaizen_nagoya/items/1d67de5e1cd207b05ef7

日本語（０）一欄
https://qiita.com/kaizen_nagoya/items/7498dcfa3a9ba7fd1e68

英語(0) 一覧
https://qiita.com/kaizen_nagoya/items/680e3f5cbf9430486c7d

転職(0)一覧
https://qiita.com/kaizen_nagoya/items/f77520d378d33451d6fe

仮説（0）一覧（目標100現在40）
https://qiita.com/kaizen_nagoya/items/f000506fe1837b3590df

音楽　一覧(0)
https://qiita.com/kaizen_nagoya/items/b6e5f42bbfe3bbe40f5d

「@kazuo_reve 新人の方によく展開している有益な情報」確認一覧
https://qiita.com/kaizen_nagoya/items/b9380888d1e5a042646b

Qiita(0)Qiita関連記事一覧（自分）
https://qiita.com/kaizen_nagoya/items/58db5fbf036b28e9dfa6

鉄道（０）鉄道のシステム考察はてっちゃんがてつだってくれる
https://qiita.com/kaizen_nagoya/items/26bda595f341a27901a0

安全（0）安全工学シンポジウムに向けて: 21
https://qiita.com/kaizen_nagoya/items/c5d78f3def8195cb2409

一覧の一覧( The directory of directories of mine.) Qiita(100)
https://qiita.com/kaizen_nagoya/items/7eb0e006543886138f39

Ethernet 記事一覧　Ethernet(0)
https://qiita.com/kaizen_nagoya/items/88d35e99f74aefc98794

Wireshark 一覧 wireshark(0)、Ethernet(48)
https://qiita.com/kaizen_nagoya/items/fbed841f61875c4731d0

線網（Wi-Fi）空中線(antenna)(0) 記事一覧(118/300目標)
https://qiita.com/kaizen_nagoya/items/5e5464ac2b24bd4cd001

OSEK OS設計の基礎　OSEK(100)
https://qiita.com/kaizen_nagoya/items/7528a22a14242d2d58a3

Error一覧 error(0)
https://qiita.com/kaizen_nagoya/items/48b6cbc8d68eae2c42b8

++ Support(0)　
https://qiita.com/kaizen_nagoya/items/8720d26f762369a80514

Coding(0) Rules, C, Secure, MISRA and so on
https://qiita.com/kaizen_nagoya/items/400725644a8a0e90fbb0

coding (101) 一覧を作成し始めた。omake:最近のQiitaで表示しない5つの事象
https://qiita.com/kaizen_nagoya/items/20667f09f19598aedb68

プログラマによる、プログラマのための、統計(0)と確率のプログラミングとその後
https://qiita.com/kaizen_nagoya/items/6e9897eb641268766909

なぜdockerで機械学習するか 書籍・ソース一覧作成中 (目標100)
https://qiita.com/kaizen_nagoya/items/ddd12477544bf5ba85e2

言語処理100本ノックをdockerで。python覚えるのに最適。:10+12
https://qiita.com/kaizen_nagoya/items/7e7eb7c543e0c18438c4

プログラムちょい替え（0）一覧:4件
https://qiita.com/kaizen_nagoya/items/296d87ef4bfd516bc394

Python(0)記事をまとめたい。
https://qiita.com/kaizen_nagoya/items/088c57d70ab6904ebb53

官公庁・学校・公的団体（NPOを含む）システムの課題、官（０）
https://qiita.com/kaizen_nagoya/items/04ee6eaf7ec13d3af4c3

「はじめての」シリーズ　 ベクタージャパン　
https://qiita.com/kaizen_nagoya/items/2e41634f6e21a3cf74eb

AUTOSAR(0)Qiita記事一覧, OSEK(75)
https://qiita.com/kaizen_nagoya/items/89c07961b59a8754c869

プログラマが知っていると良い「公序良俗」
https://qiita.com/kaizen_nagoya/items/9fe7c0dfac2fbd77a945

LaTeX(0) 一覧　
https://qiita.com/kaizen_nagoya/items/e3f7dafacab58c499792

自動制御、制御工学一覧（０）
https://qiita.com/kaizen_nagoya/items/7767a4e19a6ae1479e6b

Rust(0) 一覧　
https://qiita.com/kaizen_nagoya/items/5e8bb080ba6ca0281927

100以上いいねをいただいた記事16選
https://qiita.com/kaizen_nagoya/items/f8d958d9084ffbd15d2a

小川清最終講義、最終講義（再）計画, Ethernet(100) 英語(100) 安全(100)
https://qiita.com/kaizen_nagoya/items/e2df642e3951e35e6a53

＜この記事は個人の過去の経験に基づく個人の感想です。現在所属する組織、業務とは関係がありません。＞
This article is an individual impression based on my individual experience. It has nothing to do with the organization or business to which I currently belong.

文書履歴(document history)

ver. 0.01 初稿 20210326
ver. 0.02 書籍表紙、p487-2.f 追記 20210327
ver. 0.03 docker記事、p491.f 追記　@cure_honey さんの指摘により訂正。20210328
ver. 0.04 p491-2.f 追記 20210329
ver. 0.05 p502.f 追記 20210330
ver. 0.06 書式補正 20221211

最後までおよみいただきありがとうございました。

いいね　💚、フォローをお願いします。

Thank you very much for reading to the last sentence.

Please press the like icon 💚　and follow me for your happy life.

このエントリーをはてなブックマークに追加
https://b.hatena.ne.jp/guide/bbutton

4
5
4

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
Sign upLogin
4
5

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?