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

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

* .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 <math.h>

#include <stdio.h>

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

* .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 ");

Внимание, отключите Adblock

Вы посетили наш сайт со включенным блокировщиком рекламы!
Ссылка для скачивания станет доступной сразу после отключения Adblock!

Скачать