map214pf.py

from pylab import *

ion()

def coweb (x,g):
    plot([x,x],[x,g(x)],[x,g(x)],[g(x),g(x)], color='red')
    return g(x)
## caso concreto:
g = lambda x: (x**2+5)/(2*x)
x = linspace(1,3,100)
plot(x,g(x), x, x)
grid()
draw()
x0=1.5
while (abs(g(x0)-x0) >= 0.0001):
    x0=coweb(x0,g)
    draw()
    raw_input("tecle algo: ")
print "Acabou ..."