clear all; fun = @(x) 3*x(1)^2 + 2*x(1)*x(2) + x(2)^2 - 4*x(1) + 5*x(2) warning('off'); x0 = [1,1]; [x,fval, eflag, out] = fminunc(@(x) fun(x),x0) x0 = [1,1]; [x,fval, eflag] = fminsearch(@(x) fun(x),x0) clear all; x0 = [0.5,0]; A = [1,2]; b = 1; Aeq = [2,1]; beq = 1; [x, fval] = fmincon(@rosenbrock,x0,A,b,Aeq,beq) FunObj = @(x) 1+x(1)/(1+x(2)) - 3*x(1)*x(2) + x(2)*(1+x(1)); lb = [0,0]; ub = [10,10]; x0 = [0.5,1]; x = fmincon(FunObj,x0,[],[],[],[],lb,ub)