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

                                                                 30

                    3. ПРОГРАММА МАТЕМАТИЧЕСКОГО
                    МОДЕЛИРОВАНИЯ НА ЯЗЫКЕ БЕЙСИК

10 REM ПРОГРАММА G5,РАЗРАБОТКА ЧЕРНОГО А. А.
20 DIM F(50),H(50),L(50),Y(27),I(50),K(50),M(50),P(50)
30 DIM Q(50),U(50),V(50),О(27), В(27),Z(50), G(20), T(27)
35 DIM K6(50),K7(50),K8(50)J7(50),J8(50),J9(50)
40 PRINT"X=9"
42 J=0: X=0: F3=0: F4=0: H3=0: H4=0: L3=0: L4=0: K3=0: K4=0
44 K5=0: K6=0: K7=0: K8=0: K0=0: X0=0: Y0=0
46 J5=0:J6=0:J7=0:J8=0:J9=0:S=0
50 PRINT"IF X=0 GOTO 6830":PRINT:PRINT"ВВОД X"
60 INPUT X:PRINT"X=";X
61 РRINT"ВВОД ВЕЛИЧИН ФАКТОРОВ И ПОКАЗАТЕЛЕЙ СТЕПЕНИ"
62 FOR J=l TO X:F(J)=0:H(J)=0:L(J)=0:Z(J)=0:K5(J)=0
64 K6(J)=0:K7(J)=0:K8(J)=0:J7(J)=0:J8(J)=0:J(9)=0:NEXT J
110 IF X=0 GOTO 6830
120 PRINT "ВВОД A1 El, B1,J1,O1": INPUT A1,E1,B1,J1,O1
130 PRINT A1;E1;B1;J1;O1:A=A1:B=B1:E=E1:N=J1:R=O1
140 GOSUB 3660: V1=V0:U1=U0:Q1=Q0: PRINT V1;U1;Q1
350 PRINT "ВВОД A2,E2,B2,J2,O2": INPUT A2, E2, B2, J2, O2
360 PRINT A2;E2;B2;J2;O2:A=A2:B=B2:E=E2:N=J2:R=O2
370 GOSUB 3660: V2=V0:U2=U0:Q2=Q0: PRINT V2;U2;Q2
620 F(1)=A1:H(1)=A2:F(2)=B1:H(2)=A2:F(3)=A1:H(3)=B2
630 F(4)=B1:H(4)=B2:F(5)=A1:H(5)=E2:F(6)=B1:H(6)=E2
640 F(7)=E1:H(7)=A2:F(8)=E1:H(8)=B2:F(9)=E1:H(9)=E2
1130 PRINT "IF I0=6 GOTO 40":PRINT "IF I0=7 GOTO 1160"
1140 INPUT I0:IF I0=6 GOTO 40
1150 IF I0=7 GOTO 1160
1160 FOR J=l TO X:PRINT "ВВОД Y(":J:")": INPUT Y(J)
1170 PRINT"Y(":J;")="; Y(J): NEXT J
1180 PRINT "IF I0=1 GOTO 1160":PRINT "IF I0=2 GOTO 1210"
1190 PRINT "ВВОД I0":INPUT I0:IF I0=1 GOTO 1160
1200 IF I0=2 GOTO 1210
1210 GOTO 1340
1340 GOSUB 4150:GOSUB 4170
1410 S=0:O(1)=0:FOR J=l TO X:S=S+Y(J):O(1)=O(1)+1:NEXT J
1420 B(1)=S/O(1):S=0:O(2)=0:FOR J=l TO X:S=S+I(J)×Y(J)
1430 O(2)=O(2)+I(J)^2:NEXT J:B(2)=S/O(2):S=0:O(3)=0
1440 FOR J=l TO X:S=S+K(J) ×Y(J):O(3)=O(3)+K(J)^2:NEXT J
1450 B(3)=S/O(3)
1480 S=0:O(4)=0:FOR J=l TO X:S=S+P(J)xY(J)
1490 O(4)=O(4)+P(J)^2:NEXT J:B(4)=S/O(4):S=0:O(5)=0
1500 FOR J=l TO X:S=S+I(J)×P(J) ×Y(J):O(5)=O(5)+(I(J) ×P(J))^2
1510 NEXT J:B(5)=S/O(5):S=0:O(6)=0:FOR J=l TO X
1520 S=S+Q(J) ×Y(J):O(6)=0(6)+Q(J) ×2:NEXT J:B(6)=S/O(6)