VPF::[C++] решение интеграла a-b методом Симпсона

Rogerok

Дата 18.1.2010, 21:02 (ссылка)

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

Загрузка ...

Новичок

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

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

Доброго времени суток
Помогите пожалуйста изменить код.
Данный код ищет "интеграл от 0 до n/2 (1-1/4 sin2xdx) отрезок интегрирования 2n=20"

А мне нужно чтобы просто искал значение интеграла от а до b по методу симпсона
Помогите правильно переписать код.

Код


#include <vcl.h>
#include <math.h>
#include <iostream.h>
#include <conio.h>
#pragma hdrstop
 
using namespace std;
 
//Подынтегральная функция
double func(double x)
{
    double r = sin(x);
    r = 1. - r*r/4.;
    return r;
}
double simpson( double (*Fx)(double), double a, double b, int m);
 
#pragma argsused
int main(int argc, char* argv[])
{
   double I;
   double a,b;
   int m;
 
   cout << "\n Input integral atributes:\n" << endl;
   cout << "\t-> Enter a = ";
   cin >> a;
 
   cout << "\t-> Enter b = ";
   cin >> b;
 
   cout << "\t-> Enter m = ";
   cin >> m;
 
   // вычисление интеграла
 
   I = simpson( &func, a, b, m );
 
   // вывод результата
   cout << endl << "      I = " << I;
   getch();
}
 
//---------------------------------------------------------------------------
// Метод симпсона
double simpson( double (*Fx)(double), double a, double b, int m)
{
 
   double h;
   h=(b-a)/m;
 
   double I, I2 = 0, I4 = 0;
   I4 = Fx(a+h);
   for(int k = 2; k < m; k += 2 )
   {
       I4 += Fx(a+(k+1)*h);
       I2 += Fx(a+k*h);
   }
   I = Fx(a)+Fx(b)+4*I4+2*I2;
   I *= h/3;
 
   return I;
}

Sanaff

Дата 18.1.2010, 22:12 (ссылка)

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

Загрузка ...

Бывалый

Профиль
Группа: Участник
Сообщений: 151
Регистрация: 15.9.2009
Где: г. Северодвинск

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

а почему здесь double r = sin(x);
не так: double r = sin(2*x);

А в остальном - это и есть Симпсона. "от а до b по методу симпсона" Параметры а, b и m задаются пользователем. в принципе нормально.

--------------------

Программист - это локальный бог ©ICQ 373-628-456

Rogerok

Дата 19.1.2010, 19:10 (ссылка)

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

Загрузка ...

Новичок

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

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

спасибо

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