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

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
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
4880 PRINT #1, "МАТЕМАТИЧЕСКАЯ МОДЕЛЬ"
4902 IF X = 5 GOTO 6100
6100 PRINT #1, "Z(J)="; B(1); "+"; B(2); "*I(J)+"; B(3); "*K(J)+"
6110 PRINT #1, "+"; B(4); "*L(J)+"; B(5); "*M(J),"
6120 IF X = 5 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-й ФАКТОР "
6680 PRINT #1, "M(J)=F(J)^"; T1; "+"; G1; "*F(J)^"; P1; "+"
6690 PRINT #1, "+"; H1; "*F(J)^"; O1; "+"; K1; "*F(J)^"; I1; "+"; L1
6691 PRINT #1, "ОБОЗНАЧЕНИЕ: F(J)- 1-й ФАКТОР "
6700 IF X = 5 GOTO 6790
6790 PRINT "IF I0=18 GOTO 2660-ПЕРЕХОДЫ"
6792 PRINT "IF I0=19 GOTO 3190-ПЕРЕХОДЫ "
6793 PRINT "IF I0=35 GOTO 1160-ВВОД НОВЫХ Y(J)"
6795 PRINT "IF I0=44 GOTO 6830-КОНЕЦ"
6796 PRINT "IF I0=50 GOTO 40-НАЧАЛО"
6797 PRINT "IF I0=51 GOTO 3240-ПРОВЕРКА ТОЧНОСТИ И "
6798 PRINT "              РАСЧЕТЫ ПО МОДЕЛИ"
6799 PRINT "IF I0=52 GOTO 7000-"
6800 PRINT "           ВЫЧИСЛЕНИЯ ПОКАЗАТЕЛЕЙ Z(K5) "
6802 PRINT "              С ИСПОЛЬЗОВАНИЕМ ЦИКЛОВ "
6803 PRINT "                И ПОСТРОЕНИЕ ГРАФИКОВ"
6805 PRINT "ВВОД I0": INPUT I0
6810 IF I0 = 18 GOTO 2660
6820 IF I0 = 19 GOTO 3190
6823 IF I0 = 35 GOTO 1160
6825 IF I0 = 44 GOTO 6830




                                                     90