File: newpoly.py

package info (click to toggle)
pycode-browser 20120614%2Bgit%2Bb041dd2-7
  • links: PTS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 6,236 kB
  • ctags: 1,194
  • sloc: python: 2,220; xml: 152; makefile: 39
file content (17 lines) | stat: -rwxr-xr-x 243 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
from copy import copy

def coef(x,y):
	a = copy(y)
	m = len(x)
	for k in range(1,m):
		tmp = copy(a)
		for i in range(k,m):
			tmp[i] = (a[i] - a[i-1])/(x[i]-x[i-k])
		a = copy(tmp)
	return a

x  = [0,1,2,3]
y  = [0,3,14,39]
print coef(x,y)