function [vector] = FeederDomestic(DomesticBatterySize,SolarHouses,p)

 load('FeederSolar.mat')
 load('AllData.mat')
 %load('C:\Users\CDT6\Documents\Summer Project\MATLAB\TVV\Feeder_date_time.mat')

n=1;



AllData=cell2mat(AllData);
AllData=AllData(:,[40:94]);
AllDataMat=AllData;

while n<=SolarHouses
    House=AllDataMat(:,p(n));
    Grid=zeros(length(FeederSolar),1);
    B=zeros(length(FeederSolar),1);
    for i=1:length(FeederSolar)
         B(i)= ((FeederSolar(i)*1.14)-House(i))*0.8;
         B(1)=0.5*(DomesticBatterySize*2/SolarHouses);
         if  (0.2*(DomesticBatterySize*2/SolarHouses) <= sum(B))&&(sum(B) < 0.8*(DomesticBatterySize*2/SolarHouses) )
            Grid(i)=0;
         elseif sum(B) >= (DomesticBatterySize*2/SolarHouses)*0.8
            B(i)=0;
            Grid(i)= (FeederSolar(i)*1.14)-House(i);
         elseif sum(B)< (DomesticBatterySize*2/SolarHouses)*0.2
            B(i)=0;
            Grid(i)= (FeederSolar(i)*1.14)-House(i);
     end
     AllDataMat(:,p(n))=-Grid;
end
    n=n+1;
end

vector=sum(AllDataMat,2);

% Peak_Mean_Ratio =zeros(length(FeederSolar),1);
% GridS=sum(AllDataMat,2);
% 
% for i= 1:length(FeederSolar)
%     if i<49
%         Peak_Mean_Ratio(i)=0;
%         Peak_MRB(i)=0;
%     else
%        Peak_Mean_Ratio(i)=max(GridS(i-48:i))/mean(GridS(i-48:i));
%        Peak_MRB(i)=max(AllDataSol(i-48:i))/mean(AllDataSol(i-48:i));
%     end 
% end
% 
% PMR_Before=mean(Peak_MRB);
% 
% PMR=mean(Peak_Mean_Ratio);
% 
% Cost_Benefit=(1/(DomesticBatterySize*2*SolarHouses))*(-PMR+PMR_Before);