Вычисление интеграла - C++

@GRANDEATH · Регистрация: 13.09.2009

Считаю интеграл на промедутке от -1 до 3... Интеграл a*a*a*a. Точность вычисления 0.00001... Код выполняется 10-12 секунд. Это нормально или можно оптимизировать? Метод вычисления прямоугольники или трапеции ( в принципе без разницы.. И там и там где-то 10 сек. ) На других функциях ( y = sin x, y = x) менее чем за 150 миллисекунд. Что можете сказать? Приложу код на С++, но это не важно, важен алгоритм...

Код

C++

double IntRect(fptr f, double xn, double xk, double eps,long int* n, bool * br)
{
    
    if ( xk <= xn) 
    {
        cout << "Err";
        return 0;
    }
    double S1 = 0, S2;
    double delta = DELTA_PRESET;
 
    (*n) = 0; // Обнуляем число шагов
    DELTA_DEC  = 0.5; //Коэффициент уменьшения  
    do
    {
        S2 = S1;
        S1 = 0;
        
        (*n)++;
        
        //double delta = ( xk - xn ) / eps;
        for (double i = xn; i < xk + delta; i += delta)
        {
            if ( i > xk)
            {
                S1 = S1 * delta; // Умножаем сумму высот на дельта
                S1 = S1 + ( delta - (i - xk) ) * ( f( (i - delta + xk) /2.0)   ); // Находим интеграл последнего отрезка
                break;
            }
            S1 = S1 + f(i + delta / 2.0);
            
        }
        if ( (clock() - t1) > MAX_TIME)
        {
                (*br) = true; // Выполнение прервалось
                break;
        }
        delta *= DELTA_DEC; // Уменьшаем  дельту в DELTA_DEC раза
        
 
    } while ( fabs(S1-S2) > eps);
    return ( (S1+S2) / 2);
}

@GRANDEATH · 10.10.2009, 14:19 **[ТС]**

Спасибо за код, но она все равно считает не 2 секудны..У меня тоже 11 шагов.... Начальное дельта... Ты думаешь его стоит сделать поменьше?
Уменьшил до 0.000001... Хуже стало.. Увеличил - никак.. Практически все тоже..

@odip · 10.10.2009, 14:24

sum1=48.7999940918 count=12 sec1=2 errcode=0

У меня на компе считает 2: sec1=2

Начальное не влияет.
Все время счета уходит на маленькие delta !
Честно говоря не знаю что делать. Пробовал делать delta*= 0.6;
delta*= 0.7;
Тогда вообще за 0 сек считает.
Но подбирать delta для каждой задачи это как-то странно

Суть в том, что мы имеем плохую сходимость.
И кривой метод.

Думаю можно ничего не делать - если объяснишь почему медленно считает

Кстати возможно что метод трапеций будет лучше.

@GRANDEATH · 10.10.2009, 14:34 **[ТС]**

Я вот тоже думал что в сходимости проблема... Такое предчувствие, что надо было задавать уменьшение в зависимости от разницы между fabs(S1 - S2) и eps... Не да ладно... У меня код выполняется 4 секунды... Трапеции вроде тоже не лучше... Ай нет... Хуже... По моему начальному коду - первый метод - 3.8 сек... Трапеции (2-й метод) 7.5 секунд.. Ума не приложу, где подвох...

@odip · 10.10.2009, 14:59

http://ru.wikipedia.org/wiki/Численное_интегрирование
Тут написано что все плохо. И есть другие методы - лучшей точности при том же кол-ве вычислений f(x) !
Но если тебе нужно считать именно заданными методами (прямоугольников и трапеций) - то ничего сделать нельзя.

@GRANDEATH · 10.10.2009, 15:03 **[ТС]**

Ну вообщем понятно)) Спс за разьяснения. По крайней мере будет уверенность в том что написал нормально)))

@Sonnick · 27.10.2009, 16:59

