# coding=utf-8
from numpy import *

def newton(start,f,df):
  # Funksjon som utf?rer Newtons metode p? en funksjon 
  # f med derivert df spesifisert via de anonyme funksjonene f og df
  # Innparametre: 
  # start: startpunkt for iterasjonen 
  # f: anonym funksjon som returnerer funksjonsverdiene
  # df: anonym funksjon som returnerer verdiene til den deriverte
  maxit=30
  tol=0.0000001
  diff=tol+1
  x=start
  for t in range(maxit): 
    xny=x-f(x)/df(x)
    diff=abs(xny-x)
    x=xny
    print 'itnr=', t, ' x=', x, ' f(x)=', f(x)
    if diff<tol:
      break
  return x