אנליזה נומרית/שיטת ניוטון-רפסון/קוד מקור

#!/bin/env python
from math import *
 
eps = 1e-3
 
def NR(f, f_tag, x0):
	x = float(x0)
	while abs(f(x))>eps:
		x = x - f(x)/f_tag(x)
		print "f(" + str(x) + ") = " + str(f(x))
 
 
f = lambda x: sin(2*x)
f_tag = lambda x: 2*cos(2*x)
 
NR(f,f_tag,pi/6)