люди плиз помогите, нужно сделать прогу чтоб решала интеграл от синуса, косинуса и x в любой степени, и хотя бы 2 методами(прописать чтоб выбрать метод) трапеции и Гаусса, мне дали текст карявой проги пытался исправить(с иксом не правильно высчитывает), половину проблем убрал, но не могу понять по какому он методу сейчас считает,

C

#include "stdafx.h"
#include "iostream"
using namespace std;
#include "math.h"
double integral(double a,double b,double (*f)(double x)); \\поясните что это?
double fun(double a,double b, double x); \\и эт тоже
int n;
double x;
 
 
 
void main()
{
    double a,b,y;
    printf("\n programma integral!");
    printf("\n vvedite znachenie verhney granicy a ");
    scanf("%lf",&a);
    printf("\n vvedite znachenie nizney granicy b ");
    scanf("%lf",&b);
    printf("\n vyberite functsiu ");
    printf("\n 1-sin, \n 2-cos, \n 3-x^z ");
    scanf("%d",&n);
    printf("\n Введите степень х\n ");
    scanf("%d",&x);
    if (n==3)
    {
        double pow(x,z);
    }
    switch (n)
    {
    case 1: y=integral(a,b,sin); break;
    case 2: y=integral(a,b,cos); break;
    case 3: y=fun(a,b,x); break;
    default: printf("\n neverniy vybor "); y=false;]
    }
    printf("\n resultat = %lf",y);
    scanf("%lf",&b);
}
double integral(double a,double b,double (*f)(double x))
{
    return (f(a)+f(b))/2*(b-a);
}
double fun(double a,double b, double x)
{
    return 0;
}

@GRANDEATH · 27.10.2009, 20:35 **[ТС]**

У меня - все считает...4 функции, каждая двумя методами. Немного переделай. И все) Код я выкладывал. См эту тему стр 2 вроде..

@Sonnick · 27.10.2009, 20:42

а чего переделать?)))я имею ввиду что именно?)

@GRANDEATH · 27.10.2009, 21:17 **[ТС]**

C++

#include <math.h>
#include <iostream>
#include <iomanip>
#include <time.h>
 
 
#define DELTA_PRESET 0.0001 // ГЌГ*Г·Г*Г«ГјГ*Г»Г© Г®ГІГ°ГҐГ§Г®ГЄ
#define TIME_BORDER 300 // Г‚Г°ГҐГ¬Гї Гў Г¬ГЁГ«Г«ГЁГ±ГҐГЄГіГ*Г¤Г*Гµ Г¤Г«Гї ГўГ»ГўГ®Г¤Г* ГЇГ°ГҐГ¤ГіГЇГ°ГҐГ¦Г¤ГҐГ*ГЁГї...
#define MAX_TIME 10000 //ГЊГ*ГЄГ±ГЁГ¬Г*Г«ГјГ*Г®ГҐ ГўГ°ГҐГ¬Гї ГўГ»ГЇГ®Г«ГҐГ*ГЁГї
double  DELTA_DEC  = 0.7; //ГЉГ®ГЅГґГґГЁГ¶ГЁГЅГ*ГІ ГіГ¬ГҐГ*ГјГёГҐГ*ГЁГї
 
clock_t t1, t2; // ГЏГҐГ°ГҐГ¬ГҐГ*Г*Г»ГҐ Г¤Г«Гї ГµГ°Г*Г*ГЁГҐГ*ГЁГї ГўГ°ГҐГ¬ГҐГ*ГЁ
 
using namespace std;
 
