Call idres(x1,x2,eps,x,k)

10 STOP

END

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

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

* МЕТОДОМ ПОЛОВИННОГО ДЕЛЕНИЯ

K=0

X1=A

X2=B

1 X=(X1+X2)*0.5

K=K+1

C=F(X)

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

IF ((X2-X1).LE.EPS) GO TO 3

IF (C*F(X2).LT.0.) GO TO 2

X2=X

GO TO 1

2 X1=X

GO TO 1

3 X=(X1+X2)*0.5

4 RETURN

END

SUBROUTINE IDRES (XN,XK,EPS,X,K)

* РАСПЕЧАТКА ИСХОДНЫХ ДАННЫХ И РЕЗУЛЬТАТОВ РАСЧЕТА

C=F(X)

DO 20 I=6,7

WRITE(I,1)

1 FORMAT(//2X,16(‘-‘), ’ ИСХОДНЫЕ ДАННЫЕ’, 16(‘-‘))

WRITE(I,2)XN,XK,EPS

2 FORMAT (/2X,’НАЧАЛО ОТРЕЗКА XN=’,F6.3/2X,

* ’КОНЕЦ ОТРЕЗКА XK=’,F6.3/2X,

* ‘ДОПУСКАЕМАЯ ПОГРЕШНОСТЬ EPS=’,F8.6)

WRITE(I,3)

3 FORMAT(//2X,14(‘-‘),’РЕЗУЛЬТАТЫ РАСЧЕТА’,14(‘-‘))

WRITE(I,4)X,C,K

4 FORMAT(/2X,’ИСКОМЫЙ КОРЕНЬ ‘,6X,’X=’,F10.6/

*2X,’ЗНАЧЕНИЕ ФУНКЦИИ’,3X,’F(X)=’,F10.6/

*2X,’КОЛИЧЕСТВО ИТЕРАЦИЙ K=’,I3)

WRITE(I,5)

5 FORMAT(40('='))

20 CONTINUE

END

Примечание. Автоматически созданный файл результатов “ FT07F001”

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

МЕТОДОМ ПОЛОВИННОГО ДЕЛЕНИЯ

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

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

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

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

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

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

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

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

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

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

Примечание. Текст подпрограммы OTKOR приведен в приложении 1.

ПРИМЕР РЕШЕНИЯ УРАВНЕНИЯ МЕТОДОМ ХОРД

ПРИЛОЖЕНИЕ 3

PROGRAM RESH

* Главная программа

DO 5 N=6,7

WRITE(N,1)

1 FORMAT(//15X,'РЕШЕНИЕ НЕЛИНЕЙНЫХ УРАВНЕНИЙ'/

*15X,МЕТОДОМ ХОРД')

WRITE(N,2)

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

5 CONTINUE

CALL OTKOR(N,X1,X2)

IF(N.GT.100) GO TO 10

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.

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