#!/usr/bin/env python
# coding: utf-8

# In[32]:


from sympy import *


# In[33]:


def calcNewton(f, df, xn):
    return [N(f.subs(x, xn)), N(df.subs(x, xn)), N(xn - (f.subs(x, xn)/df.subs(x, xn)))]


# In[34]:


x0 = 1.0
x = symbols("x")
f = cos(x) - x
iteraciones = 4
df = f.diff(x)
init_printing()

display(f)
print("n \t xn \t\t f(xn) \t\t f'(xn) \t xn+1")
xn = x0
for n in range(iteraciones):
    xn1 = calcNewton(f, df, xn)
    print("%d \t %0.9f \t %0.9f \t %0.9f \t %0.9f" %(n, xn, xn1[0], xn1[1], xn1[2]))
    xn = xn1[2]