Форум программистов, компьютерный форум CyberForum.ru
Наши страницы

С++ для начинающих

Войти
Регистрация
Восстановить пароль
 
Рейтинг: Рейтинг темы: голосов - 14, средняя оценка - 4.64
af1283
0 / 0 / 0
Регистрация: 14.01.2009
Сообщений: 5
#1

Погрешность при расчете интеграла. - C++

14.01.2009, 19:56. Просмотров 1680. Ответов 0
Метки нет (Все метки)

Помогите пожалуйста, никак не пойму как вписать погрешность в такую незамысловатую программу для вычисления площади между двумя кривыми.

Программа считает эту самую площадь, но необходимо еще узнать ее с погрешностью 0,01.Есть метод Рунге,заключающийся в том что сначала интеграл считается с шагом h, а потом 2h.
вот формула ((интеграл(h)-интеграл(2h))/3)<погрешность.

Понимаю что должен осуществляться цикл while,в котором должна быть вышестоящая формула.
Дальше не понимаю.
Помогите пожалуйста.

вот мой код.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
#include <iostream.h>
#include <math.h>
#include <windows.h>
 
double func1(double);
double func2(double);
 
 
main()
{
    SetConsoleCP(1251);
   SetConsoleOutputCP(1251);
 
    double arg,warg,ord1,ord,ord2;
 
    double n = 10000;
    double S1 = 0;
    double S2 = 0;
    double h;
    double integ1,integ2;
    double R1,R2;
    
 
 
cout<<"Аргумент\t"<<"Значение"<<endl;
 
//Вычисление точки пересечения графиков
for(arg=0;arg<=1;arg+=0.0000001)
{
    ord1=sin(arg)*sin(arg)+2;
    
    ord2=exp(arg*arg);
    if(fabs(ord2-ord1)<=0.000001)
    {
        warg=arg;
        ord=ord2;
        cout<<arg<<"\t"<<ord<<endl;
    }   
}
 
    h = 2*warg/n;//длина шага
 
    //вычисление площади под экспонентой методом трапеций
    for(int i=1; i<=(n-1); i++)
 
      S1 += func1((-warg) + h*i);
 
    integ1 = h * ( ( func1(-warg) + func1(warg) )/2 + S1 );
 
    cout<<"площадь под экспонентой = "<<integ1<<endl;
 
    
 
 
 
    //вычисление площади под синусоидой методом трапеций
    for(int k=1; k<=(n-1); k++) 
 
 
      S2 += func2((-warg) + h*k);
 
    integ2 = h * ( ( func2(-warg) + func2(warg) )/2 + S2 );
 
    cout<<"площадь под синусоидой = "<<integ2<<endl;
 
 
 
 
    cout<<"площадь фигуры между дугами кривых = "<<integ2-integ1<<endl;
 
return 0;
}
 
double func1(double x )
{
    return exp(x*x);
}
 
double func2(double X)
{
   return sin(X)*sin(X)+2;
}
0
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
14.01.2009, 19:56
Здравствуйте! Я подобрал для вас темы с ответами на вопрос Погрешность при расчете интеграла. (C++):

Большая погрешность при подсчете интеграла методом трапеций - C++
#include &lt;iostream&gt; #include &lt;cmath&gt; using namespace std; double f(double arg) { double fu = cos(arg); return fu; ...

Переполнение при расчете функции - C++
рассчитываю функцию, получается переполнение, что делается в таких случаях? // raschet # include &lt;iostream.h&gt; # include...

Вывод nan при расчете - C++
При расчете codeblocks выводит nan как сделать так, чтобы он считал это ошибкой? #include &lt;iostream&gt; #include &lt;math.h&gt; using...

ошибка при расчете! найти ее - C++
#include &lt;iostream&gt; #include &lt;math.h&gt; using namespace std; int main() { double x; cin &gt;&gt; x; double sum = ((x-2)/(x-1)); ...

Зачем в целом числе писать вещественную часть при расчете триганометрических функций в библиотеке math.h - C++
Доброго времени суток. В чем собственно петрушка. При расчете функции наткнулся на ошибку вычислений. В результате было обнаружено...

Погрешность при вычислениях - C++
Вкратце условие задачи: Есть N счетов, нам нужно соединить их в один так: берется два счета и прибавляются друг к другу и при этом...

Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
14.01.2009, 19:56
Привет! Вот еще темы с ответами:

Погрешность при извлечении квадратного корня - C++
Корень вычисляется с погрешностью из-за этого получается неправильный ответ. Как можно исправить? Вот сам код: #include &lt;iostream&gt; ...

Откуда погрешность? - C++
double j=3.143; cout &lt;&lt; int(j*1000); // выдает 3142 Подскажите откуда берется погрешность?

Погрешность вычислений - C++
Читаю книгу Дейтелов &quot;Как программировать на С++&quot;, попалась следующая задача, где нужно вводить кол-во бензина, пройденный путь, исходя из...

Погрешность интерполяции на C++ - C++
написал код для интерполяционного полинома Ньютоно, он работает правильно, но вот никак не могу разобратса как обчислить его погреность ...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru