plot1.py

# coding=latin-1
import numpy as np
from pylab import *

f1 = lambda x: x+5.-x**2
f2 = lambda x: 5./x
f3 = lambda x: 0.5*(x+5/x)
y  = lambda x: 2*sqrt(5) -x
x=linspace(1.5,3.5,100)
plot(x,x, 'r--')
plot(x,f1(x))
plot(x,f2(x))
plot(x,f3(x))
plot(x,y(x), 'r--')
title(r'Métodos iterativos para resolver $x^2-5=0$')
grid()
show()