黄金分割法最优化python

python练习

import math
def f(x):
    return x**3-3*x+1
def golds(a,b,delta,epsilon):
    t = (math.sqrt(5)-1)/2
    h = b - a
    fb = f(b)
    fa = f(a)
    p = a+(1-t)*h
    q = a+t*h
    fp = f(p)
    fq = f(q)
    g = []
    k = 1
    g.append([a, p, q, b])
    while(abs(fb-fa)>epsilon or (h>delta)):
        if(fp