Discussion :

[rambc]

Bpnjour,
la question a d'être posé plein de fois (donc je m'excuse de le faire une fois de plus).

Etant donné une fonction d'une seule variable f(x), peut-on avoir un tableau de valeurs du type [f(a), f(a+2h), f(a+3h), ... , f(a+nh)] où a et h sont connus ?

Peut-on avoir une fonctionnalité équivalente pour une fonction de deux variables ?

[Matthieu Brucher]

Je pense que tu cherches la fonction vectorize pour faire f([a, a+h, a+2h, ...]), sachant que les arguments sont créés par :

Code :

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

args = a + h*arange(0, n)

[tyrtamos]

Bonjour,

Voilà une solution possible en utilisant une "list comprehension":

Code :

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

1
2
3
4
5
6
7
8
9
 
def f(x):
    return 3*x
 
a = 5
T = [f(a+h) for h in range(0,100,7)]
print T
 
[15, 36, 57, 78, 99, 120, 141, 162, 183, 204, 225, 246, 267, 288, 309]

Ça me semble assez facile à adapter à 2 variables, mais ça dépend ce que tu veux faire exactement. Si tu ne trouves pas, donne un exemple.

Tyrtamos

[rambc]

En fait je cherche un prog. à qui je donnerais la fonction et qui se chargerait ensuite des calculs. Je pensais à NumPy ou symPy. Quelqu'un connait-il la procédure ?

[Matthieu Brucher]

Pardon, je n'ai pas précisé, mais vectorize et arange sont des fonctions numpy

[rambc]

Ok Matthieu Brucher. Au temps pour moi...

Merci.

[rambc]

Ceci pourra être utile à d'autres.

Code :

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

1
2
3
4
5
6
7
8
9
# -*- coding: iso-8859-15 -*-
#!/usr/bin/env python
 
from numpy import *
 
x = arange(0,10,1)
y = pow(x,2)
 
print y

Ceci donnera les carrés de 0 à 10 en augmentant les valeurs de x de un en un :

Code :

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

[ 0  1  4  9 16 25 36 49 64 81]

[Matthieu Brucher]

Le from numpy import all, faut éviter

[rambc]

Salut,
tu dis cela à cause des risques de conflit de nom de fonctions.

[Matthieu Brucher]

Exactement, même si Numpy change les noms de fonctions pour éviter cela.

[rambc]

OK. Je vais à l'avenir changer mes importations avec des import ... as ... .