4. Используемые технические средства
Для работы программы требуется любой IBM совместимый компьютер, видео адаптор VGA, печатающее устройство. Операционная система MS DOS, драйвер для поддержки русских букв и среда программирования QuickBASIC версии 4.0.
5. Входные и выходные данные
Входными данными для расчета являются:
- исследуемая подынтегральная функция, которая должны быть описана в процедуре-функции f;
a - начальное значение интервала для поиска корней;
b - конечное значение интервала поиска корней;
n - количество расчетных точек;
e - точность расчета коней по аргументу.
Результатами расчета являются:
- таблица аргументов и значений исследуемой функции;
- значение определенного интеграла;
- точность расчета интеграла;
- количество интервалов для которого проведен расчет.
Результаты расчетов приведены в приложении 3 для заданной и тестовой функции.
Список используемой литературы:
1. Вычислительная математика: Учеб. пособие для техникумов /Данилина Н.И., Дубровская Н.С., Кваша О.П., Смирнов Г.Л. - М.: Высш.шк., 1985. - 472 с.
2. Введение в программирование на языке Microsoft BASIC: Учебное пособие / Ю.Я.Максимов, С.В.Осипов, А.В.Потемкин, Щ.С.Симоненков, А.В.Шалашов; Под ред. В.Г.Потемкина. - М.: "Диалог-МИФИ",1991. - 176 с.
3. Жигарев А.Н. и др. Основы компьютерной грамоты / А.Н.Жигарев, Н.В.Макарова, М.А.Путинцева; Под общ. ред. Н.В.Макаровой. -Л.:Машиностроение. Ленингр. отд-ние, 1987.- 255 с.
4. Зельднер Г.А. QuickBASIC для носорога. - М.: ABF, 1994.- 480с.
Приложение 3
|
|
Приложение 1
Приложение 3
|
|
Приложение 1 (продолжение)
Приложение 3
|
|
Приложение 2
DECLARE FUNCTION f! (x!)
p$ = "f = 1 / SQR(2 * pi) * EXP(-x ^ 2 / 2)"
CONST pi = 3.141593
CLS
PRINT "================================================================"
PRINT "| Программа вычисления определенных интегралов |"
PRINT "| модифицированным методом прямоугольников |"
PRINT "----------------------------------------------------------------"
PRINT "| Курсовая работа по дисциплине 'ИНФОРМАТИКА' |"
PRINT "================================================================"
PRINT " Вы изменили подынтегральную функцию y=f(x)? Да - 1, Нет - 2 ";
INPUT u%
IF u% = 2 THEN STOP
INPUT " Начальное значение диапазона изменения аргумента a=", a
INPUT " Конечное значение диапазона изменения аргумента b=", b
INPUT " Количество расчетных точек для построения таблицы n%=", n%
INPUT " Относительная точность вычисления интеграла в % e="; e
DIM x(n%), y(n%)
dx = (b - a) / (n% - 1)
FOR i% = 1 TO n%
x0 = a + px * (i% - 1)
y0 = f(x0)
x(i%) = x0
y(i%) = y0
NEXT
CLS
PRINT " Таблица значений подынтегральной функции"
PRINT p$
PRINT "============================================================="
PRINT "| n/n | Аргумент | Значение | n/n | Аргумент | Значение |"
PRINT "============================================================="
FOR i% = 1 TO n% - 1 STEP 2
PRINT USING "| ### | ###.#### |#######.####| ### | ###.#### |#######.####|"; i%; x(i%); y(i%); i% + 1; x(i% + 1); y(i% + 1) NEXT
PRINT "============================================================="
PRINT " Для продолжения нажмите любую клавишу ";
c$ = INPUT$(1)
CLS
j% = 0
ki% = 1000
DIM i1(ki%)
ee = 1
n1& = n% - 1
DO
j% = j% + 1
s = 0
x = a - dx*0.5
FOR i% = 1 TO n1&
x = x + dx
s = s + f(x)
NEXT
i1(j%) = s * dx
n1& = n1& * 2
px = (b - a) / n1&
IF j% > 1 THEN
ee = (i1(j% - 1) - i1(j%)) / i1(j%) * 100
ee = ABS(ee)
END IF
LOOP UNTIL j% = ki% OR ee <= e
PRINT " Значение определенного интеграла равно "; i1(j%)
PRINT " Относительная точность расчета равна "; ee
PRINT " Количество расчетов для достижения точности "; n1&
PRINT " Вывести на печать таблицу значений и результаты расчета - 1 ";
INPUT v%
IF v% = 1 THEN
LPRINT "================================================================"
LPRINT "| Программа вычисления определеных интегралов |"
LPRINT "| модифицированным методом прямоугольников |"
LPRINT "----------------------------------------------------------------"
LPRINT "| Курсовая работа по дисциплине 'ИНФОРМАТИКА' |"
LPRINT "================================================================"
LPRINT
LPRINT " Таблица значений подынтегральной функции"
LPRINT "============================================================="
LPRINT "| n/n | Аргумент | Значение | n/n | Аргумент | Значение |"
LPRINT "============================================================="
FOR i% = 1 TO n% - 1 STEP 2
LPRINT USING "| ### | ###.#### |#######.####| ### | ###.#### |#######.####|"; i%; x(i%); y(i%); i% + 1; x(i% + 1); y(i% + 1)
NEXT
LPRINT "============================================================="
LPRINT
LPRINT " Значение определенного интеграла равно "; i1(j%)
LPRINT " Относительная точность расчета равна "; ee
LPRINT " Количество расчетов для достижения точности "; n1&
END IF
STOP
END
FUNCTION f (x)
f = 1 / SQR(2 * pi) * EXP(-x ^ 2 / 2)
END FUNCTION
Приложение 3
|
|
Приложение 3