import optimtool as oo
from optimtool.base import np, sp, plt
pip install optimtool # any version with accessible L_BFGS
加载L_BFGS方法
import optimtool.unconstrain as ou
lbfgs = ou.newton_quasi.L_BFGS
f ( x ) = ∑ i = 1 n − 1 ( x i + x i + 1 − 3 ) 2 + ( x i − x i + 1 + 1 ) 4 , x 0 = [ 2 , 2 , . . . , 2 ] . f(x)=\sum_{i=1}^{n-1}(x_i+x_{i+1}-3)^2+(x_i-x_{i+