Вычислить интеграл

@Ирина0101 · Регистрация: 21.11.2013

Author24 — интернет-сервис помощи студентам

Помогите, пожалуйста, решить следующую задачу:

Вычислить интеграл по бесконечному промежутку с точностью 0.000001, используя правило Рунге, начал.число интервалов выбирать равным 4, при этом сначала свести задачу к вычислению интеграла по конечному отрезку(у меня вычисления происходят от 0 до 15), а затем использовать формулу прямоугольников:

$https://www.cyberforum.ru/cgi-bin/latex.cgi?\int_{0}^{\propto }{e}^{-{x}^{2}}$

Вот мой код, но как я поняла-она выдает неверное решение:

C++

#include <cmath>
#include <iomanip>
#include <iostream>
using namespace std;
double eps=0.000001,e=2.7182818;
double fun(double x)
{
    double result;
    result=pow(e,(-1)*x*x);
    return result;
}
double trap(double a,double b)
{
    return (fun(a)+fun(b))/2*(b-a);
}
bool runge (double i1,double i2)
{
    double ir;
    ir=fabs(i1-i2);
    if (ir<eps) return true;
    return false;
}
 
int main()
{
    int n=4;
    bool mark;
    mark=false;
    double a=0,b=15,h,h2,i1=0,i2=0,a1,a2,b1,b2;
    while(mark==false)
    {
        i1=0,i2=0;
        h=(b-a)/n;
        h2=h/2;
        b1=a;
        b2=a;
        for(int i=0;i<n;i++)
        {
            a1=b1;
            b1=b1+h;
            a2=b2;
            b2=a2+h2;
            i1+=trap(a1,b1);
            i2+=trap(a2,b2);
        }
        mark=runge(i1,i2);
        n+=1;
 
    }
 
    cout<<"n= "<<n<<endl;
    cout<<"Otvet= "<<i1<<endl;
    system("pause");
    return 0; 
}

	04.12.2013, 23:47

Опции темы