CALL KHORD (X1,X2,EPS,X,K)

CALL IDRES(X1,X2,EPS,X,K)

10 STOP

END

 

SUBROUTINE KHORD(A,B,EPS,X,K)

* РЕШЕНИЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙ

* МЕТОДОМ ХОРД

PRINT*,'ВВЕДИТЕ ДОПУСКАЕМУЮ ПОГРЕШНОСТЬ EPS'

READ(5,*)EPS

K=0

X1=A

X2=B

1 P=F(X1)

2 Q=F(X2)

X=X1-(X2-X1)*P/(Q-P)

K=K+1

C=F(X)

IF (ABS (C).LE.EPS) GO TO 5

IF (ABS (X2-X1).LE.2*EPS) GO TO 4

IF (C*Q.LT.0.) GO TO 3

X2=X

GO TO 2

3 X1=X

GO TO 1

4 X=(X1+X2)*0.5

5 RETURN

END

 

Результаты вычислений в файле “FT07F001”

 

РЕШЕНИЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙ МЕТОДОМ ХОРД

РЕШАЕМОЕ УРАВНЕНИЕ Y = COS(2*X)-0.25

 

-------------- ИСХОДНЫЕ ДАННЫЕ------------------

НАЧАЛО ОТРЕЗКА XN= -2.500

КОНЕЦ ОТРЕЗКА XK= -2.000

ДОПУСКАЕМАЯ ПОГРЕШНОСТЬ EPS= 0.000100

 

--------------РЕЗУЛЬТАТЫ РАСЧЕТА-----------------

ИСКОМЫЙ КОРЕНЬ X= -2.482525

ЗНАЧЕНИЕ ФУНКЦИИ F(X)= 0.00005

КОЛИЧЕСТВО ИТЕРАЦИЙ K= 2.

======================================