// Г‘ГІГ°ГіГЄГІГіГ°Г*, ГЈГ¤ГҐ Г±Г®ГµГ°Г*Г*ГҐГ*Г» Г§Г*Г*Г·ГҐГ*ГЁГї
struct dataS
{
    double method1; // Г‡Г*Г*Г·ГҐГ*ГЁГҐ Г¬ГҐГІГ®Г¤Г* 1
    double method2; // Г‡Г*Г*Г·ГҐГ*ГЁГҐ Г¬ГҐГІГ®Г¤Г* 2
    double value;  // Г€Г±ГІГЁГ*Г*Г®ГҐ Г§Г*Г*Г·ГҐГ*ГЁГҐ ( ГЄГ®ГІГ®Г°Г®ГҐ Г¤Г®Г«Г¦Г*Г» ГЇГ®Г«ГіГ·ГЁГІГј)
    long int n1; // Г—ГЁГ±Г«Г® Г®ГІГ°ГҐГ§ГЄГ®Гў, ГЇГ°ГЁ ГЄГ®ГІГ®Г°Г»Гµ Г¤Г®Г±ГІГЁГЈГ*ГіГІГ* ГІГ®Г·Г*Г®Г±ГІГј Г¤Г«Гї ГЇГҐГ°ГўГ®ГЈГ® Г¬ГҐГІГ®Г¤Г*
    long int n2; // Г—ГЁГ±Г«Г® Г®ГІГ°ГҐГ§ГЄГ®Гў, ГЇГ°ГЁ ГЄГ®ГІГ®Г°Г»Гµ Г¤Г®Г±ГІГЁГЈГ*ГіГІГ* ГІГ®Г·Г*Г®Г±ГІГј Г¤Г«Гї ГўГІГ®Г°Г®ГЈГ® Г¬ГҐГІГ®Г¤Г*
    double eps; // Г’Г®Г·Г*Г®Г±ГІГј ГўГ»Г·ГЁГ±Г«ГҐГ*ГЁГї
    long process_time_1;
    long process_time_2;
    bool break_1;
    bool break_2;
};
 
// Г”ГіГ*ГЄГ¶ГЁГЁ, ГЁГ*ГІГҐГЈГ°Г*Г«Г» ГЄГ®ГІГ®Г°Г»Гµ Г*Г*Г¤Г® ГўГ»Г·ГЁГ±Г«ГЁГІГј
double f1(double x)
{
    return (x);
};
double f2(double x)
{
    return (sin( 22 * x ));
};
double f3(double x)
{
    return ( x*x*x*x );
};
double f4(double x)
{
    return (atan(x));
};
 
 
// Г”ГіГ*ГЄГ¶ГЁГЁ, ГЄГ®ГІГ®Г°Г»ГҐ ГўГ»Г®Г§ГўГ°Г*Г№Г*ГѕГІ ГІГ®Г·Г*Г®ГҐ Г§Г*Г*Г·ГҐГ*ГЁГҐ ГЁГ*ГІГҐГЈГ°Г*Г«Г*
double f1_res(double a, double b)
{
    return ( (b*b-a*a)/2.0 );
};
double f2_res(double a, double b)
{
    return (  (cos(a*22.0)-cos(b*22.0) )/22.0 );
};
double f3_res(double a, double b)
{
    return (  (b*b*b*b*b-a*a*a*a*a)/5.0 );
};
double f4_res(double a, double b)
{
    return ( b*atan(b)-a*atan(a)-(log(b*b+1)-log(a*a+1))/2.0 );
};
 
 
 
// ГЊГ*Г±Г±ГЁГў Г*Г*ГёГЁГµ ГґГіГ*ГЄГ¶ГЁГ©
typedef double (* fptr) (double );
typedef double (* fptr_res) (double, double );
 
 
fptr f_array[] = {f1 , f2 , f3 , f4}; // ГЊГ*Г±Г±ГЁГў ГіГЄГ*Г§Г*ГІГҐГ«ГҐГ© Г*Г* ГґГіГ*ГЄГ¶ГЁГЁ
fptr_res f_array_res[] = {f1_res , f2_res , f3_res , f4_res}; // ГЊГ*Г±Г±ГЁГў ГіГЄГ*Г§Г*ГІГҐГ«ГҐГ© Г*Г* ГІГ®Г·Г*Г®ГҐ Г§Г*Г*Г·ГҐГ*ГЁГҐ ГґГіГ*ГЄГ¶ГЁГ©
 
 
 
 
 
 
 
