Напишите программу приближенного вычисления интеграла методом трапеций.
/*
Программа вычисляет приближенное значение интеграла
методом трапеций.
*/
#include "locale.h"
#include "curses.h"
#include "stdio.h"
#include "stdlib.h"
#define max(a,b) (a>b) ? a : b
main()
{
float a,b;
float dx;
float x;
float y1, y2;
float summa;
float pl;
int n, i;
char *loc;
loc = setlocale(LC_CTYPE, "ru_RU.UTF-8");
initscr();
printw("Приближенное значение интеграла\n");
printw("Введите начало отрезка -> ");
scanw("%f", &a);
printw("Конец отрезка -> ");
scanw("%f", &b);
printw("Прирощение аргумента -> ");
scanw("%f", &dx);
n = (b - a) / dx + 1;
summa = 0;
x = a;
for(i = 0; i <= n; i++)
{
y1 = 5 * x * x - x + 2;
x += dx;
y2 = 5 * x * x - x + 2;
summa += (y2 + y1) * dx / 2;
}
printw("max %5.2f\n", max(a,b));
printw("Значение интеграла равно -> %5.2f", summa);
getch();
endwin();
return 0;
}
Комментариев нет:
Отправить комментарий