extrapolation avec python
Hello !
J'essaye d'écrire un code python pour faire l'extrapolation à partir de n'importe quelle courbe (pas seulement linéaire), sans utilisation de numpy et scipy.
voici ce que j'ai écrit :
soit x la liste des valeurs de l'axe x, et y la liste des valeurs de l'axe y :
Code:
1
2
3
4
5
6
7
8
9
10
11
12
13
| k=0
if P>max(x) or P<min(x) :
for i in range (0,len(y)):
s=1
t=1
for j in range (0,len(y)) :
if (j != i) :
s=s*(y[i]-y[j])/(i-j)
t=t*(P-y[i])
k=k+y[i]+t*s
print k
Y_extrapolated=str(k)
print 'Y_extrapolated is', Y_extrapolated |
sauf que ce code ne me donne pas des valeurs acceptables.
Quelqu'un aurait-il des suggestions svp?
Merci
