Математическое моделирование при планировании экспериментов на четырех уровнях факторов. Черный А.А. - 23 стр.

3432 PRINT #1, " ФАКТОРЫ F("; S; ")="; F(S); "H("; S; ")="; H(S)
3470 IF X = 16 GOTO 3530
3530 GOSUB 4210: GOSUB 4250: GOSUB 4580: GOTO 3552
3552 PRINT #1, "Z("; S; ")="; Z(S): NEXT S: GOTO 3610
3560 FOR S = 1 TO X: F(S) = 0: H(S) = 0: L(S) = 0: Z(S) = 0
3570 PRINT "ВВОД F("; S; "),H("; S; "),L("; S; ")"
3572 INPUT F(S), H(S), L(S): PRINT #1, "ФАКТОР F("; S; ")="; F(S)
3574 PRINT #1, "ФАКТОРЫ H("; S; ")="; H(S); "L("; S; ")="; L(S)
3590 PRINT #1, "Z("; S; ")="; Z(S): NEXT S
3610 PRINT "IF I0=10 GOTO 3240-ПРОВЕРКА ТОЧНОСТИ И ";
3611 PRINT "РАСЧЕТЫ ПО МОДЕЛИ"
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
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
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
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
4400 FOR J = 1 TO X: Z(J) = B(1) + B(2) * I(J) + B(3) * K(J) + B(4) * L(J)
4410 NEXT J: RETURN
4580 FOR J = 1 TO X: N3 = B(1) + B(2) * I(J) + B(3) * K(J) + B(4) * P(J)
4590 N4 = B(5) * I(J) * P(J) + B(6) * Q(J) + B(7) * I(J) * Q(J) + B(8) * P(J) * K(J)




                                              23