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

Интеграл BC++ - C++

Восстановить пароль Регистрация
 
Хакер68
0 / 0 / 0
Регистрация: 08.12.2011
Сообщений: 22
09.09.2012, 20:48     Интеграл BC++ #1
Задание:
Вычислить определённый интеграл
http://www.cyberforum.ru/cgi-bin/latex.cgi?I=\int_{a}^{b}f(x)dx
методом прямоугольников:

http://www.cyberforum.ru/cgi-bin/latex.cgi?I\approx \frac{b-a}{n}({y}_{0}+{y}_{1}+...+{y}_{n}),
где
http://www.cyberforum.ru/cgi-bin/latex.cgi?n - количество отрезков разбиения;
http://www.cyberforum.ru/cgi-bin/latex.cgi?{y}_{0},\ {y}_{1},\ ...,\ {y}_{n} - значения функции на концах отрезков.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
#include <iostream.h>
#include <math.h>
main()
{
double x,a,b,n,y,h,f,k,t;
float x0;xn;y0;yn;
cout<<a;
cin>>a;
cout<<b;
cin>>b;
cout<<k;
cin>>k;
h=(b-a)/n;
cin>>n;
for (i=1; i<=n; i++);
f=(a+h)*(i-0,5);
y=k*f+b;
t=(y-b)/k;
return 0;
Вот по этому сайту делал ТУТ

Добавлено через 1 минуту
У меня много ошибок. Я не знаю как ее делать помогите.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
09.09.2012, 20:48     Интеграл BC++
Посмотрите здесь:

C++ Интеграл
C++ Интеграл
C++ интеграл
интеграл C++
Определенный интеграл C++
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
miriganua
129 / 100 / 4
Регистрация: 05.02.2012
Сообщений: 241
09.09.2012, 21:08     Интеграл BC++ #2
Это вроде работает:
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
#include <iostream>
#include <cmath>
 
using std::cout;
using std::cin;
using std::pow;
 
double function(double x)
{
    return 4 - x * x;
}
 
int main(int argc, char* argv[])
{
    double a;
    cout << "Enter a:";
    cin >> a;
    
    double b;
    cout << "Enter b:";
    cin >> b;
    
    double area = 0;
    for (double x = a; x < b; x += 0.0001)
    {
        area += function(x + 0.00005) * 0.0001;
    }
 
    cout << area << '\n';
 
    return 0;
}
Хакер68
0 / 0 / 0
Регистрация: 08.12.2011
Сообщений: 22
09.09.2012, 22:01  [ТС]     Интеграл BC++ #3
Цитата Сообщение от Хакер68 Посмотреть сообщение
Задание:
Вычислить определённый интеграл
http://www.cyberforum.ru/cgi-bin/latex.cgi?I=\int_{a}^{b}f(x)dx
методом прямоугольников:

http://www.cyberforum.ru/cgi-bin/latex.cgi?I\approx \frac{b-a}{n}({y}_{0}+{y}_{1}+...+{y}_{n}),
где
http://www.cyberforum.ru/cgi-bin/latex.cgi?n - количество отрезков разбиения;
http://www.cyberforum.ru/cgi-bin/latex.cgi?{y}_{0},\ {y}_{1},\ ...,\ {y}_{n} - значения функции на концах отрезков.

C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
#include <iostream.h>
#include <math.h>
main()
{
double x,a,b,n,y,h,f,k,t;
float x0;xn;y0;yn;
cout<<a;
cin>>a;
cout<<b;
cin>>b;
cout<<k;
cin>>k;
h=(b-a)/n;
cin>>n;
for (i=1; i<=n; i++);
f=(a+h)*(i-0,5);
y=k*f+b;
t=(y-b)/k;
return 0;
Вот по этому сайту делал ТУТ

Добавлено через 1 минуту
У меня много ошибок. Я не знаю как ее делать помогите.
Как делать?

Добавлено через 55 секунд
#include <cmath> не проходил мне нужен math.h
miriganua
129 / 100 / 4
Регистрация: 05.02.2012
Сообщений: 241
09.09.2012, 22:53     Интеграл BC++ #4
Цитата Сообщение от Хакер68 Посмотреть сообщение
#include <cmath> не проходил мне нужен math.h
Это одно и то же.
А насчет алгоритма: фактически тебе нужно найти площадь криволинейной трапеции(геометрический смысл определенного интеграла). А метод прямоугольников заключается в разбиении ее на n-ое кол-во прямоугольников, нахождения их площадей и суммирование их.

Для начала тебе нужна сама математическая функция(желательно вынести ее в отдельную функцию). А сам интеграл находим следующим образом(исходя из формулы): в цикле (например, for(double x = a; x <= b; x+= (a+b)*n)) проводим суммирование произведений (a+b)/n умноженное на значение функции в данной точке f(a + (a+b)*2n).

Это я реализовал в своем коде только вместо (a+b)*n я принял 0.0001
Хакер68
0 / 0 / 0
Регистрация: 08.12.2011
Сообщений: 22
14.09.2012, 21:24  [ТС]     Интеграл BC++ #5
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
#include <iostream.h>
#include <math.h>
main()
{
int i, x;
int n;
float a,b,s,y,h;
cout«"\n Input integral:";
cout«"a";
cin»a;
cout«"b";
cin»b;
cout«"n";
cin»n;
cout«"y";
cin»y;
h=(b-a)/n;
x=a;
for(i=0; i<=n-1; i++)
{
x=x+h;
s=h+y*(x+h/2);
}
cout«"\n Result"«s;
return 0;
}
ВОТ Я СДЕЛА КОРОЧЕ. НО ЕСТЬ ОДНО НО. Как мне теперь сюда вписать yi=yi-1(xi-xi-1) и y1=y0(x1-X0) и это еще y2=y1(x2-x1). Все цифры и буква i они находятся внизу формулы. НЕ степень, а наоборот. Как мне в эту задачу это вписать. Цикл нам дан. А как вписать?
-=ЮрА=-
Заблокирован
Автор FAQ
15.09.2012, 08:47     Интеграл BC++ #6
Хакер68 i - означает индекс элемента массива. Всё делается элементарно просто
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
#include <cmath>
#include <iostream>
using namespace std;
 
int main()
{
    //Âîò Г*Г*ГёГҐ ГІГ*áëè÷Г*îå Г§Г*Г¤Г*Г*ГЁГҐ ГґГіГ*êöèè
    double x[] = {-3,-2.2,0,3.4,5.2,6.7,12.8};
    double y[] = {20,-6.7,1,1.2,4.3,-8,3.6};
 
    double s   = 0;//ÑîáñòâåГ*Г*Г® ГЁ áóäåò Г*Г*ГЄГ*ïëèâГ*ГІГј Г§Г*Г*Г·ГҐГ*ГЁГҐ ñóììû
    //ïëîùГ*äåé ýëåìåГ*ГІГ*Г°Г*ûõ ïðÿìîóãîëüГ*èêîâ, 
    //ГІ.ГЄГҐ ГўГЄГ®Г*öå áóäåò ñîäåðæГ*ГІГј Г§Г*Г*Г·ГҐГ*ГЁГҐ ГЁГ*ГІГҐГЈГ°Г*Г«Г*
    
    int i;//Ñ÷¸ò÷èê
    //Îïðåäåëÿåì Г°Г*çìåðГ*îñòü Г¬Г*Г±Г±ГЁГўГ*
    int n = sizeof(x)/sizeof(x[0]);
    //Âû÷èñëåÿì ГЁГ*ГІГҐГЈГ°Г*Г« ìåòîäîì ïðÿìîóãîëüГ*èêîâ
    for(i = 1; i < n - 1; i++)
        s += y[i - 1]*(x[i] - x[i - 1]);
    cout<<"Integral = "<<s<<endl;
    return 0;
}
Отработка http://codepad.org/purx9XdX

PS:Предыдущее задание где не было массивом, правильней делать вот так вычислить интеграл т.е задавать интервал и значение шага, а не брать заранее заданное значение шага (есть случае когда шаг должен быть ещё меньше)
Цитата Сообщение от Игорь Миронюк Посмотреть сообщение
function(x + 0.00005) * 0.0001;
- причём почему умножаем на шаг, либо его половину, а к аргументу x + 0.00005 прибавляем в 5-ть раз большее значение
Миниатюры
Интеграл BC++  
Yandex
Объявления
15.09.2012, 08:47     Интеграл BC++
Ответ Создать тему
Опции темы

Текущее время: 00:14. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru