Расчет двойного интеграла при помощи метода Симпсона

.FILE numeric.c .TITLE Расчет двойного интеграла при помощи методаСимпсона парабол .DESCR Курсовой проект по численным методам 1994год Выполнил Петренко В.С. гр. ПС-301 Проверил Панюков А.В. .NOTE NOT FORRENTAL OR SALE. FEDERALLAW PROVIDES SEVERE CIVIIMINAL PENALTIES FOR UNAUTHORIZED

DUPLIISTRIBUTION. C 94 by P P include lt math.h gt include lt stdio.h gt .NAME m Simpson .TITLE Расчетинтеграла методом Симпсона парабол .DESCR .PARAMS double m Simpson double func double,double , double t fix, doublet limit, int N double func double, double -подынтегральная ф-я double t fix - фиксированный первый аргумент double t limit - верхний пределинтегрирования, нижний равен -t limit int

N - число точек разбиения .RETURN Значение вычисленного интеграла double m Simpson double func double, double , double t fix, double t limit, int N double sum1 0 - double sum2 0 временные переменные double sum3 0 double sum конечный результат double h 2 t limit N шаг сетки int i временная sum1 func t fix, -t limit func t fix, t limit for i 1 i lt N - 1 i sum2 func t fix, -t limit i h sum2 2 for i 1 i lt

N i sum3 func t fix, -t limit i-0.5 h sum3 4 sum sum1 sum2 sum3 sum h 6 sum return sum Глобальные переменные define PI 3.1415926536 число П double k параметр функции - задается пользователем int N MAX число узлов сетки разбиения double currFunc double, double выбраннаяпользователем функция double f1 double x, double y double f2 double x, double y double f3 double x, double y .

NAME doubleF double dummy, double t .TITLE Вычисляет внутренний интеграл G t . См.текст курсового проекта. .DESCR П 2 - перваявычисляемая функция G t П 2 sin t П 2 dt t - - -П 2 где G t currFunc t, tau dtau - -t .PARAMS double F double dummy, double t double dummy - фиктивный первыйаргумент, при вызове этой функции он неиспользуется, т.к. она функция одного аргумента double - действительный второй аргумент .

RETURN Значение функции G t sin k t double F double dummy, double t double G t t PI 2 сдвижка начала координат, чтобы пределы были симметричны в нашем случае - на П 2 G m Simpson currFunc, t, t, N MAX return G sin k t .NAME main .TITLE Основная диалоговая функция. .DESCR Запрашивается интересующая пользователя функция, параметр k и число узлов сетки

N MAX. Выводитна экран вычисленное значение интеграла и два справочных значения - П и П 2. .PARAMS void main void .RETURN void main void doubleintegral значение вычисленногоинтеграла int selection номер выбранной функции массивдоступных функций double functions double, double f1, f2, f3 printf n Вычисление интеграла методомСимпсона парабол printf n printf n printf n I sin k x y f x, y dx dy printf n printf n D printf n где

D x, y x, y gt 0 x y lt П , f Е C D printf n printf nДля какой функции рассчитывать printf n printf n 1 f x, y 0.5 cos y printf n - 0 k ! 1 printf n sin x sin kx dx gt printf n - L- П 2 k 1 printf n printf n 2 f x, y 0.5 - sin y printf n - printf n x sin kx dx gt П k 1 printf n - printf n printf n 3 f x, y sqrt x x y y printf n do printf Ваш выбор scanf d , amp selection while ! 1 lt selection amp amp selection lt 3 printf

Параметр k scanf lg , amp k do printf Число узлов сетки N scanf d , amp N MAX while ! N MAX gt 0 printf n printf n Расчет интеграла currFunc functions selection - 1 текущая функция integral m Simpson F, 0, PI 2, N MAX вычисляем интеграл printf n Значение интеграла равно .12lg , integral вывод printf n

Величины П .12lg П 2 .12lg , PI, PI 2 .FILE func.c .TITLE Содержит функции пользователя, которые можно изменятьбез перекомпиляции основной программы .DESCR После изменения этого модуля его необходимоперекомпилировать ислинковать с numeric.obj .NOTE NOT FORRENTAL OR SALE. FEDERALLAW PROVIDES SEVERE CIVIIMINAL PENALTIES FOR UNAUTHORIZED DUPLICATION

OR DISTRIBUTION. C 94by P P include lt math.h gt выбираемая пользователем функция No.1 double f1 double x, double y return 0.5 cos y выбираемая пользователем функция No.2 double f2 double x, double y return 0.5 - sin y выбираемая пользователем функция No.3 double f3 double x, double y return sqrt x x y y