from numpy import *

x = array([[-0.7],[-2.5]])
for k in range(12):
    f = array([[x[0,0]**2*x[1,0] + 1],[exp(x[0,0]) + x[1,0]]])
    fder = array([[ 2*x[0,0]*x[1,0] , x[0,0]**2],[exp(x[0,0]) , 1]])
    x = x - linalg.inv(fder)@f
    print(f)
    print(x)