/*
f - ГЁГ¬Гї ГЁГ*ГІГҐГЈГ°ГЁГ°ГіГҐГ¬Г®Г© ГґГіГ*ГЄГ¶ГЁГЁ,
xn, xk – ГЈГ°Г*Г*ГЁГ¶Г» ГЁГ*ГІГҐГ°ГўГ*Г«Г* ГЁГ*ГІГҐГЈГ°ГЁГ°Г®ГўГ*Г*ГЁГї,
eps – ГІГ°ГҐГЎГіГҐГ¬Г*Гї ГІГ®Г·Г*Г®Г±ГІГј,
n – Г·ГЁГ±Г«Г® Г®ГІГ°ГҐГ§ГЄГ®Гў, ГЇГ°ГЁ ГЄГ®ГІГ®Г°Г®Г¬ Г¤Г®Г±ГІГЁГЈГ*ГіГІГ* ГІГ°ГҐГЎГіГҐГ¬Г*Гї ГІГ®Г·Г*Г®Г±ГІГј (ГўГ»ГµГ®Г¤Г*Г®Г©). 
*/
 
double IntRect(fptr f, double xn, double xk, double eps,long int* n, bool * br)
{
    
    if ( xk <= xn) 
    {
        cout << "Err";
        return 0;
    }
    double S1 = 0, S2;
    double delta = DELTA_PRESET;
 
    (*n) = 0; // ГЋГЎГ*ГіГ«ГїГҐГ¬ Г·ГЁГ±Г«Г® ГёГ*ГЈГ®Гў
    //DELTA_DEC  = 0.5; //ГЉГ®ГЅГґГґГЁГ¶ГЁГҐГ*ГІ ГіГ¬ГҐГ*ГјГёГҐГ*ГЁГї   
    do
    {
        S2 = S1;
        S1 = 0;
        
        (*n)++;
        
        //double delta = ( xk - xn ) / eps;
        for (double i = xn; i < xk + delta; i += delta)
        {
            if ( i > xk)
            {
                S1 = S1 * delta; // Г“Г¬Г*Г®Г¦Г*ГҐГ¬ Г±ГіГ¬Г¬Гі ГўГ»Г±Г®ГІ Г*Г* Г¤ГҐГ«ГјГІГ*
                S1 = S1 + ( delta - (i - xk) ) * ( f( (i - delta + xk) /2.0)   ); // ГЌГ*ГµГ®Г¤ГЁГ¬ ГЁГ*ГІГҐГЈГ°Г*Г« ГЇГ®Г±Г«ГҐГ¤Г*ГҐГЈГ® Г®ГІГ°ГҐГ§ГЄГ*
                break;
            }
            S1 = S1 + f(i + delta / 2.0);
 
            
        }
        if ( (clock() - t1) > MAX_TIME)
        {
                (*br) = true; // Г‚Г»ГЇГ®Г«Г*ГҐГ*ГЁГҐ ГЇГ°ГҐГ°ГўГ*Г«Г®Г±Гј
                break;
        }
        delta *= DELTA_DEC; // Г“Г¬ГҐГ*ГјГёГ*ГҐГ¬  Г¤ГҐГ«ГјГІГі Гў DELTA_DEC Г°Г*Г§Г*
        cout << "#" << *n <<" Delta = " << delta << endl;
 
 
    } while ( fabs(S1-S2) > eps);
    return ( (S1+S2) / 2);
}
 
 
int print(dataS* data)
{
    cout << endl;
    cout << setw(20) << left <<"Г‡Г*Г*Г·ГҐГ*ГЁГҐ 1:" << setw(12) << right << fixed << data->method1 << endl;
    cout << setw(20) << left <<"Г—ГЁГ±Г«Г® ГёГ*ГЈГ®Гў 1:" << setw(12) << right << data->n1 << endl;
    cout << setw(20) << left <<"Г‚Г°ГҐГ¬Гї ГўГ»ГЇГ®Г«ГҐГ*ГЁГї 1-ГЈГ®" << setw(12) << right << data->process_time_1;
    if (data->process_time_1 > TIME_BORDER) cout << "  !!!";
    if (data->break_1)  cout << "  break";
    cout<< endl;
 
    cout << setw(20) << left <<"Г‡Г*Г*Г·ГҐГ*ГЁГҐ 2:" << setw(12) << right << fixed << data->method2 << endl;
    cout << setw(20) << left <<"Г—ГЁГ±Г«Г® ГёГ*ГЈГ®Гў 2:" << setw(12) << right << data->n2 << endl;
    cout << setw(20) << left <<"Г‚Г°ГҐГ¬Гї ГўГ»ГЇГ®Г«ГҐГ*ГЁГї 2-ГЈГ®" << setw(12) << right  << data->process_time_2;
    if (data->process_time_2 > TIME_BORDER) cout << "  !!!";
    if (data->break_2)  cout << "  break";
    cout<< endl;
 
    cout << setw(20) << left <<"Г’Г®Г·Г*Г®Г±ГІГј" << setw(12) << right << data->eps << endl;
    cout << setw(20) << left <<"Г’Г®Г·Г*Г®ГҐ Г§Г*Г*Г·ГҐГ*ГЁГҐ" << setw(12) << right  << fixed <<data->value << endl;
 
    
    return 0;
}
 
