# coding=utf-8
from numpy import *

def newtonfler(start,F,J):
  # Funksjonsprosedyre som utf?rer Newtons metode p? en 
  # vektorvaluert funksjon F med Jacobi-matrise J
  # Innparametre:
  # 
  # start: s?ylevektor som inneholder startpunktet for iterasjonen 
  # F: funksjon som regner ut funksjonsverdiene
  # J: Funksjon som regner ut Jacobideterminanten
  maxit=30
  tol=0.0000001
  diff=tol+1
  x=start
  for t in range(maxit):
    xny=x-linalg.solve(J(x),F(x))
    diff=linalg.norm(xny-x)
    x=xny
    print 'itnr=', t, 'x=', x, 'F(x)=', F(x) 
    if diff<tol:
      break
  return x