function [Min_fx] = yourquestion(Bgi,Pgi,MP,Pgrid)
NT = 24; %no. of hours ie.24
Ng = 5; %no. of DG which is 5
evec = [];
for t = [1:1:NT]
for i = [1:1:Ng]
B = [];
A = (Bgi(i,1).*(Pgi(i,1).^t)) + ((MP(i,1).^(t))*(Pgrid(i,1).^(t)));
A = A';
evec = [evec A];
B = evec;
end
evec = [evec B];
end
Min_fx = sum(evec,'all');
end
