#!/usr/bin/env python
# -*- coding: iso-8859-1 -*-
# interp - Programa para interpolar datos usando
# el polinomio de Lagrange cuadratico para tres puntos dados.
# Importa un modulo que permite graficar.

from pylab import *

def fun(x):
    return x*x

def dfun(x):
    return 2*x

def err(x,h):
    return abs(dfun(x)-(fun(x+h)-fun(x))/h)

def main():
    x_min=-10
    x_max=1

    xo=1
    
    ho=10.0
    x1=[]
    x2=[]
    y1=[]
    for i in range(20):
        ho=ho/10.0;
        x1.append(ho)
        x2.append(-i)
        y1.append(log10(err(xo,ho)))
        

    print x1
    print x2
    print y1

    plot(x2,y1,'ro',label='line 1', linewidth=2)
    plot(x2,y1)
    title("ERROR en derivada")
    axis([-20, 0, -10, 1])
    xlabel(r'$\log\; h$')
    ylabel(r'$\log\; \Delta(h)$')
    grid(True)
    show()

    ##plot([1,2,3], [1,2,3],  label='line 1', linewidth=2)
##    plot([1,2,3], [1,4,9], 'rs',  label='line 2')
##    axis([0, 4, 0, 10])
##    legend()
##    title("Polinomio de interpolacion de Lagrange")
##    grid(True)
##    show()

if __name__=='__main__':
    main()