clear?all;?
clc;?
tic?
global?mvab?kvb?busdata?linedata?nf?ng?Ntiesw?Lsw?Lswcnt?Fbest;?
data16;?
Xb=[5?11?16];?
[TPLVolt]=loadflow（Xb）;?
TPLold=TPL;?
Fbest=TPL;?
for?i=1:ng?%maxno.of?iterations?????
????[SolutSet]=searchdomain（Xb）;?
?????u=length（SolutSet（:1））;?
?????for?j=1:u?%length（SolutSet（:1））?
?????????clear?TPL?Volt;??
?????????[TPLVolt]=loadflow（SolutSet（j:））;?
?????????g1（j）=TPL;?
?????????g2（j:）=Volt;?
?????end?
?????g1;?
?????fbest=min（g1）;?
?????k=min（find（fbest==g1））;?
?????Xbest=SolutSet（k:）;?
?????if（fbest?????????Fbest=fbest;?
?????????Xb=Xbest;?
??