clear all
clc

meanPg=0;
test=10000
conv=0;
kk=0;

for k=1:test


% x=0:0.01:20;
% y=(x.^2).*sin(x);
% figure(1);plot(x,y)
% 
% [Min,Cmin]=min(y);
% Pmin=x(Cmin);
 
w=0.7;c1=0.7;c2=0.7;gen=20;nbrP=50;
 
a=0;b=20;
 
xi=a + (b-a).*rand(1,nbrP);

vi=zeros(1,nbrP);

fi1=rand(1,nbrP);

Pi=xi;

fi2=rand(1,nbrP);

Fpi=(Pi.^2).*sin(Pi);


[Min,Cmin]=min(Fpi);
Pg=Pi(Cmin);

vi=(vi*w)+((c1*fi1).*(Pi-xi))+ ((c2*fi2).*(Pg-xi));

xi=xi+vi;

for j=1:gen-1

Fxi=(xi.^2).*sin(xi);

for i=1:nbrP
if Fxi(i)<Fpi(i)
    Pi(i)=xi(i);
    Fpi(i)=Fxi(i);
end
end

[Min,Cmin]=min(Fpi);
Pg=Pi(Cmin);



fi1=rand(1,nbrP);
fi2=rand(1,nbrP);
vi=(vi*w)+((c1*fi1).*(Pi-xi))+ ((c2*fi2).*(Pg-xi));
xi=xi+vi;


for i=1:nbrP
if xi(i)<a || xi(i)>b 
    xi(i)=a + (b-a).*rand(1,1);
end
end



end

Fxi=(xi.^2).*sin(xi);

for i=1:nbrP
if Fxi(i)<Fpi(i)
    Pi(i)=xi(i);
    Fpi(i)=Fxi(i);
end
end

[Min,Cmin]=min(Fpi);
Pg=Pi(Cmin);



if Fpi(Cmin)<-290
    conv=conv+1;
    meanPg=meanPg+Pg;
    kk=kk+1;
    FPg=Fpi(Cmin);
    Vpg(kk)=Pg;
    Vfpg(kk)=FPg;
end

end

meanPg=meanPg/kk
conv=(conv/test)*100
POURcopt=meanPg/17.3932