Планы проведения экспериментов и компьютерные программы эффективного математического моделирования на языках Бейсик и Турбо Паскаль. Черный А.А. - 208 стр.

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
4210 FOR J = 1 TO X: I(J) = F(J) ^ J1 + V1
4220 K(J) = F(J) ^ O1 + U1 * F(J) ^ J1 + Q1
4230 L(J) = F(J) ^ P1 + I1 * F(J) ^ O1 + M1 * F(J) ^ J1 + F1: NEXT J
4240 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
4630 FOR J = 1 TO X: N3 = B(1) + B(2) * I(J) + B(3) * K(J) + B(4) * P(J)
4640 N4 = B(5) * I(J) * P(J) + B(6) * Q(J) + B(7) * I(J) * Q(J) + B(8) * P(J) * K(J)
4650 N5 = B(9) * K(J) * Q(J) + B(10) * U(J) + B(11) * I(J) * U(J)
4660 N6 = B(12) * K(J) * U(J) + B(13) * V(J) + B(14) * I(J) * V(J) + B(15) * K(J) * V(J)
4670 N7 = B(16) * L(J) + B(17) * P(J) * L(J) + B(18) * Q(J) * L(J) + B(19) * L(J) * U(J)
4680 Z(J) = N3 + N4 + N5 + N6 + N7 + B(20) * L(J) * V(J): NEXT J: RETURN
4880 PRINT #1, "МАТЕМАТИЧЕСКАЯ МОДЕЛЬ"
4906 IF X = 20 GOTO 6370
6370 PRINT #1, "Z(J)="; B(1); "+"; B(2); "*I(J)+"; B(3); "*K(J)+"
6380 PRINT #1, "+"; B(4); "*P(J)+"; B(5); "*I(J)*P(J)+"
6390 PRINT #1, "+"; B(6); "*Q(J)+"; B(7); "*I(J)+Q(J)+"
6400 PRINT #1, "+"; B(8); "*P(J)*K(J)+"; B(9); "*K(J)*Q(J)+"
6410 PRINT #1, "+"; B(10); "*U(J)+"; B(11); "*I(J)*U(J)+"
6420 PRINT #1, "+"; B(12); "*K(J)*U(J)+"; B(13); "*V(J)+"
6430 PRINT #1, "+"; B(14); "*I(J)*V(J)+"; B(15); "*I(J)*V(J)+"
6440 PRINT #1, "+"; B(16); "*L(J)+"; B(17); "*P(J)*L(J)+"
6450 PRINT #1, "+"; B(18); "*Q(J)*L(J)+"; B(19); "*L(J)*U(J)+"
6460 PRINT #1, "+"; B(20); "*L(J)*V(J),"
6470 IF X = 20 GOTO 6600
6600 PRINT #1, "ГДЕ"
6610 PRINT #1, "I(J)=F(J)^"; J1; "+"; V1; ";"
6620 PRINT #1, "K(J)=F(J)^"; O1; "+"; U1; "*F(J)^"; J1; "+"; Q1
6621 PRINT #1, "ОБОЗНАЧЕНИЕ: F(J)- 1-й ФАКТОР "
6650 PRINT #1, "L(J)=F(J)^"; P1; "+"; I1; "*F(J)^"; O1; "+"
6660 PRINT #1, "+"; M1; "F(J)^"; J1; "+"; F1
6661 PRINT #1, "ОБОЗНАЧЕНИЕ: F(J)- 1-й ФАКТОР "
6675 IF X = 20 GOTO 6710
6710 PRINT #1, "P(J)=H(J)^"; J2; "+"; V2; ";"
6720 PRINT #1, "Q(J)=H(J)^"; O2; "+"; U2; "*H(J)^"; J2; "+"; Q2; ";"




                                                     208