Главная -> Книги (0) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68) (69) (70) (71) (72) (73) (74) (75) (76) (77) (78) (79) (80) (81) (82) ( 83 ) (84) (85) (86) (87) (88) (89) (90) (91) (92) (93) (94) (95) (96) (97) (98) (99) (100) (101) (102) (103) (104) (105) (83)
С ***************************************************************** 0002 C0MM0N/DATA1/DAT1(16) 0003 DIMENSION TEX1(16) 0004 ITOT=IPA+IC+IR 0005 DO 200,1 = 1,ITOT 0006 TYPE 10,TEX1(I),DAT1(I) 0007 200 CONTINUE 0008 RETURN 0009 10 FORMAT(X,A4,T15,=,1PE10.3) 0010 END 0001 SUBROUTINE LPLQ С ***************************************************************** 0002 COMMON/STO/STOP 0003 C0MM0N/DATA1/XFR,XQP,XK,XGSP,XP,XR11,XR12,XC2, XR3 XC4 0004 DIMENSION TEXl (16),IN (6) 0005 LOGICAL*! STOP 0006 DATA TEX1/FREQ,QP,K,GSP/P,RI r,R12,C2,R3,C4, <""""" "/IN/1,2,3,0,8,10/ 0007 DATA IFR,IV,IC,IRIN.IR0UT,IPA,ICH/1 3,2,0,3,5,0/ 0008 CALL INPUT(IFR,IV,IC,IRIN,ICH,TEX1,IN) 0009 100 XKIN = XK 0010 A=XC2/(XC4*2.*XQP*XQP) ООН IF(A.GE.2) GOTO 150 0013 TYPE 5 0014 5 FORMAT (/,X,C2>=C4*4*QP**2) 0015 CALL TEST(IFR,IV,IC,1RIN,TEX1,IN) 0016 GOTO 100 0017 150 XCHI = XC4/XC2 0018 XQHI=XQP*XQP 0019 XR= (l./(XCHI*2.*XQHI)-l.)+SQRT((l./(2*yQHI*XCM)- 1)**2-1.) 0020 160 XR1 = 1/(SQRT(XP*XC2*XC4)*(6.2831853*XFR)) 0021 XR3=XP*XR1 0022 XGSP=XQP*SQRT(1./(XCHI*XP)) 0023 IF (XKIN EQ.O) GOTO 200 0025 IF ((1-XKIN).LE.O) GOTO 200 0027 XK=XKIN 0028 XR11=XR1/XK 0029 XR12=XR1/(1-XK) 0030 GOTO 300 0031 200 XR11=XR1 0032 XR12=1.E38 0033 XK=1. 0034 300 CALL OUTPUT(IPA,IC,IROUT,TEXI) 0035 XK=XKIN 0036 CALL TEST(IFR,IV,IC,IRIN,TEX1.IN) 0037 IF(.NOT.STOP) GOTO 100 0039 RETURN 0040 END 0001 SUBROUTINE APLQS С ***************************************************************** COMMON/STO/STOP COMMON/DATA 1 /XFR,XQP,XK,XRO.XC 1 ,XC2,XC3,XR4,XR5, XR6 DIMENSION TEX1(16),IN(7) LOGICAL*! STOP DATA TEXl/FREQ,QPVK,R0,Cr.C2VC3,R4,R5,R6, Г/ V V V V 7IN/1,2,0,5,6,7,4/ DATA IFR,IV,IC,IRIN,IROUT,IPA,ICH/1,2,3,I,4,3,0/ TYPE 10 FORMAT (X,K=1,T20,(TYPE 1), ACCEPT 11,XK FORMAT (FIO.O) CALL INPUT(IFR,IV,IC,IRIN,ICH,TEX1,IN) XH= 1.-((XCl +XC2) * (XC2+XC3) *4.*XQP*XQP/(XCl*XC3)) IF(XH.GE.O.) GOTO 150 TYPE 5 F0RMAT(/.X,C1=C3> = C2*2*QP/(1-2*QP)) CALL TEST(IFR,IV,IC,IRIN,TEX1,IN) GOTO 100 XW=XFR*6 2831853 IF(XK.EQ.l) XH1 = 1.+SQRT(XH) IF(XK.EQ.-l)XHl = l.-SQRT(XH) XR4=XH1/ (2.*XW*XQP* (XC2+XC3)) XR5 = (XC1+XC2) / (XR4* XW* XW* XC1 *XC2*XC31
0027 XR6 = l./(XW*XW*XCrxC2*(XR4+XR5)) 0028 CALL 0UTPUT(IPA,IC,IR0UT,TEX1) 0029 CALL TEST(IFR,IVJC,IRIN,TEXLIN) 0030 IF(.NOT.STOP) GOTO 100 0032 RETURN 0033 END 0001 SUBROUTINE BPLQR С ***************************************************************** 0002 COMMON/STO/STOP 0003 COMMON/DATA 1/XFR,XQP,XK,XGSP,XP,XR11,XR12,XC2, XC3,XR4 0004 DIMENSION TEXl (16),IN(6) 0005 LOGICAL*! STOP 0006 DATA TEX1/FREQ,QP,K/GSPVP,R11,R12,C2,C3, lR4V, V V , 7IN/1,2,3,0,8,9/ 0007 DATA IFR,IV,IC,IRIN,IROUT,IPA,ICH/1,3,2,0,3,5,0/ 0008 CALL INPUT/IFR,IV,IC,IRIN,ICH,TEX1,IN/ 0009 150 XKIN=XK 0010 XCHI=XC3/XC2 ООП XQHI = XQP*XQP 0012 XP=XQHI*/2.+XCHH-l./XCHI/ 0013 160 XR1 = 1/(SQRT(XP*XC2*XC3)*(62831853*XFR) 0014 XR4=XP*XR1 0015 XGSP = XQHI*(1.+XCHI) 0016 XKO=XGSP 0017 IF (XKIN.EQ.OOR (XKO-XKIN).LE0) GOTO 200 0019 XK=XKIN 0020 XR11 = FR1*XK/XK 0021 XR12=XK*XR11/(XK0-XK) 0022 GOTO 300 0023 200 XR11 = XR1 0024 XR12==1.E38 0025 XK=XKO 0026 300 CALL OUTPUT (IPA,IC,IR0UT,TEX1) 0027 XK=XKIN 0028 CALL TEST(IFR,IV,IC,IRIN,TEXI,IN) 0029 IF(NOT.STOP) GOTO 150 0031 RETURN 0032 END 0001 SUBROUTINE BPLQC С ***************************************************************** 0002 COMMON/STO/STOP 0003 C0MM0N/DATA1/XFR,XQP,XK,XGSP,XP,XC11,XC12,XR2, XR3 XC4 0004 DIMENSION TEXl (16),IN(6) 0005 LOGICAL*! STOP 0006 DATA TEX1/FREQ.QP,KVGSPVPVC11.C12VR2/R3, lC4V,,,,, 7IN/1,2,0,6,7,10/ 0007 DATA IFR,IV,IC,IRIN,IROUT,IPA,ICH/1,2,3,0,2,5,0/ 0008 CALL INPUT(IFR,IV,IC,IRIN,ICH,TEX1,IN) 0009 150 XC!=XC!!+XC!2 0010 A=XC1/(XC4*2.*XQP*XQP) ООП IF(A.GE.2.) GOTO 170 0013 TYPE 5 0014 CALL TEST(IFR,IV.IC,IRIN,TEX1,IN) (0) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68) (69) (70) (71) (72) (73) (74) (75) (76) (77) (78) (79) (80) (81) (82) ( 83 ) (84) (85) (86) (87) (88) (89) (90) (91) (92) (93) (94) (95) (96) (97) (98) (99) (100) (101) (102) (103) (104) (105) |
|