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

* : Выводит на экран вычисленное значение интеграла и

* : два справочных значения - П и П/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 - ");