Derivata och beräkningar i Python

Jag använder python 2.7.2

Jag är välldigt osäker på vad dom menar med uppgifterna 5a och 5b ock kan inte fråga min lärare fören tisdag.

Uppgiften är 5a och 5b från http://www.sm.luth.se/csee/courses/d0009e/labbar/lab2.html

och här är det jag har codat
problemet är att den retunerar 0 när den ska retunera ca 0.9911

import math
def derivative(f, x, h):

    t=float((1/(2*h))*(f(x+h)-f(x-h)))
    print t
    raw_input()
    return t

def solve(f, x0, h):
    '''
    x0 är startvärdet
    f är den funktion vars nollställe söks
    h är nogranheten
    '''
    f=float(f)
    x0=float(x0)
    t=0
    z=float(h/h/(10**h))
    while   round(t, h)!= round(f, h):
        #t = (f-((f**2-x0)/(2*f)))  #Det här är formeln jag får att fungera utan att anropa derivative() om jag har förstått uppgiften rätt
        t = (f-((f**2-x0)/(derivative(math.sin,x0,h))))
        print t
        raw_input()
        if round(t, h)!= round(f, h):
            
            #f = (t-((t**2-x0)/(2*t)))       #Det här är formeln jag får att fungera utan att anropa derivative() om jag har förstått uppgiften rätt
            f = (t-((t**2-x0)/(derivative(math.sin,x0,h))))
            print f
            raw_input()

    print round(t,h)

Fråga görna om jag har förklarat dåligt eller om ni undrar hur jag tengt.

det är derivative(math.sin,x0,h) som returnar 0 när jag skickar solve(12, 1395, 5)
dock om jag använder funktionen derivative(math.sin,1395,0.00001) så får jag tillbaka 0.9911868435075187

Anledningen till att round inte kraschar är för att round(t, h) har ett värde som första gången den testar är 0.00000 och round(f, h) har ett värde som första gången den testar är 1395.00000

-EDIT-

Nu har jag hittat det felet

z=float(h/h/(10**h))
retunerar 0.0

p=float(h)
z=(p/p/(10**p))
retunerar 0.00001

dock så får jag ett annat fel

Traceback (most recent call last):
File "<pyshell#1>", line 1, in <module>
solve(12,1395,5)
File "G:\Dropbox\Pyton 2.7\Testprogram\Loop.py", line 86, in solve
f = (t-((t**2-x0)/(derivative(math.sin,x0,z))))
OverflowError: (34, 'Result too large')

där line 86 är f = (t-((t**2-x0)/(derivative(math.sin,x0,z))))

Jag kollade på http://sv.wikipedia.org/wiki/Newtons_metod för att försöka förstå formeln och det är därifrån jag tagit 12 och 1395.
Som sagt så förstår jag inte rigtigt uppgiften så jag har fått göra som jag tror att dom menar. Om någon kan läsa uppgiften 5a och 5b och förklara i lite andra ord så att jag förstår så skulle det kanske gå bättre.

Sen om jag ska skicka in en formel tex x**2=1395 så får jag felmedelandet att x inte är definerat.