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

מתוך ויקיספר, אוסף הספרים והמדריכים החופשי
קפיצה לניווט קפיצה לחיפוש

שימו לב כי לצורך פתרון נומרי, עליכם לדעת לא רק את הפונקציה אלא גם את נגזרתה.

שפת פייתון[עריכה]

#!/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)