Objective function

%===========================================%
%Power Electronics in Distribution Networks %
%Coursework Task 2 Objective function       %
%===========================================%
function f=Psop(x)
global V1;
global G12;
global B12;
global G13;
global B13;
global G14;
global B14;
global G22;
global B22;
global G23;
global B23;
global G33;
global B33;
global G24;
global B24;
global G34;
global B34;
global G44;
global B44;
global PG;
global PSOP;
%Psop=PG-P2
%When calculate maximum value of Psop, change signe of f
f=(PG-(V1*x(1)*(G12*cos(x(4))+B12*sin(x(4)))+x(1)^2*G22+x(1)*x(2)*(G23*cos(x(4)-x(5))+B23*sin(x(4)-x(5)))+x(3)*x(1)*(G24*cos(x(4)-x(6))+B24*sin(x(4)-x(6)))));
end