double IntSum(fptr f, double xn, double xk, double eps,long int* n, bool *br)
{
 
    //DELTA_DEC  = 0.5; //ГЉГ®ГЅГґГґГЁГ¶ГЁГЅГ*ГІ ГіГ¬ГҐГ*ГјГёГҐГ*ГЁГї   
    if ( xk <= xn) 
    {
        cout << "Err";
        return 0;
    }
 
    double S1 = 0, S2;
    double f1,f2;
    double delta = DELTA_PRESET;
    (*n) = 0; // ГЋГЎГ*ГіГ«ГїГҐГ¬ Г·ГЁГ±Г«Г® ГёГ*ГЈГ®Гў
    do
    {
        (*n)++;
        S2 = S1;
        S1 = 0;
        
        //double delta = ( xk - xn ) / eps;
        for (double i = xn; i < xk + delta; i += delta)
        {
            if ( i > xk)
            {
                S1 = S1 * delta; // Г“Г¬Г*Г®Г¦Г*ГҐГ¬ Г±ГіГ¬Г¬Гі ГўГ»Г±Г®ГІ Г*Г* Г¤ГҐГ«ГјГІГ*
                S1 = S1 + ( delta - (i - xk) ) * ( f(i-delta) + f(xk)) / 2.0; // ГЌГ*ГµГ®Г¦ГЁГ¬ ГЁГ*ГІГҐГЈГ°Г*Г« ГЇГ®Г±Г«ГҐГ¤Г*ГҐГЈГ® Г®ГІГ°ГҐГ§ГЄГ*
                break;
            }
            f1 = f(i);
            f2 = f(i + delta);
            S1 = S1 +  ( ( f1 + f2) / 2.0 );
            
 
        }
        if ( (clock() - t1) > MAX_TIME)
        {
                if ( fabs(S1-S2) < eps) break; // Г‚Г±ГҐ Г*Г®Г°Г¬Г*Г«ГјГ*Г®
 
                (*br) = true; // Г‚Г»ГЇГ®Г«Г*ГҐГ*ГЁГҐ ГЇГ°ГҐГ°ГўГ*Г«Г®Г±Гј
                break;
        }
        delta *= DELTA_DEC; // Г“Г¬ГҐГ*ГјГёГ*ГҐГ¬  Г¤ГҐГ«ГјГІГі Гў DELTA_DEC Г°Г*Г§Г*
        cout << "#" << *n <<" Delta = " << delta << endl;
 
 
    } while ( fabs(S1-S2) > eps);
    cout << endl;
 
    return ( (S1+S2) / 2);
 
 
};

Это funct.h

C++

// Lab_7.cpp : Defines the entry point for the console application.
//
 
#include "stdafx.h"
#include "funct.h"
#include <Windows.h>
 
using namespace std;
 
