Page 12 - 4824
P. 12
30 " ЧАСТКОВІ ПОХІДНІ G(l), G(2),.... G(N)
ОБЧИСЛЮЮТЬСЯ,
35 • ПОЧИНАЮЧИ 3 3000 РЯДКА
40 " ПОШУК ДОВЖИНИ КРОКУ ПРОВОДИТЬСЯ
МЕТОДОМ КУБІЧНОЇ ІНТЕРПОЛЯЦІЇ
50 INPUT "ВВЕДІТЬ ЧИСЛО ЗМІННИХ N = ";N —60 PRINT
"ЧИСЛО ЗМІННИХ N = ";N
70 DIM X(N), Y(N), P(N), Q(N), D(N), G(N)
80 PRINT "ВВЕДІТЬ ПОЧАТКОВУ ТОЧКУ"
90 FOR i = 1 TO N: PRINT "U";I;" = "; : INPUT X(I):NEXT I
350 SV=1 : TV=0 360 " ПОМІЖНИИ ВВІД 380 PRINT
"ПРОМІЖНІ ЗНАЧЕННЯ"
550 FOR i = 1 TO N:P(I)=X(I) : PRINT "U",i = "X(I) : NEXT I
560 GOSUB 2000: FP = Z: PRINT "R(U)="Z
570 GOSUB 3000 : G1 = GO : GK=GO
575 " ЗА ПЕРШИЙ НАПРЯМОК ВЗЯТИ НАПРЯМОК
НАЙШВИДШОГО СПУСКУ
580 FOR I=1 ТО N:D(I)=—G(I) :NEXT I
585 "К—ЛІЧИЛЬНИК ІТЕРАЦІЙ
590 К=1
600 GP = 0
610 FOR I=1 ТО N:GP = GP + G(I)*D(I) : NEXT I
620 IF GP<=0 THEN 680
625 "ВИЗНАЧИТИ ПОЧАТКОВИЙ КРОК І, ЯКЩО
НЕОБХІДНО,
626 "ЗМHИТИ НАПРЯМОК СПУСКУ НА ПРОТИЛЕЖНИЙ
630 QX = ABS(2*FP/GP) : IF QX>1 THEN QX=1
640 FOR I = 1 TO N
650 X(I)=P(I)—QX*D(I) : P(I)=X(I) : NEXT I
660 GOSUB 2000: FP = Z: PRINT "НЕСТАБІЛЬНІСТЬ"
670 GPSUB 3000 : G1 = GO : GOTO 600
680 QX = ABS(2*FP/GP):IF QX>1 THEN QX=1
690 HH==QX
700 " ЗНАХОДИМО ЧЕРГОВУ ТОЧКУ
710 BB = HH
720 FOR 1=1 ТО N
730 Q(I)=P(I)+BB*D(I) :X(I)=Q(I) : NEXT I