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

3612 PRINT "IF I0=11 GOTO 4880 - МАТЕМАТИЧЕСКАЯ МОДЕЛЬ"
3615 PRINT "IF I0=14 GOTO 7000-ВЫЧИСЛЕНИЯ ПОКАЗАТЕЛЕЙ Z(K5)"
3616 PRINT "             С ИСПОЛЬЗОВАНИЕМ ЦИКЛОВ И "
3617 PRINT "              ПОСТРОЕНИЕ ГРАФИКОВ"
3620 PRINT "IF I0=12 GOTO 6830-КОНЕЦ"
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
3660 N0 = (A ^ N + B ^ N + E ^ N) / 3: R0 = (A ^ R + B ^ R + E ^ R) / 3
3670 L2 = 2 * N: N3 = (A ^ L2 + B ^ L2 + E ^ L2) / 3: N4 = N + R
3680 N5 = (A ^ N4 + B ^ N4 + E ^ N4) / 3: V0 = -N0
3690 U0 = (N0 * R0 - N5) / (N3 - N0 ^ 2): Q0 = -(R0 + U0 * N0)
3700 PRINT #1, "КОЭФФИЦИЕНТЫ ОРТОГОНАЛИЗАЦИИ": RETURN
3710 N0 = (A ^ N + B ^ N + C ^ N + D ^ N) / 4
3720 R0 = (A ^ R + B ^ R + C ^ R + D ^ R) / 4
3730 S0 = (A ^ S + B ^ S + C ^ S + D ^ S) / 4: L2 = 2 * N
3740 N3 = (A ^ L2 + B ^ L2 + C ^ L2 + D ^ L2) / 4: K2 = 2 * R
3750 R3 = (A ^ K2 + B ^ K2 + C ^ K2 + D ^ K2) / 4: N4 = N + R
3760 N5 = (A ^ N4 + B ^ N4 + C ^ N4 + D ^ N4) / 4: N6 = N + S
3770 N7 = (A ^ N6 + B ^ N6 + C ^ N6 + D ^ N6) / 4: R4 = R + S
3780 R5 = (A ^ R4 + B ^ R4 + C ^ R4 + D ^ R4) / 4: V0 = -N0
3790 U0 = (N0 * R0 - N5) / (N3 - N0 ^ 2): Q0 = -(R0 + U0 * N0)
3800 P0 = (N0 * S0 - N7) / (N3 - N0 ^ 2): Z1 = R0 * S0 - R5 + P0 * (N0 * R0 - N5)
3810 Z2 = U0 * (N0 * S0 - N7) + U0 * P0 * (N0 ^ 2 - N3)
3820 Z3 = R3 - R0 ^ 2 + 2 * U0 * (N5 - N0 * R0)
3830 I0 = (Z1 + Z2) / (Z3 + (N3 - N0 ^ 2) * U0 ^ 2): M0 = I0 * U0 + P0
3840 F0 = -(S0 + I0 * R0 + M0 * N0)
3850 PRINT #1, "КОЭФФИЦИЕНТЫ ОРТОГОНАЛИЗАЦИИ": RETURN
4150 FOR J = 1 TO X: I(J) = F(J) ^ J1 + V1
4160 K(J) = F(J) ^ O1 + U1 * F(J) ^ J1 + Q1: NEXT J: RETURN
4250 FOR J = 1 TO X: P(J) = H(J) ^ J2 + V2
4260 Q(J) = H(J) ^ O2 + U2 * H(J) ^ J2 + Q2
4270 U(J) = H(J) ^ P2 + I2 * H(J) ^ O2 + M2 * H(J) ^ J2 + F2: NEXT J
4280 RETURN
4490 FOR J = 1 TO X: N3 = B(1) + B(2) * I(J) + B(3) * K(J) + B(4) * P(J)
4500 N4 = B(5) * I(J) * P(J) + B(6) * Q(J) + B(7) * I(J) * Q(J) + B(8) * P(J) * K(J)
4510 N5 = B(9) * K(J) * Q(J) + B(10) * U(J) + B(11) * I(J) * U(J)
4520 Z(J) = N3 + N4 + N5 + B(12) * K(J) * U(J): NEXT J: RETURN
4880 PRINT #1, "МАТЕМАТИЧЕСКАЯ МОДЕЛЬ"
4903 IF X = 12 GOTO 6130
6130 PRINT #1, "Z(J)="; B(1); "+"; B(2); "*I(J)+"; B(3); "*K(J)+"
6140 PRINT #1, "+"; B(4); "*P(J)+"; B(5); "*I(J)*P(J)+"
6150 PRINT #1, "+"; B(6); "*Q(J)+"; B(7); "*I(J)*Q(J)+"
6160 PRINT #1, "+"; B(8); "*P(J)*K(J)+"; B(9); "*K(J)*Q(J)+"
6170 PRINT #1, "+"; B(10); "*U(J)+"; B(11); "*I(J)*U(J)+"
6180 PRINT #1, "+"; B(12); "*K(J)*U(J),"
6190 IF X = 12 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-й ФАКТОР "
6630 IF X = 12 GOTO 6710
6710 PRINT #1, "P(J)=H(J)^"; J2; "+"; V2; ";"
6720 PRINT #1, "Q(J)=H(J)^"; O2; "+"; U2; "*H(J)^"; J2; "+"; Q2; ";"
6730 PRINT #1, "U(J)=H(J)^"; P2; "+"; I2; "*H(J)^"; O2; "+"




                                                      199