Discussion :

[medphy]

Bonjours a tous:

Je débute en python et j'ai quelque petit probleme.
J'ai l'intégrale suivante : P(x)=Intégrale (4x**6+4**4+9)**0.5
je voudrais savoir comment faire cette intégrale de 0 a X, avec X qui varie de 0 a 6 et stocké toute ces valeurs dans un tableau (avec numpy). (avec un pas petit)

Merci de votre aide, par avance.

[Matthieu Brucher]

Tu calcules dans un tableau l'intégrale partielle, et tu fais un cumsum dessus.

Code :

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

1
2
t = valeurs[0:-1] + valeurs[1:]
i = numpy.cumsum(t) / dt

ça, c'est pour une intégration par trapèze.

[vikki]

Bonjours a tous:

Je débute en python et j'ai quelque petit probleme.
J'ai l'intégrale suivante : P(x)=Intégrale (4x**6+4**4+9)**0.5
je voudrais savoir comment faire cette intégrale de 0 a X, avec X qui varie de 0 a 6 et stocké toute ces valeurs dans un tableau (avec numpy). (avec un pas petit)

Merci de votre aide, par avance.

Hello,
Tu cherche une méthode existante dans les bibliothèques Python ou à coder un algo de méthode numérique?
Dans le premier cas, regarde du coté de scipy, ca existe peut être (à voir).
Sinon un tit coup de Wikipedia sur le sujet va grandement t'aider

[medphy]

voila je crois que je tient un debut de reponse:
#Dans ce petit programme on calcule p(x,y) qui correspond a la longueur de la courbe en fonction de x et y.
from pylab import *
import math as M
import numpy as N

#Correspond a mon pas d'integration
dt=0.01

# Correspond aux valeurs prise par t, valeur du parametrage
T=N.arange(-3,3,dt)

#ns = nbre de case du tableau
ns= 6/0.01
P=N.zeros(ns)

#dans le tableau remplie de 0 de p on fait une boucle pour y metre toute nos valeurs de l'integrale
for i in range(0, ns-1):

P[i+1]=P[i]+ ( ( (1/3)*(4*(T[i]**6))+49*(T[i]**4)+9-28*(T[i]**4) )**0.5 )*dt


print P

[Matthieu Brucher]

Ta boucle for est à éviter. Avec Numpy (comme avec MAtlab), il faut penser vectoriel/matriciel.

Tu devrais commencer par calculer les valeurs temporaires basées sur les T et ensuite, utilise cumsum comme indiqué dans mon message précédent.

[vikki]

Si le but est simplement d'avoir le résultat, regarde à http://docs.scipy.org/doc/scipy-0.7....integrate.html, y'a des exemples d'utilisation assez claires.