Реферат: Расчет двойного интеграла при помощи метода Симпсона - текст реферата. Скачать бесплатно.
Банк рефератов, курсовых и дипломных работ. Много и бесплатно. # | Правила оформления работ | Добавить в избранное
 
 
   
Меню Меню Меню Меню Меню
   
Napishem.com Napishem.com Napishem.com

Реферат

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

Банк рефератов / Математика

Рубрики  Рубрики реферат банка

закрыть
Категория: Реферат
Язык реферата: Русский
Дата добавления:   
 
Скачать
Архив Zip, 12 kb, скачать бесплатно
Обойти Антиплагиат
Повысьте уникальность файла до 80-100% здесь.
Промокод referatbank - cкидка 20%!

Узнайте стоимость написания уникальной работы




/*****************************************************************************

* .FILE : numeric.c

* .TITLE : Расчет двойного интеграла при помощи метода Симпсона (парабол)

* .DESCR : Курсовой проект по численным методам (1994 год)

* :

* : Выполнил: Петренко В.С. (гр. ПС-301)

* : Проверил: Панюков А.В.

* :

* .NOTE : NOT FOR RENTAL OR SALE.

* : FEDERAL LAW PROVIDES SEVERE CIVIL & CRIMINAL PENALTIES FOR

* : UNAUTHORIZED DUPLICATION OR DISTRIBUTION.

* :

* : (C) '94 by P$P

*****************************************************************************/

#include

#include

/*****************************************************************************

* .NAME : m_Simpson

* .TITLE : Расчет интеграла методом Симпсона (парабол)

* .DESCR :

* :

* .PARAMS : double m_Simpson (double (*func) (double, double),

* : double t_fix, double t_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 <= N - 1; i++) sum2 += (*func) (t_fix, -t_limit + (i * h));

sum2 *= 2;


for (i = 1; i <= 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 : double F (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)

{

double integral; /* значение вычисленного интеграла */

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 >= 0; x + y <= П }, 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 => ¦ ");

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 =====> П; k = 1 ");

printf ("\n - ");

printf ("\n ");

printf ("\n 3) f (x, y) = sqrt (x * x + y * y)");

printf ("\n");

do

{

printf ("Ваш выбор: ");

scanf ("%d", &selection);

} while (!(1 <= selection && selection <= 3));


printf ("Параметр k: ");

scanf ("%lg", &k);


do

{

printf ("Число узлов сетки N: ");

scanf ("%d", &N_MAX);

} while (!(N_MAX > 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 FOR RENTAL OR SALE.

* : FEDERAL LAW PROVIDES SEVERE CIVIL & CRIMINAL PENALTIES FOR

* : UNAUTHORIZED DUPLICATION OR DISTRIBUTION.

* :

* : (C) '94 by P$P

*****************************************************************************/

#include


/* выбираемая пользователем функция 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); }




1Авиация и космонавтика
2Архитектура и строительство
3Астрономия
 
4Безопасность жизнедеятельности
5Биология
 
6Военная кафедра, гражданская оборона
 
7География, экономическая география
8Геология и геодезия
9Государственное регулирование и налоги
 
10Естествознание
 
11Журналистика
 
12Законодательство и право
13Адвокатура
14Административное право
15Арбитражное процессуальное право
16Банковское право
17Государство и право
18Гражданское право и процесс
19Жилищное право
20Законодательство зарубежных стран
21Земельное право
22Конституционное право
23Конституционное право зарубежных стран
24Международное право
25Муниципальное право
26Налоговое право
27Римское право
28Семейное право
29Таможенное право
30Трудовое право
31Уголовное право и процесс
32Финансовое право
33Хозяйственное право
34Экологическое право
35Юриспруденция
36Иностранные языки
37Информатика, информационные технологии
38Базы данных
39Компьютерные сети
40Программирование
41Искусство и культура
42Краеведение
43Культурология
44Музыка
45История
46Биографии
47Историческая личность
 
48Литература
 
49Маркетинг и реклама
50Математика
51Медицина и здоровье
52Менеджмент
53Антикризисное управление
54Делопроизводство и документооборот
55Логистика
 
56Педагогика
57Политология
58Правоохранительные органы
59Криминалистика и криминология
60Прочее
61Психология
62Юридическая психология
 
63Радиоэлектроника
64Религия
 
65Сельское хозяйство и землепользование
66Социология
67Страхование
 
68Технологии
69Материаловедение
70Машиностроение
71Металлургия
72Транспорт
73Туризм
 
74Физика
75Физкультура и спорт
76Философия
 
77Химия
 
78Экология, охрана природы
79Экономика и финансы
80Анализ хозяйственной деятельности
81Банковское дело и кредитование
82Биржевое дело
83Бухгалтерский учет и аудит
84История экономических учений
85Международные отношения
86Предпринимательство, бизнес, микроэкономика
87Финансы
88Ценные бумаги и фондовый рынок
89Экономика предприятия
90Экономико-математическое моделирование
91Экономическая теория

 Анекдоты - это почти как рефераты, только короткие и смешные Следующий
На чемпионате мира по самым частым поводам для секса победила дружба.
Anekdot.ru

Узнайте стоимость курсовой, диплома, реферата на заказ.

Обратите внимание, реферат по математике "Расчет двойного интеграла при помощи метода Симпсона", также как и все другие рефераты, курсовые, дипломные и другие работы вы можете скачать бесплатно.

Смотрите также:


Банк рефератов - РефератБанк.ру
© РефератБанк, 2002 - 2017
Рейтинг@Mail.ru