function f(x)
z=(x[1]-2)^2+(x[2]-1)^2+0.04/(-0.25x[1]^2-x[2]^2+1)+5(x[1]-2*x[2]+1)^2
return z
end
using TyOptimization
x0=[0,0]
opt1 = optimoptions(:fminunc, Algorithm="quasi-newton",HessianApproximation="bfgs")
x1,=fminunc(f,x0,opt1)
opt2 = optimoptions(:fminunc,Algorithm="quasi-newton",HessianApproximation="dfp")
x2,=fminunc(f,x0,opt2)
opt3 = optimoptions(:fminunc,Algorithm="trust-region",SpecifyObjectiveGradient="false")
x3,=fminunc(f,x0,opt3)
println(x1)
println(x2)
println(x3)