2.4.1.3 Реализация программного модуля.

Форма обращения к модулю

USE source_data

Спецификация идентификаторов параметров модуля

Идентиф.	Тип 18перем.	Инженерное обознач.	Значение	Размерность	Содержательное наименование
Cy	Real				Коэффициент подъемной силы
V	Real	V	250	м/с	Скорость полёта
Cx	Real				Коэффициент лобового сопротивления
G0	Real		120000	Н	Начальный вес самолета
Gk	Real		70000	Н	Конечный вес самолета
Ro	Real.		0,705	кг/м3	плотность воздуха на высоте полета
S	Real	S	40	м2	Площадь крыла

3) Текст модуля

module source_data

real,dimension(5)::Cy=(/0.,0.1,0.2,0.3,0.4/),&

Cx=(/0.018,0.018,0.019,0.021,0.027/)

real,parameter::G0=120000,Gk=70000,V=250,S=40

real,parameter::ro=0.705

end module

2.4.2 Описание подпрограммы «Trap».

2.4.2.1 Общие характеристики.

1) Наименование подпрограммы

интегрирование методом трапеций

2) Назначение модуля

модуль предназначен для численного интегрирования заданной функции в заданных пределах

3) Язык программирования

Фортран Compaq Visual Fortran

4) Тип и версия операционной системы

Microsoft Windows XP

Home Edition

Версия 2002

Service Pack 2

2.4.2.2 Описание модуля.

1)Математическая постановка задачи

2)Описание используемого алгоритма

блок – схема:

см. приложение № 1

Реализация программного модуля

1. Форма обращения к модулю

trap(Gk,G0,n,f)

2)Спецификация идентификаторов входных параметров модуля

Идентиф.	Тип перем.	Инженерное обознач.	Значение	Размерность	Содержательное наименование
G0	Real		120000	H	Нижний предел интегрирования
Gk	Real		70000	H	верхний предел интегрирования
nt	Integer		20		количество шагов

3) Спецификация идентификаторов выходных параметров модуля

Идентиф.	Тип 22перем.	Инженерное обознач.	Значение	Размерность	Содержательное наименование
trap	Real				Значение подынтегрального выражения

Текст модуля

! Подпрограмма интерирования

function trap(Gk,G0,nt,f)

h=(G0-Gk)/(nt-1)

x=Gk

y=G0

s=f(x)+f(y)/2

do i=2,nt-1

x=x+h

s=s+f(x)

end do

trap=s*h

end function

Тестовая задача модуля

1)Описание тестовой задачи

Вычислить интеграл:

a=0

b=1

nt=51

2. описание исходных данных тестовой задачи

3. текст вызывающей программы

program integ

external f

integer, parameter:: N=51

a=0

b=1

S=trap(a,b,N,f)

print*, ‘S=’,S

end

4) текст вызываемых прграмм

Function f(x)

f=X**2

End function

5) Результаты решения тестовой задачи

S= 0.3333998

2.4.3 Описание подпрограммы «f».

2.4.3.1 Общие характеристики.

1) Наименование подпрограммы

Подынтегральная функция

2) Назначение модуля

модуль предназначен для задания подынтегральной функции

3) Язык программирования

Фортран Compaq Visual Fortran

4) Тип и версия операционной системы

Microsoft Windows XP

Home Edition

Версия 2002

Service Pack 2

2.4.3.2 Описание модуля.

Подпрограмма задающая подынтегральную функцию

2.4.3.3 Реализация программного модуля.

Алгоритм.

См. раздел 2.2.3

Блок – схема.

см. приложение №1

Таблица идентификаторов.

Идентиф.	Тип 25перем.	Инженерное обознач.	Значение	Размерность	Содержательное наименование
Cyz	Real				Текущее значение коэффициента подъемной силы
V	Real	V	250	м/с	Скорость полёта
Cxz	Real				Текущее значение коэффициента лобового сопротивления
G	Real			Н	Вес самолета

Листинг.

function f(G)

use source_data

common/cml/ Cxz,Cyz

f=(V*Cyz)/(G*Cxz)

end function

Обращение к подпрограмме.

f(x)

Тестирование подпрограммы.

Подсчет контрольной точки на калькуляторе показал, что расчеты верны:

f (x)= 0,0156 для первой итерации.

2.4.4 Описание подпрограммы «Inter».