codepad
[
create a new paste
]
login
|
about
Language:
C
C++
D
Haskell
Lua
OCaml
PHP
Perl
Plain Text
Python
Ruby
Scheme
Tcl
def computeRoot(poly, x_0, epsilon): ctr=-1 y=0 dy=42 epsilon=abs(epsilon) while (y>epsilon or y<-epsilon) or ctr==-1: ctr=ctr+1 print ctr x_0=x_0-(y/dy) print "x_0=",x_0 y=0 dy=0 for thing in range(len(poly)): y+=poly[thing]*(x_0**thing) dy+=(thing)*poly[thing]*(x_0**(thing-1)) print 'y=',y return [x_0,ctr] print computeRoot([1, -1, 1, -1], 2, .001)
Private
[
?
]
Run code
Submit