![]() |
Главная Случайная страница Контакты | Мы поможем в написании вашей работы! | |
|
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
Дата публикования: 2014-11-18; Прочитано: 273 | Нарушение авторского права страницы | Мы поможем в написании вашей работы!