function f=cost(x) f=1.10471*x(1)^2*x(2)+0.04811*x(3)*x(4)*(14.0+x(2)); subject to

g(1)=tau-tau_max; g(2)=sigma-sigma_max; g(3)=x(1)-x(4); g(4)=0.10471*x(1)^2+0.04811*x(3)*x(4)*(14+x(2))-5.0; g(5)=0.125-x(1); g(6)=delta-delta_max; g(7)=P-PcX; all g[1..7] should be <= 0

I have to keep any two of them as a constant and plot the f(x) The upper and lower bounds are respectively (x1..x4)