%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%??????????????????????????????????????????????????????????????????????%
%?SEP:?A?Stable?Election?Protocol?for?clustered????????????????????????%
%??????heterogeneous?wireless?sensor?networks??????????????????????????%
%??????????????????????????????????????????????????????????????????????%
%?（c）?Georgios?Smaragdakis?????????????????????????????????????????????%
%?WING?group?Computer?Science?Department?Boston?University???????????%
%??????????????????????????????????????????????????????????????????????%
%?You?can?find?full?documentation?and?related?information?at:??????????%
%?http://csr.bu.edu/sep????????????????????????????????????????????????%
%??????????????????????????????????????????????????????????????????????%??
%?To?report?your?comment?or?any?bug?please?send?e-mail?to:?????????????%
%?gsmaragd@cs.bu.edu???????????????????????????????????????????????????%
%??????????????????????????????????????????????????????????????????????%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%??????????????????????????????????????????????????????????????????????%
%?This?is?the?SEP?[1]?code?we?have?used.???????????????????????????????%
%??????????????????????????????????????????????????????????????????????%
%?[1]??Georgios?Smaragdakis?Ibrahim?Matta?and?Azer?bestavros?????????%??????
%??????“SEP:?A?Stable?Election?Protocol?for?clustered??????????????????%
%??????heterogeneous?wireless?sensor?networks“????????????????????????%
%??????Second?International?Workshop?on?Sensor?and?Actor?Network???????%
%??????Protocols?and?Applications?（SANPA?2004）Boston?MA?August???????%???????????????					
%??????2004.???????????????????????????????????????????????????????????%
%??????????????????????????????????????????????????????????????????????%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%?PARAMETERS?%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Field?Dimensions?-?x?and?y?maximum?（in?meters）
xm=100;
ym=100;

%x?and?y?Coordinates?of?the?Sink
sink.x=0.5*xm;
sink.y=0.5*ym;

%Number?of?Nodes?in?the?field
n=100

%Optimal?Election?Probability?of?a?node
%to?become?cluster?head
p=0.1;

%Energy?Model?（all?values?in?Joules）
%Initial?Energy?
Eo=0.2;
%Eelec=Etx=Erx
ETX=50*0.000000001;
ERX=50*0.000000001;
%Transmit?Amplifier?types
Efs=10*0.000000000001;
Emp=0.0013*0.000000000001;
%Data?Aggregation?Energy
EDA=5*0.000000001;

%Values?for?Hetereogeneity
%Percentage?of?nodes?than?are?advanced
m=0.1;
%\alpha
a=1;

%maximum?number?of?rounds
rmax=1100

%%%%%%%%%%%%%%%%%%%%%%%%%?END?OF?PARAMETERS?%%%%%%%%%%%%%%%%%%%%%%%%

%Computation?of?do
do=sqrt（Efs/Emp）;?

%Creation?of?the?random?Sensor?Network
figure（1）;
for?i=1:1:n
????S（i）.xd=rand（11）*xm;
????XR（i）=S（i）.xd;
????S（i）.yd=rand（11）*ym;
????YR（i）=S（i）.yd;
????S（i）.G=0;
????%initially?there?are?no?cluster?heads?