from math import exp # Risolve l'equazione f(x)=0 # definizione della funzione f def f(x): return exp(x)+x # derivata prima def f1(x): return exp(x)+1 a=-1.0 # stima per difetto b=0 # stima per eccesso xn = (a+b)/2 # punto di partenza per il metodo di newton (x1,x2)=(a,b) # punto di partenza per il metodo delle secanti p = (a+b)/2 # punto di partenza per il metodo del punto fisso eps = 0.001 # precisione while abs(x2-x1) > eps: # il calcolo della precisione viene fatto sul metodo delle secanti, se lo facessi sul metodo di bisezione convergerebbe troppo lentamente e avrei un errore di divisione per zero # bisezione m = (a+b)/2 if f(m)*a > 0: # se f(m) e a sono concordi, sposto a a destra a = m else: b = m # Newton xn = xn - f(xn)/f1(xn) # Secanti x3 = x2 - (x2 - x1)/(f(x2) - f(x1))*f(x2) (x1,x2) = (x2,x3) # Punto fisso c = -1/f1(p) p = c*f(p)+x # Stampa della tabella print "{0:.10f 1:.10f 2:.10f 3:.10f}".format(m,xn,x3,p)
Run
Reset
Share
Import
Link
Embed
Language▼
English
中文
Python Fiddle
Python Cloud IDE
Follow @python_fiddle
Browser Version Not Supported
Due to Python Fiddle's reliance on advanced JavaScript techniques, older browsers might have problems running it correctly. Please download the latest version of your favourite browser.
Chrome 10+
Firefox 4+
Safari 5+
IE 10+
Let me try anyway!
url:
Go
Python Snippet
Stackoverflow Question