#!/usr/bin/python3
import numpy as np
import math
import matplotlib.pyplot as plt
import mpl_toolkits.mplot3d as Axes3D
import mysearch as mys

def f(x,y):
    return ((x-1)**2+y**2)


xn = np.array([-2,-1])
pn = np.array([1,0])

def fline(t):
    x=xn+t*pn
    return f(x[0],x[1])
opt = []
blist = mys.minbracket(fline,0,4,25)
for b in blist:
    tmin = mys.golden(fline,b[0],b[1],1e-5)
    print(xn+tmin*pn)