File: newraph_plot.py

package info (click to toggle)
pycode-browser 1%3A1.03-1
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 9,860 kB
  • sloc: python: 2,919; xml: 152; makefile: 76
file content (25 lines) | stat: -rwxr-xr-x 310 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
 
#Newton Raphson method

from pylab import *

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

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

vf = vectorize(f)
x = linspace(-2, 5, 100)
y = vf(x)

x1 = 4
tg1 = df(x1)*(x-x1) + f(x1)
x1 = -3
tg2 = df(x1)*(x-x1) + f(x1)

grid(True)
plot(x,y)
plot(x,tg1)
plot(x,tg2)
ylim([-20,40])
show()