Cómo se hace?
Publicado por Carlos (11 intervenciones) el 06/06/2012 19:40:00
Alguien tiene algún código para poder sacar el polinomio de interpolacion de Lagrange? He encontrado este:
Pero no entiendo que significa tmp = scipy.poly1d([0])
result=scipy.poly1d([0])
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
def lagrange(x):
tmp = scipy.poly1d([0])
result=scipy.poly1d([0])
for i in x.keys():
numerator=scipy.poly1d([1])
denom = 1.0
for j in x.keys():
if (i != j):
tmp = scipy.poly1d([1,-j])
numerator = numerator * tmp
denom = denom * (i - j)
tmp = (numerator/denom) * x.get(i)
result = result + tmp
return result
Pero no entiendo que significa tmp = scipy.poly1d([0])
result=scipy.poly1d([0])
Valora esta pregunta


0