円計算式x² + y² = 1
をx
とy
から-1
から1
に変更し、各ステップで0.1
を増やす円形のグリッドを作成する必要があります。Fortran 90で循環グリッドを作成するにはどうすればよいですか?
私は多くの異なることを試みましたが、失敗しました。誰かが私を助けることができたら嬉しいです。
私の試みの一つであった:attempt one
PROGRAM cir
IMPLICIT real*8(a-h, o-z)
real :: x(20,20), py(20,20), ny(20,20)
delx = -1.1
dely = -1.1
DO i=1,20
delx = delx + 0.1
DO j=1.20
x(i,j) = delx
ny(i,j)= -1.0 * SQRT (1.0 - (x(i,j)**2))
py(i,j)= SQRT (1.0 - (x(i,j)**2) )
PRINT*, x(i,j), ny(i,j), py(i,j)
END DO
END DO
STOP
END PROGRAM cir
PROGRAM circle
IMPLICIT real*8(a-h,o-z)
real*8 :: x(20,20), y(20,20), z(20,20)
delx = -1.0
dely = -1.0
DO i = 1, 20
delx = delx + 0.1
DO j = 1 , 20
dely = dely + 0.1
x(i,j) = delx
y(i,j) = dely
z(i,j) = x(i,j)**2 + y(i,j)**2 -1.0
END DO
delx = -1.0
dely = -1.0
PRINT*, x(i,j), y(i,j)
END DO
STOP
END PROGRAM circle
円形のグリッドはかなり曖昧です。あなたが望む結果のタイプは何ですか? – innoSPG
最終結果として何をしたいですか?ポイントがライン上にある場合はブール値の配列ですか?あなたはいくつかのポイントだけを持っているので、彼らはほとんど決してライン上にあることはありません。ポイントが内部にある場合はブール値の配列ですか?あなたはx²+y²<1と書かれていたはずです。あなたは何を達成したいですか? –
はあなたが試したことを示すために有益なことがあります。 – agentp