#include <cmath>
#include <iostream>
#define PI 3.1415926535
/*----подинтегральная функция---*/
double INTEGR(double x)
{
return exp(-(x*x)/2) ;
}
/*----Интегрирование методом трапеции---*/
double trapez(double left, double right, double h)
{
double sum = 0;
double runner;
/*Формула трапеции*/
for(runner = left + h; runner < right; runner += h)
{
sum += INTEGR(runner) ;
}
sum = (sum + 0.5*(INTEGR(left) + INTEGR(right)) ) * h;
return sum * 1 / sqrt(2 * PI);
}
int main(int argc, char ** argv)
{
setlocale(LC_ALL, "Russian");
double a, b;
double h;
std::cout << "Нижнее значение интеграла: " ;
std::cin >> a;
std::cout << "Верхнее значение интеграла: ";
std::cin >> b;
std::cout << "Шаг интегрирования: ";
std::cin >> h;
double trap = trapez(a, b, h);
std::cout <<"Ответ по методу трапеции: " << trap << std::endl;
system("pause");
return 0;
}
I2luY2x1ZGUgPGNtYXRoPgojaW5jbHVkZSA8aW9zdHJlYW0+CiNkZWZpbmUgUEkgMy4xNDE1OTI2NTM1CgovKi0tLS3Qv9C+0LTQuNC90YLQtdCz0YDQsNC70YzQvdCw0Y8g0YTRg9C90LrRhtC40Y8tLS0qLwpkb3VibGUgSU5URUdSKGRvdWJsZSB4KQp7CglyZXR1cm4gZXhwKC0oeCp4KS8yKSA7Cn0KCi8qLS0tLdCY0L3RgtC10LPRgNC40YDQvtCy0LDQvdC40LUg0LzQtdGC0L7QtNC+0Lwg0YLRgNCw0L/QtdGG0LjQuC0tLSovCmRvdWJsZSB0cmFwZXooZG91YmxlIGxlZnQsIGRvdWJsZSByaWdodCwgZG91YmxlIGgpCnsKCWRvdWJsZSBzdW0gPSAwOwoJZG91YmxlIHJ1bm5lcjsKCgkvKtCk0L7RgNC80YPQu9CwINGC0YDQsNC/0LXRhtC40LgqLwoJZm9yKHJ1bm5lciA9IGxlZnQgKyBoOyBydW5uZXIgPCByaWdodDsgcnVubmVyICs9IGgpCgl7CgkJc3VtICs9IElOVEVHUihydW5uZXIpIDsKCX0KCglzdW0gPSAoc3VtICsgIDAuNSooSU5URUdSKGxlZnQpICsgSU5URUdSKHJpZ2h0KSkgKSAqIGg7CgoJcmV0dXJuIHN1bSAqIDEgLyBzcXJ0KDIgKiBQSSk7Cn0KCmludCBtYWluKGludCBhcmdjLCBjaGFyICoqIGFyZ3YpCnsKCXNldGxvY2FsZShMQ19BTEwsICJSdXNzaWFuIik7CgoJZG91YmxlIGEsIGI7Cglkb3VibGUgaDsKCglzdGQ6OmNvdXQgPDwgItCd0LjQttC90LXQtSDQt9C90LDRh9C10L3QuNC1INC40L3RgtC10LPRgNCw0LvQsDogIiA7CglzdGQ6OmNpbiA+PiBhOwoJc3RkOjpjb3V0IDw8ICLQktC10YDRhdC90LXQtSDQt9C90LDRh9C10L3QuNC1INC40L3RgtC10LPRgNCw0LvQsDogIjsKCXN0ZDo6Y2luID4+IGI7CglzdGQ6OmNvdXQgPDwgItCo0LDQsyDQuNC90YLQtdCz0YDQuNGA0L7QstCw0L3QuNGPOiAiOwoJc3RkOjpjaW4gPj4gaDsKCglkb3VibGUgdHJhcCA9IHRyYXBleihhLCBiLCBoKTsKCXN0ZDo6Y291dCA8PCLQntGC0LLQtdGCINC/0L4g0LzQtdGC0L7QtNGDINGC0YDQsNC/0LXRhtC40Lg6ICIgPDwgdHJhcCA8PCBzdGQ6OmVuZGw7IAoJc3lzdGVtKCJwYXVzZSIpOwoJcmV0dXJuIDA7Cn0=