VPF::Интегрирование методом трапеций

karre

Дата 16.5.2010, 19:50 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 1
Регистрация: 16.5.2010

Репутация: нет
Всего: нет

В этом коде при интегрировании методом трапеций требуется точность 10^(-8), но при сравнении ответов, полученных в этой программе и в Mathematica точность есть лишь 10^(-4). Как это исправить?

Код


#include "stdafx.h"
#include <stdio.h>
#include "math.h"


double f(double x)
{
    return x*exp(x)*sin(x);
}
double Integrate(double a, double b, double eps,double (*f)(double))
{
  double x=a;
  double I=0; 
  while (x<b-eps)
  {
    I += (f(x)+f(x+eps))*(b-a)*eps/2;
    x += eps;
  }
  return I;
}
int main(int argc, char* argv[])
{
  printf("I=%f\n",Integrate(1,2,1e-8,f));
  return 0;
 
}

Фантом

Дата 16.5.2010, 20:11 (ссылка)

(нет голосов)

Загрузка ...

Вы это прекратите!

Профиль
Группа: Участник Клуба
Сообщений: 1516
Регистрация: 23.3.2008

Репутация: 4
Всего: 49

Ничего исправлять не надо, программа считает достаточно точно. По-видимому, Вы просто как-то неудачно воспользовались Mathematica'ой.

P.S. Другое дело, что программа просто-таки на редкость неэффективна. Но результат от этого неправильным не становится.

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Для новичков \| Следующая тема »