justify">FORMAT(6X,«X»,9X,5H[T1]«,10X,5H[T2]»,10X,4H[T1],11X,4H[T2]/)(3,64)(COSU(I),RM1(I),RM1(I+N),RM2(I),RM2(I+N),I=1,N)
FORMAT (1X, F9.6, 2X, E13.6, 2X, E13.6, 2X, E13.6, 2X, E13.6) (3,4) 33 I=1, N2
R0 (I)=R (I) 34 K=1, NG=ACOS (XG (K))=PI-UK=0.=0.35 I=1, N=0.36 J=1 , M
S2=S2 + COS (J * U (I)) * SIN (J * UK) / J=PUKR - 2. * S2=S11 + R (I) * PUS
S12=S12 + R (I + N) * PUSG (K)=S11 / N
T2G (K)=S12 / N=1 (IPECH.EQ.0) GO TO 37 (3,38)
FORMAT(5X,«XG»,9X,«WG»,13X,«T1G»,12X,«T2G»/)(3,39)(XG(I),WG(I),T1G(I),T2G(I),I=1,NG)
FORMAT (1X, F9.6, 2X, E13.6, 2X, E13.6, 2X, E13.6)
CONTINUE (3,4) (3,4) PT (TXY)! ПРОГРАМА ВИВЕДЕННЯ Т1 У ПЕРШОМУ КВАДРАТІ
DIMENSION TM (300) 41 J=1,30=0.=0.1 * J (J)=TXY (X1, Y1)
CONTINUE (3,42) (TM (I), I=1,30)
FORMAT (10 (1X, E11.4)) (3,43)
FORMAT (/) 45 K=1,3=044 J=1,30=31-J=10 * (K - 1) +1=10 * (K - 1) +1044 I =JMIN, JMAX=0.1 * JY=0.1 * I=IN +1 (IN)=TXY (X1, Y1)
CONTINUE (3,42) (TM (I), I=1,300) (3,43)
CONTINUET1XY (X1, Y1)! РОЗРАХУНОК Т1 В ОДНІЙ ТОЧКЕHR1,HR2,HL,BIP,BIMXG(10),WG(10),T1G(20),T2G(20)/CTCH/NG,XG,WG/CT/AHL,BHL,A2,B2/CTCHT/T1G,T2G/CM/CM1,CM2,CM2Z,CM0,CABRAR1,REZ(AR1,AR),(REZ,RR)PI/3.1415926/=(0.,0.)=(0.,0.)=0.1 I=1, NG=SQRT (Y1 * Y1 + (X1-XG (I)) ** 2)=BHL * RGNVMN (1,0, AR1, REZ)
SL1=(CM1 * B2-3. * CM0) * RR=B2 * RR
AR=AHL*RGNVMN(1,0,AR1,REZ)=(CM1*A2-3.*CM0)*RR=A2*RR=-Y1/R**2+HL**2*Y1/CABR*(SL1-SLG1)=(HL**2*Y1*CM2Z/CABR)*(SL2-SLG2)=S+WG(I)*(T1G(I)*S1+T2G(I)*S2)
CONTINUEXY=- 1. * S/PI/2.
CIMQ
SUBROUTINE SIMQ (A, B, N, KS) A (1), B (1)=0.0=0 =-N65 J=1, N=J +1=JJ + N +1 =0=JJ-J30 I=J, N=IT + I (ABS (BIGA)-ABS (A (IJ))) 20,30,30
BIGA=A (IJ)=I
CONTINUE (ABS (BIGA)-TOL) 35,35,40
KS=1
II=J + N * (J - 2)=IMAX-J50 K=J, N=II + N=II + IT=A (II) (II)=A (I2) ( I2)=SAVE
A (II)=A (II) / BIGA=B (IMAX) (IMAX)=B (J) (J)=SAVE / BIGA (JN) 55,70,55
IQS=N * (J - 1) 65 IX=JY, N=IQS + IX=J-IX60 JX=JY, N=N * (JX - 1) + IX=IXJX + IT
A (IXJX)=A (IXJX) - (A (IXJ) * A (JJX))
B (IX)=B (IX) - (B (J) * A (IXJ))
NY=N - 1=N * N80 J=1, NY=IT-J=NJ=N80 K=1, J (IB)=B (IB)-A (IA) * B (IC)=IA-N
IC=IC - 1
IGNV
SUBROUTINE IGNV (IN, IV, ZI, YK)
INTEGER V
FORMAT (5X, «ARGUMENT <0.0»)
N=IN=IV=ZI=0.1E - 06=0.1E +01=0.0=0.0=0.1E +01=0.0=- 0.1E +01=N + V (NV) 1 , 2,2
N=0=1
IF (Z) 3,4,5
WRITE (3,300) TO 64
YK=0.0TO 64
IF (Z - 5.5) 6,7,7
IF (Z-NV - 2) 6,6,8
IF (5.5-NV - 2) 9,9,10