Составители:
Рубрика:
26
6.
Макет печати (по C-программе)
Как видно, результаты ручного и машинного расчётов совпа-
дают.
7. Листинг Fortran-программы
Программный код
! Табулирование функции, зависящей от параметра
PROGRAM
PARAMETR
! Объявление типа и задание данных
INTEGER
:: K = 0
REAL
:: A,A1 = 0.5,A2 = 1,XCP = 0.4
REAL
:: X,Y,XMIN = 0,XMAX = 1,DX = 0.2
X
= XMIN
WRITE
(*,*)' ТАБЛИЦА' ! Вывод заголовка
DO
WHILE (X <= XMAX + DX/3) ! Начало цикла
! Выбор параметра A
IF (X <= XCP) THEN ! Если X ≤ XCP
A = A1
ELSE ! Иначе, если X > XCP
A = A2
END IF
Y = A*X - COS(A*X)**2 ! Вычисление значения Y
! Вывод значений
X,Y
WRITE(*,2) X,Y
2 FORMAT
(1X,' X = ',F3.1,2X,' Y = ',F7.4)
IF(Y < 0) K = K + 1 ! Вычисление K
X = X + DX
END DO
! Конец цикла
WRITE
(*,"( A,I2)")'Результат K = ',K ! Вывод K
END
Результаты представлены в макете печати в п.6.
6. Макет печати (по C-программе) Как видно, результаты ручного и машинного расчётов совпа- дают. 7. Листинг Fortran-программы Программный код ! Табулирование функции, зависящей от параметра PROGRAM PARAMETR ! Объявление типа и задание данных INTEGER :: K = 0 REAL :: A,A1 = 0.5,A2 = 1,XCP = 0.4 REAL :: X,Y,XMIN = 0,XMAX = 1,DX = 0.2 X = XMIN WRITE(*,*)' ТАБЛИЦА' ! Вывод заголовка DO WHILE (X <= XMAX + DX/3) ! Начало цикла ! Выбор параметра A IF (X <= XCP) THEN ! Если X ≤ XCP A = A1 ELSE ! Иначе, если X > XCP A = A2 END IF Y = A*X - COS(A*X)**2 ! Вычисление значения Y ! Вывод значений X,Y WRITE(*,2) X,Y 2 FORMAT (1X,' X = ',F3.1,2X,' Y = ',F7.4) IF(Y < 0) K = K + 1 ! Вычисление K X = X + DX END DO ! Конец цикла WRITE(*,"( A,I2)")'Результат K = ',K ! Вывод K END Результаты представлены в макете печати в п.6. 26
Страницы
- « первая
- ‹ предыдущая
- …
- 24
- 25
- 26
- 27
- 28
- …
- следующая ›
- последняя »