Вычислить приближенное значение интеграла по формулам трапеций и Симпсона

@mkkrzm · Регистрация: 15.06.2015

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

Код ниже.
Слишком сильная разница между двумя формулами. При n=70 в первой формуле выдает 0.03, то во второй 8.3, чего (наверное) не должно быть.
Ткните мне пальцем в ошибку, буду благодарен.

Задание: вычислить приближенное значение интеграла по формулам трапеций и Симпсона для n = 70.

Комментарий модератора

Интересно, для чего нужен редактор формул внизу страницы
$https://www.cyberforum.ru/cgi-bin/latex.cgi?\int_{0}^{2}\sqrt{{(4-{x}^{2})}^{3}}dx$

Вычислить приближенное значение интеграла по формулам трапеций и Симпсона

C++

#include "stdafx.h"
#include <iostream>
#include <locale.h>
#define _USE_MATH_DEFINES
#include <math.h>
using namespace std;
 
void main(){
    setlocale(LC_ALL,".1251");
    float a,b,h,x,s1=0,s2=0,s,s3,rs,res;
    int i,n;
    cout<<"Введите нижний предел:";
    cin>>a;
    cout<<"Введите верхний предел:";
    cin>>b;
    cout<<"\nВведите n=";
    cin>>n;
    x=a;
    h=(b-a)/2/n;
    for (i=1;i>=2*n-2;i+=2)
    {
        s1+=sqrt(pow(4-x*x,3));
        x=2*h;
    }
    s2*=4;
    s=s1+s2+sqrt(pow(4-(a*a),3))+sqrt(pow(4-(b*b),3));
    res=(b-a)*s/(6*n);
    cout<<"Интеграл Симпсона:"<<res<<endl;
    rs=(b-a)/n;
    x=a;
    s3=0;
    for (i=1;i<=n-1;i++)
    {
        s3+=(sqrt(pow(4-(x*x),3)+sqrt(pow(4-(a*a),3))))/2;
        x+rs;
    }
    s3*=rs;
    cout<<"Интеграл Трапеции="<<s3<<endl;
    system("PAUSE");
}

Опции темы