BatterySizeRange=(0:1000:100*10^3);
HouseNumberRange=5;
p=randi(96, HouseNumberRange, 1);    


for i = 1:length(BatterySizeRange)
    [Cost_Benefit_d(i),PAR_d(i)] = FeederDomestic(BatterySizeRange(i),HouseNumberRange,p);
end

plotyy([BatterySizeRange,BatterySizeRange] , [PAR,PAR_d],[BatterySizeRange,BatterySizeRange] , [Cost_Benefit.',Cost_Benefit_d])