Математическое моделирование. Черный А.А. - 215 стр.

3625 PRINT "ВВОД I0": INPUT I0
3630 IF I0 = 10 GOTO 3240
3640 IF I0 = 11 GOTO 4880
3650 IF I0 = 12 GOTO 6830
3653 IF I0 = 14 GOTO 7000
3860 N0 = (A ^ N + B ^ N + C ^ N + D ^ N + E ^ N) / 5
3870 R0 = (A ^ R + B ^ R + C ^ R + D ^ R + E ^ R) / 5

3880 S0 = (A ^ S + B ^ S + C ^ S + D ^ S + E ^ S) / 5
3890 W0 = (A ^ W + B ^ W + C ^ W + D ^ W + E ^ W) / 5
3900 L2 = 2 * N: N3 = (A ^ L2 + B ^ L2 + C ^ L2 + D ^ L2 + E ^ L2) / 5
3910 K2 = 2 * R: R3 = (A ^ K2 + B ^ K2 + C ^ K2 + D ^ K2 + E ^ K2) / 5
3920 M2 = 2 * S: S3 = (A ^ M2 + B ^ M2 + C ^ M2 + D ^ M2 + E ^ M2) / 5
3930 N4 = N + R: N5 = (A ^ N4 + B ^ N4 + C ^ N4 + D ^ N4 + E ^ N4) / 5
3940 N6 = N + S: N7 = (A ^ N6 + B ^ N6 + C ^ N6 + D ^ N6 + E ^ N6) / 5
3950 N8 = N + W: N9 = (A ^ N8 + B ^ N8 + C ^ N8 + D ^ N8 + E ^ N8) / 5
3960 R4 = R + S: R5 = (A ^ R4 + B ^ R4 + C ^ R4 + D ^ R4 + E ^ R4) / 5
3970 R6 = R + W: R7 = (A ^ R6 + B ^ R6 + C ^ R6 + D ^ R6 + E ^ R6) / 5
3980 S4 = S + W: S5 = (A ^ S4 + B ^ S4 + C ^ S4 + D ^ S4 + E ^ S4) / 5
3990 V0 = -N0: U0 = (N0 * R0 - N5) / (N3 - N0 ^ 2): Q0 = -(R0 + U0 * N0)
4000 P0 = (N0 * S0 - N7) / (N3 - N0 ^ 2): Z1 = R0 * S0 - R5 + P0 * (N0 * R0 - N5)
4010 Z2 = U0 * (N0 * S0 - N7) + U0 * P0 * (N0 ^ 2 - N3)
4020 Z3 = R3 - R0 ^ 2 + 2 * U0 * (N5 - N0 * R0)
4030 I0 = (Z1 + Z2) / (Z3 + (N3 - N0 ^ 2) * U0 ^ 2): M0 = I0 * U0 + P0
4040 F0 = -(S0 + I0 * R0 + M0 * N0): Z4 = R0 + U0 * N0
4050 Z5 = Z4 * N0 - N5 - U0 * N3: Z6 = R3 + U0 * N5 - Z4 * R0 - Z5 * U0
4060 Z7 = Z4 * S0 + Z5 * P0 - R5 - U0 * N7: Z0 = (N0 * W0 - N9) / (N3 - N0 ^ 2)
4070 Z8 = Z5 * Z0 + Z4 * W0 - R7 - U0 * N9: Z9 = S3 + I0 * R5 + M0 * N7
4080 T7 = R5 + I0 * R3 + M0 * N5: T8 = N7 + I0 * N5 + M0 * N3
4090 T9 = S0 + I0 * R0 + M0 * N0: G3 = S5 + I0 * R7 + M0 * N9
4100 G4 = T9 * N0 - T8: G5 = Z9 - T9 * S0 - G4 * P0
4110 G6 = T9 * R0 - T7 + G4 * U0: G7 = G4 * Z0 + T9 * W0 - G3
4120 G0 = (Z6 * G7 + Z8 * G6) / (Z6 * G5 - Z7 * G6): H0 = (G0 * Z7 + Z8) / Z6
4130 K0 = G0 * P0 + H0 * U0 + Z0: L0 = -(W0 + G0 * S0 + H0 * R0 + K0 * N0)
4140 PRINT #1, "КОЭФФИЦИЕНТЫ ОРТОГОНАЛИЗАЦИИ": RETURN
4290 FOR J = 1 TO X: I(J) = F(J) ^ J1 + V1
4300 K(J) = F(J) ^ O1 + U1 * F(J) ^ J1 + Q1
4310 L(J) = F(J) ^ P1 + I1 * F(J) ^ O1 + M1 * F(J) ^ J1 + F1
4320 M(J) = F(J) ^ T1 + G1 * F(J) ^ P1 + H1 * F(J) ^ O1 + K1 * F(J) ^ J1 + L1
4330 NEXT J: RETURN
4340 FOR J = 1 TO X: P(J) = H(J) ^ J2 + V2
4350 Q(J) = H(J) ^ O2 + U2 * H(J) ^ J2 + Q2
4360 U(J) = H(J) ^ P2 + I2 * H(J) ^ O2 + M2 * H(J) ^ J2 + F2
4370 V(J) = H(J) ^ T2 + G2 * H(J) ^ P2 + H2 * H(J) ^ O2 + K2 * H(J) ^ J2 + L2
4380 NEXT J: RETURN
4420 FOR J = 1 TO X
4430 Z(J) = B(1) + B(2) * I(J) + B(3) * K(J) + B(4) * L(J) + B(5) * M(J)
4440 NEXT J: RETURN
4690 FOR J = 1 TO X: N3 = B(1) + B(2) * I(J) + B(3) * K(J) + B(4) * P(J)
4700 N4 = B(5) * I(J) * P(J) + B(6) * Q(J) + B(7) * I(J) * Q(J) + B(8) * P(J) * K(J)
4710 N5 = B(9) * K(J) * Q(J) + B(10) * U(J) + B(11) * I(J) * U(J)
4720 N6 = B(12) * K(J) * U(J) + B(13) * V(J) + B(14) * I(J) * V(J) + B(15) * K(J) * V(J)
4730 N7 = B(16) * L(J) + B(17) * P(J) * L(J) + B(18) * Q(J) * L(J) + B(19) * L(J) * U(J)
4740 R3 = B(20) * L(J) * V(J) + B(21) * M(J) + B(22) * P(J) * M(J)
4750 R4 = B(23) * Q(J) * M(J) + B(24) * U(J) * M(J) + B(25) * M(J) * V(J)
4760 Z(J) = N3 + N4 + N5 + N6 + N7 + R3 + R4: NEXT J: RETURN
4880 PRINT #1, "МАТЕМАТИЧЕСКАЯ МОДЕЛЬ"
4902 IF X = 5 GOTO 6100
4907 IF X = 25 GOTO 6480
6100 PRINT #1, "Z(J)="; B(1); "+"; B(2); "*I(J)+"; B(3); "*K(J)+"



                                                  215