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

def f(x):
	return 2.0 * x**2 - 3*x - 5

def df(x):
	return 4.0 * x - 3

def nr(x, tol = 1.0e-9):
	for i in range(30):
		dx = -f(x)/df(x)
		#print x
		x = x + dx
		if abs(dx) < tol: 
			return x

print nr(4)
print nr(1)
print nr(-4)