Q&A
つまり格子点の数は16個。
だと思うのであれば、以下を先にお読みください。
https://alpha-katekyo.jp/tips/tips100/
縦横長さ３cmの正方形の中に1cm間隔で格子点が配置されているとします。つまり格子点の数は16個。その配置された格子点の中の左上のものの番号を１とします。
この番号1のついた格子点とそのほかの格子点の距離の合計は
do i = 1,16
r(1) = r(1) + sqrt((x(i)-x(1))**2 + (y(i)-y(1))**2)
enddo
とあらわされます。(x(i),y(i)は各格子点のx,y座標の値を格納したもの。例えば格子点1の時はx(1),y(1)とあらわされます。)これを1から16のすべての場所について計算し、r(k)（k番目の格子点とそのほかの格子点との距離の合計）に格納したいときどのようにしたらいいのかお聞きしたいです。もちろん
do k = 1,16
do i = 1,16
r(k) = r(k) + sqrt((x(i)-x(k))**2 + (y(i)-y(k))**2)
enddo
enddo
は試してみましたが、うまくいきません。
上記のプログラムを実行すると結果は256(16*16通り）出てきてしまいます。
回答1件
0
ベストアンサー
これを1から16のすべての場所について計算し、r(k)（k番目の格子点とそのほかの格子点との距離の合計）に格納したい
そのまま書けば良いと思います。
Fortran
1do k = 1,16 2r(k) = 0 3do i = 1,16 4r(k) = r(k) + sqrt((x(i)-x(k))**2 + (y(i)-y(k))**2) 5enddo 6enddo
投稿2022/01/29 15:22
総合スコア84533
otnさんの指摘のように回答のr(k) = r(k) + sqrt((x(i)-x(k))**2 + (y(i)-y(k))**2)の右辺のr(k)の部分を追加して実行しましたが、結果は16*16の256通り出てしまいます。
解決しました。r(i)の出力位置(write(*,*)の位置)を
do k = 1,16
r(k) = 0
do i = 1,16
r(k) = r(k) + sqrt((x(i)-x(k))**2 + (y(i)-y(k))**2)
write(*,*)r(i)
enddo
enddo
としていたことが原因でした。
do k = 1,16
r(k) = 0
do i = 1,16
r(k) = r(k) + sqrt((x(i)-x(k))**2 + (y(i)-y(k))**2)
enddo
enddo
do i = 20
write(*,*)r(i)
enddo
とすると
7.6986150401005802E-004
6.4590329555585017E-004
6.4590329555585017E-004
7.6986150401005802E-004
6.4590329555585006E-004
5.0859106568475470E-004
5.0859106568475470E-004
6.4590329555585038E-004
6.4590329555585017E-004
5.0859106568475481E-004
5.0859106568475470E-004
6.4590329555585027E-004
7.6986150401005824E-004
6.4590329555585027E-004
6.4590329555585038E-004
7.6986150401005835E-004
0.0000000000000000
0.0000000000000000
0.0000000000000000
0.0000000000000000
となり16個の結果が出力できました。
あなたの回答
tips
プレビュー
