wolfe 算法_最优化课课后作业笔记

from numpy import *
import sys
MAXINT=sys.maxint
def gradient_(X):
    x1=X[0];x2=X[1]
    return array([(-400*(x2-x1**2)*x1-2*(1-x1)),200*(x2-x1**2)])
def f(X):
    x1=X[0];x2=X[1]
    r