Численное интегрирование. Проблема.

[nikozavr]

Вот программа, находит определенный интеграл методом левых прямоугольников.
Вроде все работает, но считает неверно. Кто может подсказать?
Пределы интегрирования a=1 b=3, подинтегральное выражение в коде.

#include <cstdlib>
#include <iostream>
#include <math.h>

using namespace std;



int main(int argc, char *argv[])
{


int i,n; float a,b,x,f,d,s,y,h;
cout <<"vvedite a ";cin>>a;
cout <<"vvedite b ";cin>>b;
cout <<"vvedite n ";cin >>n;

h=(b-a)/n;

s=0;
x=a;

for (i=0;i<=n-1;i++)

{

f = pow(x,x)*(1+log10(x));
s=s+f;
cout<<" f= "<<f;

x=x+h;
}
y=s*h;
cout<<" y= "<<y;
system("PAUSE");
return EXIT_SUCCESS;
}




В матлабе, в интернете и в ехеле ответ 26.
У меня около 18. При любом количестве узлов.

[Алсунчик]

Дан двумерный массив размерностью Nx(в квадрате) увеличить все элемены находящиеся выше главной диагонали в пать раз.Вывести полученный массив.

[nikozavr]

Сообщение от Алсунчик

Дан двумерный массив размерностью Nx(в квадрате) увеличить все элемены находящиеся выше главной диагонали в пать раз.Вывести полученный массив.

Что за бред!

[Д_М]

В матлабе, в интернете и в ехеле ответ 26.
У меня около 18. При любом количестве узлов.

Ошибок не вижу.
Ответ 26 получается, если десятичный логарифм заменить натуральным. Вы не перепутали случаем log и log10 ?

[nikozavr]

Блин, ну я и дятел. Конечно перепутал. Спасибо.

Может еще кто подскажет, как учесть в данной программе точность 0.001?

[Д_М]

Для монотонно возрастающей функции

I1 <= I <= I2

I - точное значение интеграла
I1 - приближение "левыми" прямоугольниками
I2 - приближение "правыми" прямоугольниками

[nikozavr]

Извините не догнал. Как это в коде реализовать?

[Д_М]

...
y=s*h;
float z = (s - f(a) + f(b))*h;
float precision = z -y;
cout<<" y= "<<y << " precision= " << precision << endl;
system("PAUSE");
return EXIT_SUCCESS;

[nikozavr]

А где же тут точность?