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

![41Hb9INzhYL._SL500_.jpg](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/51423/10c2e4c3-04a9-18f2-cf22-aa7f776cf64f.jpeg)
演習 大学院入試問題[数学]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

```ubuntu: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
```
ソースコードは下記。

```fortran: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
```

単純な打ち誤り。

```ubuntu/bash
# 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)だった。

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

```ubuntu/bash
./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にしてみる。

```ubuntu/bash
# 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()に変えて、コンパイル。実行。

```ubuntu: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
```

現状のソースは

```fortran: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

```ubuntu/bash
# 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
```

ソースは

```fortran: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

再度コンパイル

```ubuntu/bash
# 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
# 
```
この時点のソース

```fortran: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

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

現在のソース

```fortran: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

```ubuntu/bash
# 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
```

最初の入力は

```fortran: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

```ubuntu: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
```
文字の誤植を直し、コンパイル、実行。

```ubuntu: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
```

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

```fortran: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

```utunbu/bash
# 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)
```

```fortran: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:

  
```
最初のソースは

```fortran: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
```


```ubuntu: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
```

現状のソース

```fortran: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.

<a href="https://b.hatena.ne.jp/entry/s/qiita.com/kaizen_nagoya/items/67e877bbead1e53b5c24" class="hatena-bookmark-button" data-hatena-bookmark-layout="basic-label-counter" data-hatena-bookmark-lang="ja" title="このエントリーをはてなブックマークに追加"><img src="https://b.st-hatena.com/images/v4/public/entry-button/button-only@2x.png" alt="このエントリーをはてなブックマークに追加" width="20" height="20" style="border: none;" /></a><script type="text/javascript" src="https://b.st-hatena.com/js/bookmark_button.js" charset="utf-8" async="async"></script>
https://b.hatena.ne.jp/guide/bbutton