-- Oct 18 In-Class Exercise
Xn+1 = Xn - f(xn)/f'(xn)
f'(xn) = 2*(x-2)*(x-3)^2 + 2*(x-2)^2*(x-3)
x1, where x0=5, => 5- 3^2*2^2 / (60) = 4.4
xs = [5]
def f(x):
return ((x-2)**2)*(x-3)**2
def f_prime(x):
return (2*(x-2))*(x-3)**2 + 2*((x-2)**2)*(x-3)
print "x0", xs[0]
for i in range(0,4):
temp = xs[i] - f(xs[i])*1.0 /f_prime(xs[i])
print "x" + str(i+1), temp
xs.append(temp)
- x0 5
- x1 4.4
- x2 3.95789473684
- x3 3.63629108873
- x4 3.40722182808
(
Edited: 2017-10-18)
Xn+1 = Xn - f(xn)/f'(xn)
f'(xn) = 2*(x-2)*(x-3)^2 + 2*(x-2)^2*(x-3)
x1, where x0=5, => 5- 3^2*2^2 / (60) = 4.4
xs = [5]
def f(x):
return ((x-2)**2)*(x-3)**2
def f_prime(x):
return (2*(x-2))*(x-3)**2 + 2*((x-2)**2)*(x-3)
print "x0", xs[0]
for i in range(0,4):
temp = xs[i] - f(xs[i])*1.0 /f_prime(xs[i])
print "x" + str(i+1), temp
xs.append(temp)
*x0 5
*x1 4.4
*x2 3.95789473684
*x3 3.63629108873
*x4 3.40722182808