Точность в методе симпсона

[nikozavr]

Вот мой код. Программа считает определенный интеграл методом симпсона.
Как в этот код вставить точность 0.001?

#include <cstdlib>
#include <iostream>
#include <math.h>
#include <windows.h>
char textRus[128];

using namespace std;

float f(float x)
{
return (pow(x,x)*(1+log(x)));
}


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


int i,n; float a,b,x,s2,s1,sym,h;//
CharToOem("Ïðèâåò Console!",textRus);
cout<<textRus<<endl;
cout <<"vvedite a ";cin>>a;
cout <<"vvedite b ";cin>>b;
cout <<"vvedite n ";cin >>n;

h=(b-a)/n;

s1=0;
s2=0;

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

{
x=a+i*h;

if (i%2==0)
s1=s1+f(x);
else
s2=s2+f(x);
}
sym=h/3*(f(a)+f(b)+4*s1+2*s2);




cout<<" sym= "<<sym;
system("PAUSE");
return EXIT_SUCCESS;
}

[nikozavr]

Кто нибудь, АУ?

[Hemul]

используй double переменные

[nikozavr]

Как это? Не понимаю.

[Hemul]

int i,n; 
double a,b,x,s2,s1,sym,h;

double f(double x)

тип данных double - числа с плавающей точкой,которые обрабатываются с двойной точностью (на float выделяется 4 байта, на double 8)

[Obey-Kun]

это не относится к теме вопроса
у float точно и так гораздо выше 0.001

[Vago]

Сообщение от nikozavr

Как в этот код вставить точность 0.001?

Так же, как и в практически любую другую задачу итеративных вычислений:

Полагаем значение_интеграла_на_предыдущем_шаге = 0 ;
[По методу Симпсона] для заданного h вычисляем значение_интеграла_на_следующем_шаге;
while ( | значение_интеграла_на_предыдущем_шаге - значение_интеграла_на_следующем_шаге | > 0.001 ) {
    значение_интеграла_на_предыдущем_шаге = значение_интеграла_на_следующем_шаге ;
    уменьшаем h в два раза ;
    для нового h вычисляем [по методу Симпсона] значение_интеграла_на_следующем_шаге ;
}
Радуемся полученному значению_интеграла_на_следующем_шаге ;

[alexcoder]

В методе Симпсона Погрешность контролируется по правилу Рунге.
http://ru.wikipedia.org/wiki/%D0%9F%...BD%D0%B3%D0%B5
Я делал так:

#include <stdio.h>
#include <math.h>
#include <conio.h>
double f(double x)
{
return sqrt(1*pow(x,5));
}


double simpson(double a,double b,int n)
{
  double h;
   h=(b-a)/n;

   double I, I2 = 0, I4 = 0;
   I4 = f(a+h);
   for(int k = 2; k < n; k += 2 )
   {
       I4 += f(a+(k+1)*h);
       I2 += f(a+k*h);
   }
   I = f(a)+f(b)+4*I4+2*I2;
   I *= h/3;

   return I;

}
double integral(double aa,double bb,int n)
{
double locI = 0.0, xi, h;
int j;
h = double((bb - aa) / n);

xi = aa + h / 2;
for (j=0;j<=n;j++)
{
locI = locI + f(xi);
xi = xi + h;
}
locI = locI * h;
return (locI);
}
int main()
{
double a=0,b=0.5, eps=0.001;
double I1,I2;
int n=1;
do{
I1 = integral(a, b, n);
I2 = integral(a, b, n*2);
printf("n = %d integral = %.5lf\n",n, I2);
n++;
}
while (fabs(I1 - I2) > eps);
printf("integral prug= %.5lf\n", I2);
n=1;
do{
I1 = simpson(a, b, n);
I2 = simpson(a, b, n*2);
printf("n = %d integral simpson = %.5lf\n",n, I2);
n++;
}
while (fabs(I1 - I2) > eps);
printf("integral simpson= %.5lf\n", I2);


getch();
return 0;
}

[nikozavr]

А как это интегрировать в мой код?

[nikozavr]

Ну помогите, А?