File: taylor.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 (15 lines) | stat: -rwxr-xr-x 280 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
from numpy import *

p = poly1d([1,1,1,0])
dp = p.deriv()
dp2 = dp.deriv()
dp3 = dp2.deriv()

a = 0  # The known point
x = 0
while x < .5:
	tay = p(a) + (x-a)* dp(a) + \
         (x-a)**2 * dp2(a) / 2 + (x-a)**3 * dp3(a)/6
	print '%5.1f  %8.5f\t%8.5f'%(x, p(x), tay)
	x = x + .1