Page 58 - 4729

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

57