From aa0d37540f5bc4bb6327b43c7213b46d0e261cf0 Mon Sep 17 00:00:00 2001 From: justanothercatgirl Date: Fri, 4 Oct 2024 17:37:14 +0300 Subject: Пока не закончил логику выбора программы, но готовы все алгоритмы MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- src/integral.c | 38 ++++++++++++++++++++++++++++++++++++++ 1 file changed, 38 insertions(+) create mode 100644 src/integral.c (limited to 'src/integral.c') diff --git a/src/integral.c b/src/integral.c new file mode 100644 index 0000000..d2ae11f --- /dev/null +++ b/src/integral.c @@ -0,0 +1,38 @@ +#include "integral.h" + +// эта переменная позволяет контролировать точность интегрирования +// и может быть установлена в файле main +double accrc = 10000.0; + +// интеграл прямоугольниками +double int_rect(f_t f, double a, double b) { + // f(x + dx/2) * dx + const double dx = (b-a) / accrc; + const double dx_2 = dx / 2.0; + double sum = 0.0; + for (double i = a + dx; i < b; i += dx) + sum += f(i - dx_2) * dx; + return sum; +} + +// интеграл трапециями +double int_trap(f_t f, double a, double b) { + // ( f(x) + f(x+dx) )/2 * dx + const double dx = (b-a) / accrc; + double sum = 0.0; + for (double i = a; i < b; i += dx) + sum += dx * (f(i) + f(i+dx))/2 ; + return sum; +} + +// интегрирование параболами +double int_simp(f_t f, double a, double b) { + const double dx = (b-a) / accrc; + const double dx_2 = dx / 2; + const double dx_6 = dx / 6; + double sum = 0.0; + for(double i = a; i < b; i += dx) { + sum += dx_6 * (f(i) + 4 * f(i + dx_2) + f(i + dx)); + } + return sum; +} -- cgit v1.2.3-70-g09d2