Forum de mathématiques - Bibm@th.net
Vous n'êtes pas identifié(e).
- Contributions : Récentes | Sans réponse
#1 30-05-2006 16:25:17
- SANDY
- Membre
- Inscription : 30-05-2006
- Messages : 2
[Maple] Faire une procédure de méthode de trapèze
Bonjour j' ai voulu créer une procéture maple qui calcul une intégral à partir de la méthode des trapèzes
J'ai des erreurs pourrer vous m'aider ?
merci
code:
Trap:=proc(a,b,n)
> local resultat,intervalle;
> intervalle:=(b-a)/n;
> resultat:=0;
> initialisation:=(f(a)+f(b))/2;
> for i from 1 to n do
>
> resultat=f(((a+(n-1)*(b-a))/n));
>
> end do;
> resulat:=resultat*(intervalle+initialisation);
> print(evalf(resultat));
> end proc;
Hors ligne
#2 31-05-2006 07:37:04
- MasterJ
- Membre
- Inscription : 09-05-2006
- Messages : 10
Re : [Maple] Faire une procédure de méthode de trapèze
A quoi te sert ta boucle for ? Tu n'utilises même pas le i de la boucle ?
A chaque itération tu écrases la valeur de réusltat par une autre... au lieu d'en faire la somme... et au final je comprends plus ce que tu calcules....
Si tu vas voir : http://www.bibmath.net/dico/index.php3? … ntnum.html
tu t'apperçois que l'intégrale entre a et b est la SOMME des aires de tous les trapèzes... donc à chaque itération de ta boucle, tu rajoutes l'aire du trapèzes suivant, et au final tu as ton intégrale ^^
Je suis pas un pro de Mapple (j"utilise plutôt Matlab ^^) mais je ferais plutôt comme ça :
Trap:=proc(a,b,n)
> local resultat,intervalle;
> intervalle: = (b-a)/n;
> resultat:=0;
> for i from 1 to n do
>
> resultat:=resultat + intervalle*( f(a+(i-1)*intervalle) + f(a + i*intervalle) ) ;
>
> end do;
> print(evalf(resultat));
> end proc;
^^
Dernière modification par MasterJ (31-05-2006 07:39:19)
Hors ligne
#3 09-06-2013 21:49:24
- thaddée
- Invité
Re : [Maple] Faire une procédure de méthode de trapèze
ce qui serait super dans cette procédure c'est aussi de tracer le graphe de la fonction + les trapezes quelqu'un aurait une idée?
#4 09-06-2013 21:56:45
- thaddée
- Invité
Re : [Maple] Faire une procédure de méthode de trapèze
si ça peut aider quelqu'un... http://alamanya.free.fr/cours/chapn71.php#3