function F = CostFcn1(x,C,ub)
F = f_t(x,C,ub) + OM_DG(C) + TCPD_BES(x,C);
function ft = f_t(x,C,ub)
N = cell(1,2);
ft = x(2)*C.B_MT+ub;
for i = 1:2
N{i} = C.B_MT*x(i+2)+ub;
end
ft = ft(x) + sum(arrayfun(@(N_q) N_q{1}(x), N));
end
function om_dg = OM_DG(C)
om_dg = C.T*C.OM_MT;
end
function tcpd = TCPD_BES(x,C)
tcpd = x(1)*C.P_WTmax;
end
end
function F = CostFcn1(x,C,ub)
om_dg = C.T*C.OM_MT;
tcpd = x(1)*C.P_WTmax;
N = zeros(1,2);
ft = x(2)*C.B_MT+ub;
for i = 1:2
N(i) = C.B_MT*x(i+2)+ub;
end
ft = ft + sum(N);
F = ft + om_dg + tcpd;
end