Вывести результат (Вычисление интегралов)

[prikolist]

Здравствуйте программисты! Подскажите пожалуйста,как вывести результат в этом коде, когда-то брался решать,но не дорешал,подскажите пож.как вывести результат. Вот код:

#include <vcl>
#include <iostream>
#include <conio>
using namespace std;

double Simpson(double (*Fx)(double))
{
        double I;
        double h;
        double x;
        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;
 
        h=(b-a)/(m-1);
        x = a;
 
        I = Fx(a) + Fx(b);
        int n=0;
 
        while (n < m-2)
        {
                x = x+h;
                if (n%2==0) I=I+4*Fx(x);
                        else  I = I+2*Fx(x);
                n++;
        }
 
        return I*h/3;
}

int main()
{
system("pause");
}

Мне нужно чтоб вычисляло результат 2 интегралов:
1 -ый интеграл какой вычислить нужно (1/sqrt(0.5*pow(x,2)+1.5))
2 - ой интеграл (x+1)*cos(x*x)

[Sazary]

Объявляете функции:

double F1(double x)
{
 return (1/sqrt(0.5*pow(x,2)+1.5));
}

double F2(double x)
{
return (x+1)*cos(x*x);
}

И в главной функции:

 cout<<Simpson(&F1)<<endl;
 cout<<Simpson(&F2)<<endl;

[prikolist]

Ясно,спасибо, попробывал вроде работает, осталось последнее,методом трапеций, не знаю как сделать. У вас Сазари где-то видел такой код, по трапециям:

h=(b-a)/n;
 
for(i=0,x=a; i<n; i++,x+=h)
 F[i] = 1/((0.5 + 0.1*pow(x,3))*sqrt(x));
 
I1 = 0;
for(i=1;i<n-1;i++)
 I1 = I1+F[i];
I1 = I1+ F[0]/2 + F[n-1]/2;
I1 =I1 * h;

Но не знаю, как сделать так,чтоб был похожее, как выше у Симпсона,только по трапециям.

[Sazary]

Что-то вроде этого:

for(i=1,x=a; i<m-1; i++,x+=h)
         I += Fx(x);
        I += Fx(a)/2 + Fx(b)/2;
        I *= h;
        return I;

[MaTBeu]

Вычисление интегралов:
1. Метод Левых прямоугольников

double M_leftRectangle( pointFunc F, double a, double b,unsigned N)
{
        double Res = 0.0;
        double h=(b-a)/N;
        for (double x=a; x<b;x+=h)
                Res+= (*F)(x);
        return Res*h;
}

2. Метод Правых прямоугольников

double M_rightRectangle( pointFunc F, double a, double b,unsigned N)
{
        double Res = 0.0;
        double h=(b-a)/N;
        for (double x=a+h; x<=b;x+=h)
                Res+= (*F)(x);
        return Res*h;
}

3. Метод средних прямоугольников

double M_middleRectangle( pointFunc F, double a, double b,unsigned N)
{
        double Res = 0.0;
        double h=(b-a)/N;
        for (double x=a+h/2; x<=b; x+=h)
                Res+= (*F)(x);
        return Res*h;
}

4. Метод Трапеций

double M_trapet( pointFunc F, double a, double b,unsigned N)
{
        double Res = 0.0;
        double h=(b-a)/N;
        for (double x=a+h; x<b; x+=h)
                Res+= (*F)(x);
        Res=h*((*F)(a)+(*F)(b))/2+Res*h;
        return Res;
}

5. Метод Симпсона

double M_Simpson( pointFunc F, double a, double b,unsigned N)
{
        double Res = 0.0;
        double h=(b-a)/N;
        for (double x=a+h; x<b; x+=2*h)
                Res+= 2*(*F)(x);
        for (double x=a+2*h; x<b; x+=2*h)
                Res+= 4*(*F)(x);
        Res=h*((*F)(a)+(*F)(b)+Res)/3;
        return Res;
}

[prikolist]

Попробывал так, что-то ерунда получилась.

#include <vcl>
#include <iostream>
#include <math>
#include <conio>
using namespace std;


double Trapesia(double (*Fx)(double))
{
        double I;
        double h;
        double x;
        double a,b;
        int m,i;


        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;

        for(i=1,x=a; i<m-1; i++,x+=h)
         I += Fx(x);
        I += Fx(a)/2 + Fx(b)/2;
        I *= h;
        return I;
 }

double F1(double x)
{
return (1/sqrt(0.5*pow(x,2)+1.5));
}

double F2(double x)
{
return (x+1)*cos(x*x);
}


int main()
{
cout<<Trapesia(&F1)<<endl;
cout<<Trapesia(&F2)<<endl;
cout<<"\n Integral = "<<I;
system("pause");
}

[Sazary]

Так а подсчет h зачем убрали?

h=(b-a)/m;

А вообще, вон вам MaTBeu какую подборку сделал.. Берите и используйте )

[MaTBeu]

Только под мои функции в начале программы нужно написать

typedef double (*pointFunc)(double);

так вы определите указатель на функцию.

[prikolist]

Попробывал так,согласно инструкции Матвея останавливается на первой строке функции мэин и дальше не идёт.

#include <vcl>
#include <iostream>
#include <math>
#include <conio>
using namespace std;
double F1(double x);
double F2(double x);
typedef double (*pointFunc)(double);


double M_Simpson( pointFunc F, double a, double b,unsigned N)
{
        double Res = 0.0;
        double h=(b-a)/N;
        for (double x=a+h; x<b; x+=2*h)
                Res+= 2*(*F)(x);
        for (double x=a+2*h; x<b; x+=2*h)
                Res+= 4*(*F)(x);
        Res=h*((*F)(a)+(*F)(b)+Res)/3;
        return Res;
}

double F1(double x)
{
return (1/sqrt(0.5*pow(x,2)+1.5));
}

double F2(double x)
{
return (x+1)*cos(x*x);
}


int main()
{
cout<<Trapesia(&F1)<<endl;
cout<<Trapesia(&F2)<<endl;

system("pause");
}

[Sazary]

Функция

double M_Simpson( pointFunc F, double a, double b,unsigned N)

А вызываете:

cout<<Trapesia(&F1)<<endl;
cout<<Trapesia(&F2)<<endl;

Ничего не замечаете?

И еще на количество параметров обратите внимание.