function [j] = fun( x)
% Objective function
% 
%Parameters
global p;
if p.case_I==1
    j=x(p.NT+2);
else
    j=-x(p.NT+8);
end
