Leach_Protocol.m

function Result=Leach_Protocol(Model,Network)
Iteration=Model.Iteration;
Packet.ToCH=zeros(Iteration,1);
Packet.ToBS=zeros(Iteration,1);
Dead.All=zeros(Iteration,1);
Dead.Normal=zeros(Iteration,1);
Dead.Advance=zeros(Iteration,1);
HeadClusters=zeros(Iteration,1);
Networks=[];
%% Main loop
for i=1:Iteration
    disp(['It = ' , num2str(i)]);
    % netwrok snapshot
    if(mod(i-1,Model.SaveIteration)==0)
        Networks=[Networks,Network];
    end            
   
    if((mod(i-1, round(1/Model.P) )==0))
        Temp=num2cell(zeros(1,length(Network)));
        [Network.G]=Temp{:};
        [Network.cl]=Temp{:};                                                        
    end        
    [Dead.Normal(i),Dead.Advance(i)]=DeadCounter(Network);    
    Dead.All(i)=Dead.Normal(i)+Dead.Advance(i);
    
    Network=Leach_SelectCluster(Model,Network);
    HeadClusters(i)=sum([Network.IsCluster]);    
    Packet.ToBS(i)=HeadClusters(i);
        
    NotCluster=Model.NodeNumber-Dead.All(i)-HeadClusters(i);
    if((NotCluster)>0)
    Network=Leach_UsedEnergy(Model,Network);       
    Packet.ToCH(i)=NotCluster;
    end       
end
%% Result
Result.Dead=Dead;
Result.Networks=Networks;
Result.Packet=cumsum(Packet.ToCH+Packet.ToBS);
Dead.All(Dead.All==0)=Inf;
[~,Location]=min(Dead.All);
Result.FirstDead=Location;
end