int _tmain(int argc, _TCHAR* argv[])
{
 
    setlocale(0,"russian");
    cout << "Г‚ГўГҐГ¤ГЁГІГҐ ГІГ®Г·Г*Г®Г±ГІГј ГўГ»Г·ГЁГ±Г«ГҐГ*ГЁГ©" << endl;
    
    double eps = 0;
    cin >> eps;
    
    
    dataS *data = new dataS[4]; // ГЊГ*Г±Г±ГЁГў Г§Г*Г*Г·ГҐГ*ГЁГ©
    
    // Г‚Г»Г·ГЁГ±Г«ГїГҐГ¬ Г§Г*Г*Г·ГҐГ*ГЁГї ГЁГ*ГІГҐГЈГ°Г*Г«Г®Гў
    for (int i = 0; i < 4; i++)
    {
        
        
        data[i].break_1 = false; //// Г‚Г»ГЇГ®Г«Г*ГҐГ*ГЁГҐ Г*ГҐ ГЇГ°ГҐГЄГ°Г*Г№ГҐГ*Г®
        data[i].n1 = 0;
        Sleep(1);
        t1 = clock();
        data[i].method1 = IntRect(f_array[i], -1, 3, eps, &(data[i].n1), &(data[i].break_1));
        t2 = clock();
        data[i].process_time_1 = t2 - t1;
 
        
        data[i].break_2 = false; // Г‚Г»ГЇГ®Г«Г*ГҐГ*ГЁГҐ Г*ГҐ ГЇГ°ГҐГЄГ°Г*Г№ГҐГ*Г®
        data[i].n2 = 0;
        Sleep(1);
        t1 = clock();
        data[i].method2 = IntSum(f_array[i], -1, 3, eps, &(data[i].n2), &(data[i].break_2));
        t2 = clock();
        data[i].process_time_2 = t2 - t1;
 
        data[i].eps = eps; 
        data[i].value = f_array_res[i](-1,3);
    }
 
 
    //cout << IntRect(f_array[0], 01, 0,1, NULL);
 
    for (int i = 0; i < 4; i++)
    {
        cout << endl << endl << endl;
        cout << "Funct#" << i << endl;
        print(&(data[i]));
        
    }
 
 
    // ГЋГ·ГЁГ±ГІГЄГ* ГЇГ*Г¬ГїГІГЁ
    delete data;
    return 0;
}

- Это сам cpp

Добавлено через 4 минуты

C++

 
// Функции, интегралы которых надо вычислить
double f1(double x)
{
        return (x);
};
double f2(double x)
{
        return (sin( 22 * x ));
};
double f3(double x)
{
        return ( x*x*x*x );
};
double f4(double x)
{
        return (atan(x));
};

- это 4 функции, интегралы которых считает код.. Указывай вместо них свои...

C++

// Г”ГіГ*ГЄГ¶ГЁГЁ, ГЄГ®ГІГ®Г°Г»ГҐ ГўГ»Г®Г§ГўГ°Г*Г№Г*ГѕГІ ГІГ®Г·Г*Г®ГҐ Г§Г*Г*Г·ГҐГ*ГЁГҐ ГЁГ*ГІГҐГЈГ°Г*Г«Г*
double f1_res(double a, double b)
{
    return ( (b*b-a*a)/2.0 );
};
double f2_res(double a, double b)
{
    return (  (cos(a*22.0)-cos(b*22.0) )/22.0 );
};
double f3_res(double a, double b)
{
    return (  (b*b*b*b*b-a*a*a*a*a)/5.0 );
};
double f4_res(double a, double b)
{
    return ( b*atan(b)-a*atan(a)-(log(b*b+1)-log(a*a+1))/2.0 );
};

Это точное значение интеграла... Указывай точное значение F(x) = интеграл f(x)
то есть пусть надо найти интеграл от y = x*x. В первой

C++

double f1(double x)
{
        return (x*x); // Это х квадрат
};

Во второй

C++

