ΠΡΠ·ΠΎΠ²: CALL lndKat_plosk_q (R, Z, ndimrz, nwtk, NACHW, KONW, shr, shz, rindkat).
! SUBROUTINE lndKat_plosk_q (R, Z, ndimrz, nwtk, NACHW, KONW, shr, shz, rindkat).
!ΠΠΏ ΡΠ°ΡΡΡΠΈΡΡΠ²Π°Π΅Ρ ΠΈΠ½Π΄ΡΠΊΡΠΈΠ²Π½ΠΎΡΡΡ ΠΏΠ»ΠΎΡΠΊΠΎΠΉ ΠΊΠ°ΡΡΡΠΊΠΈ Ρ ΠΠΠΠΠΠΠΠΠ«ΠΠ ΡΠ΅ΡΠ΅Π½ΠΈΡΠΌΠΈ Π²ΠΈΡΠΊΠΎΠ²! ΠΈ Ρ ΠΠΠΠΠΠΠΠΠΠ ΡΠ°Π·Π±ΠΈΠ²ΠΊΠΎΠΉ ΡΠ΅ΡΠ΅Π½ΠΈΠΉ ΠΏΠΎ ΡΠΎΡΠΌΡΠ»Π΅.
! ΠΏ ΠΏ-1 ΠΏ.
! Lk=SUM Ljj + 2*SUM SUM Ljm, Π³Π΄Π΅.
! j =1 j =1 m=j + 1.
!Ljj — ΡΠΎΠ±ΡΡΠ²Π΅Π½Π½Π°Ρ ΠΈΠ½Π΄ΡΠΊΡΠΈΠ²Π½ΠΎΡΡΡ ΠΌΠ°ΡΡΠΈΠ²Π½ΠΎΠ³ΠΎ Π²ΠΈΡΠΊΠ°;
!Ljm — Π²Π·Π°ΠΈΠΌΠ½Π°Ρ ΠΈΠ½Π΄ΡΠΊΡΠΈΠ²Π½ΠΎΡΡΡ j-ro ΠΈ m-Π³ΠΎ ΠΌΠ°ΡΡΠΈΠ²Π½ΡΡ
Π²ΠΈΡΠΊΠΎΠ² ΠΊΠ°ΡΡΡΠΊΠΈ;
!ΠΏΠ²Π΅ΡΡ
Π½ΠΈΠΉ ΠΏΡΠ΅Π΄Π΅Π» ΡΡΠΌΠΌΡ — ΠΊΠΎΠ»ΠΈΡΠ΅ΡΡΠ²ΠΎ Π²ΠΈΡΠΊΠΎΠ² Π² ΠΊΠ°ΡΡΡΠΊΠ΅;
!j, mΠΈΠ½Π΄Π΅ΠΊΡΡ ΡΡΠΌΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΡ Π² ΡΡΠΌΠΌΠ°Ρ
.
!ΠΠΈΡΠΊΠΈ ΠΌΠΎΠ³ΡΡ ΠΈΠΌΠ΅ΡΡ ΠΊΠ°ΠΊ ΠΎΠ΄ΠΈΠ½Π°ΠΊΠΎΠ²ΡΠ΅, ΡΠ°ΠΊ ΠΈ Π½Π΅ΠΎΠ΄ΠΈΠ½Π°ΠΊΠΎΠ²ΡΠ΅ ΠΏΠΎΠΏΠ΅ΡΠ΅ΡΠ½ΡΠ΅ ΡΠ΅ΡΠ΅Π½ΠΈΡ.
1Π’Π ΠΠΠ£ΠΠΠ«Π ΠΠ: lndKat_plosk_q, Vzi_WW_hor_qq, REAVZI {ΠΠ‘ΠΠΠΠ¬ΠΠΠΠΠ Π£Π‘ΠΠΠ ΠΠΠΠ«Π ΠΠΠΠΠ ΠΠ’Π! ΠΡΠ°Π²ΠΈΠ»ΡΠ½ΠΎΠ΅ ΡΠ°ΡΠΏΠΎΠ»ΠΎΠΆΠ΅Π½ΠΈΠ΅ Π²ΠΈΡΠΊΠΎΠ² ΠΏΠ»ΠΎΡΠΊΠΎΠΉ ΠΊΠ°ΡΡΡΠΊΠΈ.
! |Z.
! |__.
! I I I I I I I II.
! I1_1 Π« Π« 11
! |_R.
!
! ΠΠ₯ΠΠ:
!R, Z — ΠΎΠ΄Π½ΠΎΠΌΠ΅ΡΠ½ΡΠ΅ ΠΌΠ°ΡΡΠΈΠ²Ρ, ΡΠΎΠ΄Π΅ΡΠΆΠ°ΡΠΈΠ΅ ΡΠ°Π΄ΠΈΡΡΡ ΠΈ z-ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ Π½ΠΈΡΠ΅Π²ΠΈΠ΄Π½ΡΡ
.
! ΠΊΠΎΠ½ΡΡΡΠΎΠ² Π² ΡΠ΅Π½ΡΡΠ°Ρ
ΡΠ»Π΅ΠΌΠ΅Π½ΡΠ°ΡΠ½ΡΡ
ΠΏΠ»ΠΎΡΠ°Π΄ΠΎΠΊ.
Indimrz — ΡΠ°Π·ΠΌΠ΅ΡΠ½ΠΎΡΡΡ ΠΌΠ°ΡΡΠΈΠ²ΠΎΠ² R ΠΈ Z Inwtk — ΠΊΠΎΠ»ΠΈΡΠ΅ΡΡΠ²ΠΎ Π²ΠΈΡΠΊΠΎΠ² Π² ΠΊΠ°ΡΡΡΠΊΠ΅.
!NACHW, KONW — ΠΎΠ΄Π½ΠΎΠΌΠ΅ΡΠ½ΡΠ΅ ΠΌΠ°ΡΡΠΈΠ²Ρ ΡΠ°Π·ΠΌΠ΅ΡΠ½ΠΎΡΡΠΈ nwtk, ΡΠΎΠ΄Π΅ΡΠΆΠ°ΡΠΈΠ΅ Π½ΠΎΠΌΠ΅Ρ ΠΏΠ΅ΡΠ²ΠΎΠΉ ΠΈ.
! ΠΊΠΎΠ½Π΅ΡΠ½ΠΎΠΉ ΡΠ»Π΅ΠΌΠ΅Π½ΡΠ°ΡΠ½ΡΡ
ΠΏΠ»ΠΎΡΠ°Π΄ΠΎΠΊ Π² ΠΎΡΠ΅ΡΠ΅Π΄Π½ΠΎΠΌ Π²ΠΈΡΠΊΠ΅.
!nrr, nrz — ΡΠΈΡΠ»Π° ΡΠ°Π·Π±ΠΈΠ΅Π½ΠΈΠΉ ΡΠ΅ΡΠ΅Π½ΠΈΡ Π²ΠΈΡΠΊΠ° Π²Π΄ΠΎΠ»Ρ ΡΠ°Π΄ΠΈΡΡΠ° ΠΈ Π²Π΄ΠΎΠ»Ρ z-ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ.
!shr, shz — ΡΠ°Π³ΠΈ ΡΠ°Π·Π±ΠΈΠ΅Π½ΠΈΡ Π²ΠΈΡΠΊΠΎΠ² ΠΏΠΎ ΡΠ°Π΄ΠΈΡΡΡ ΠΈ Π²Π΄ΠΎΠ»Ρ z-ΠΊΠΎΠΎΡΠ΄ΠΈΠ½Π°ΡΡ.
! ΠΠ«Π₯ΠΠ:
Irindkat — ΡΠ°ΡΡΡΠΈΡΠ°Π½Π½Π°Ρ ΠΈΠ½Π΄ΡΠΊΡΠΈΠ²Π½ΠΎΡΡΡ ΠΊΠ°ΡΡΡΠΊΠΈ Π² ΠΠ½.
SUBROUTINE lndKat_plosk_q (R, Z, ndimrz, nwtk, NACHW, KONW, nrr, nrz, shr, shz, rindkat) integer: ndimrz, nwtk, NACHW (nwtk), KONW (nwtk), j, m, nachl, konl, nrr, nrz real (8): R (ndimrz), Z (ndimrz), shr, shz, rj, zj, & rindkat, soblnd, s, vziww.
rindkat=OdO; s=0d0 DO j=l, nwtk.
rj=R (j); zj=Z (j); nachl=NACHW (j); konl=KONW (j).
!Π Π°ΡΡΠ΅Ρ ΡΠΎΠ±ΡΡΠ²Π΅Π½Π½ΠΎΠΉ ΠΈΠ½Π΄ΡΠΊΡΠΈΠ²Π½ΠΎΡΡΠΈ Π²ΠΈΡΠΊΠΎΠ² ΡΠ»ΠΎΡ:
CALL lndMWtk_qq (R/Z/ndimrz/nrr/nrz/nachl/konl/shr, shz/ soblnd) rindkat=rindkat+Soblnd Π‘ΡΠΌΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ ΡΠΎΠ±ΡΡΠ². ΠΈΠ½Π΄ΡΠΊΡΠΈΠ²Π½ΠΎΡΡΠ΅ΠΉ Π²ΠΈΡΠΊΠΎΠ² if (j==nwtk)exit do m=j+l, nwtk.
CALL Vzi_WW_hor_qq (R, Z, ndimrz, nachl, konl, NACHW (m), KONW (m), nrr, nrz, vziww) s=s+vziww ! Π‘ΡΠΌΠΌΠΈΡΠΎΠ²Π°Π½ΠΈΠ΅ Π²ΠΎ Π²Π½ΡΡΡ. ΠΈ Π²Π½Π΅ΡΠ½. ΡΠΈΠΊΠ»Π°Ρ
.
enddo enddo.
rindkat=rindkat+2d0*s.
END.