Discussion :

[le_voisin]

salut,

je suis debutant en c++, et je veux faire un fonction me calculant une integrale par la methode des trapezes...
J'ai donc defini une fonction dont un des parametres est une fonction au sens mathématique.
Dans ma fonction (au sens c++) j ai une boucle dans laquelle ma focntion mathématique doit prendre une nouvelle valeurà chaque increment.
Voilà donc comment j ai programme la chose...:

ma fonction : 4 parametres : la fonction mathematique, la borne inf, la borne sup et le nombre de segments avec lequel on discretise le domaine.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
float Integration_Trapeze(float f, float a, float b, float n)
	{
		float I=0; /*valeur de l'intégrale*/
		float h=(b-a)/n; /*pas d'intégration*/
		float x=a;
		float fa=0; //initialisation
		float fb=0; //initialisation
 
		for (int i=0; i<n; i=i+1)
		{
			fa=f;
			x=x+h; fb=f;
			I=0.5*h*(fa+fb);
		};
		return(I);
}

Ma classe main :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
#include <iostream> 
#include <cmath> 
using namespace std;
 
float Integration_Trapeze(float, float, float, float);
 
int main(void)
{
float x=0; //initialisation
float f=x*x; //definition de la fonction
float I=Integration_Trapeze(f, 0, 2, 10);
cout<<"integrale = "<<I<<"\n";
}

Voilà, si quelqu un peut me filer un coup de main ou un endroit ou je pourrais trouver l info ca serait sympa...
Tant que j y suis autre question, comment appeler le nombre Pi en c++ ?

Merci

[le_voisin]

oui j ai un peu oublie d expliciter clairement le problème : c est tout simplement que ma fonction ne change pas de valeur dans la boucle !

[khayyam90]

Bien le bonjour,

le premier argument de Integration_Trapeze est un float. Lors de l'appel, ce float vaut 0. Oui, il vaut 0*0.
Ce que tu passes à Integration_Trapeze est donc la valeur 0. Ce n'est pas une fonction au sens mathématique (qui soit dit en passant n'existe pas en C++).

Ensuite dans la boucle dans Integration_Trapeze tu fais toutes les opérations que tu veux, mais tu as toujours f qui vaut 0.

Et pour pi, un simple #define PI 3.14159 peut suffire. Sinon tu as la constante flottante M_PI qui vaut 3.14159 déclarée dans <cmath>

[bigboomshakala]

pas mieux, tu auras toujours 0 comme ça.
sinon c'est plus de l'algorithmie qu'un problème de prog je pense (cf. forum algorithmes)

[Laurent Gomila]

Une solution à ton problème serait de passer en paramètre une fonction (au sens C++), plutôt qu'un flottant.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
// Définition du type des fonctions que tu vas utiliser
typedef float (*TypeFonction)(float);
 
float Integration_Trapeze(TypeFonction f, float a, float b, float n)
	{
		float I=0; /*valeur de l'intégrale*/
		float h=(b-a)/n; /*pas d'intégration*/
		float x=a;
		float fa=0; //initialisation
		float fb=0; //initialisation
 
		for (int i=0; i<n; i=i+1)
		{
			fa=f(x);
			x=x+h; fb=f(x);
			I=0.5*h*(fa+fb);
		};
		return(I);
}
 
float Carre(float x)
{
    return x * x;
}
 
int main(void)
{
    float I=Integration_Trapeze(&Carre, 0, 2, 10);
    cout<<"integrale = "<<I<<"\n";
}

[le_voisin]

ok je vous remercie la dernière reponse me semble interessante mais je n'ai pas tout compris !
Est ce que quelqu'un pourrait m expliquer ces lignes de commandes ou du moins me dire dans quelle rubrique je pourrais trouver des infos me permettant de les utiliser correctement:

// Définition du type des fonctions que tu vas utiliser
typedef float (*TypeFonction)(float);
et

int main(void)
{
float I=Integration_Trapeze( &Carre, 0, 2, 10);
cout<<"integrale = "<<I<<"\n";
}