from math import *
import numpy  as np

# Function definition 
def g(x):
  #return (1/2)*(x**2)- sin(x)
  return x**3 - 12.2*x**2 + 7.45*x + 42

# Implementation of Secant method

def secant(x0,x1,eps):
  
  i=0

  x2 = ((g(x1)*x0) - (g(x0)*x1))/(g(x1)-g(x0))

  x2 = np.float("{:.2f}".format(x2))

  arr = np.empty((5), dtype=object)

  arr1 = np.empty((5), dtype=object)

  arr[0] = i; arr[1] = x0; arr[2] = x1; arr[3] = x2; arr[4] = g(x2)

  while abs(x2-x1) > abs(x1)*eps:

    i = i+1 

    x0 = x1

    x1 = x2

    x2 = ((g(x1)*x0) - (g(x0)*x1))/(g(x1)-g(x0))

    x2 = np.float("{:.2f}".format(x2))

    arr1[0] = i;arr1[1] = x0; arr1[2] = x1; arr1[3] = x2; arr1[4] = g(x2)
    
    arr = np.vstack((arr,arr1))
  
  arr1 = np.array(["Iter-k","x0","x1","x2","g(x)"])
   
  arr = np.vstack((arr1,arr))

  for i in arr:

      print(i)


x0 = 13
x1 = 12
eps = 0.00001
secant(x0,x1,eps)

# Function definition 
def g1(x):
  
  return ((2*x-1)**2) + (4*(4-1024*x)**4)

# Implementation of Secant method

def secant(x0,x1,eps):
  
  i=0

  x2 = ((g1(x1)*x0) - (g1(x0)*x1))/(g1(x1)-g1(x0))

  arr = np.empty((5), dtype=object)

  arr1 = np.empty((5), dtype=object)

  arr[0] = i; arr[1] = x0; arr[2] = x1; arr[3] = x2; arr[4] = g1(x2)

  while abs(x2-x1) > abs(x1)*eps:

    i = i+1 

    x0 = x1

    x1 = x2

    x2 = ((g1(x1)*x0) - (g1(x0)*x1))/(g1(x1)-g1(x0))

    arr1[0] = i;arr1[1] = x0; arr1[2] = x1; arr1[3] = x2; arr1[4] = g1(x2)
    
    arr = np.vstack((arr,arr1))
  
  arr1 = np.array(["Iter-k","x0","x1","x2","g(x)"])
   
  arr = np.vstack((arr1,arr))

  for i in arr:

      print(*i)


x0 = 0
x1 = 1
eps = 0.00001
secant(x0,x1,eps)