double f1_res(double a, double b)
{
    return ( (b*b*b-a*a*a)/3.0 ); // Это определенный интеграл от x*x
};

@DevilVergil · 31.03.2010, 01:03

GRANDEATH, в принципе, выполнение кода может занять 10-12 сек., в зависимости от того, как написать функцию. Хочешь оптимизировать - попробуй поменять класс памяти для наиболее значимых элементов, но при такой точности...Ведь обращение идет не к одной функции. При твоем подходе вряд ли что то поможет, нужно коренным образом перестроить программу.

@GRANDEATH 39 / 39 / 1 Регистрация: 13.09.2009 Сообщений: 108
	10.10.2009, 14:19 [ТС]	#16
	Спасибо за код, но она все равно считает не 2 секудны..У меня тоже 11 шагов.... Начальное дельта... Ты думаешь его стоит сделать поменьше? Уменьшил до 0.000001... Хуже стало.. Увеличил - никак.. Практически все тоже.. 0

@odip 7161 / 3223 / 58 Регистрация: 17.06.2009 Сообщений: 14,164
	10.10.2009, 14:24	#17
	sum1=48.7999940918 count=12 sec1=2 errcode=0 У меня на компе считает 2: sec1=2 Начальное не влияет. Все время счета уходит на маленькие delta ! Честно говоря не знаю что делать. Пробовал делать delta= 0.6; delta= 0.7; Тогда вообще за 0 сек считает. Но подбирать delta для каждой задачи это как-то странно Суть в том, что мы имеем плохую сходимость. И кривой метод. Думаю можно ничего не делать - если объяснишь почему медленно считает Кстати возможно что метод трапеций будет лучше. 1

@GRANDEATH 39 / 39 / 1 Регистрация: 13.09.2009 Сообщений: 108
	10.10.2009, 14:34 [ТС]	#18
	Я вот тоже думал что в сходимости проблема... Такое предчувствие, что надо было задавать уменьшение в зависимости от разницы между fabs(S1 - S2) и eps... Не да ладно... У меня код выполняется 4 секунды... Трапеции вроде тоже не лучше... Ай нет... Хуже... По моему начальному коду - первый метод - 3.8 сек... Трапеции (2-й метод) 7.5 секунд.. Ума не приложу, где подвох... 0

@odip 7161 / 3223 / 58 Регистрация: 17.06.2009 Сообщений: 14,164
	10.10.2009, 14:59	#19
	http://ru.wikipedia.org/wiki/Численное_интегрирование Тут написано что все плохо. И есть другие методы - лучшей точности при том же кол-ве вычислений f(x) ! Но если тебе нужно считать именно заданными методами (прямоугольников и трапеций) - то ничего сделать нельзя. 1

@GRANDEATH 39 / 39 / 1 Регистрация: 13.09.2009 Сообщений: 108
	10.10.2009, 15:03 [ТС]	#20
	Ну вообщем понятно)) Спс за разьяснения. По крайней мере будет уверенность в том что написал нормально))) 0

@GRANDEATH 39 / 39 / 1 Регистрация: 13.09.2009 Сообщений: 108
	27.10.2009, 20:35 [ТС]	#22
	У меня - все считает...4 функции, каждая двумя методами. Немного переделай. И все) Код я выкладывал. См эту тему стр 2 вроде.. 0

@Sonnick 0 / 0 / 0 Регистрация: 18.11.2015
	27.10.2009, 20:42	#23
	а чего переделать?)))я имею ввиду что именно?) 0

@DevilVergil 0 / 0 / 0 Регистрация: 27.03.2010 Сообщений: 1
	31.03.2010, 01:03	#25
	GRANDEATH, в принципе, выполнение кода может занять 10-12 сек., в зависимости от того, как написать функцию. Хочешь оптимизировать - попробуй поменять класс памяти для наиболее значимых элементов, но при такой точности...Ведь обращение идет не к одной функции. При твоем подходе вряд ли что то поможет, нужно коренным образом перестроить программу. 0

С++ для начинающих

Вычисление интеграла - C++