
figure;
z=[5 20 100];
for w=1:length(z)
A=POP{1,z(w)};
A1=cell2mat(A);
solu=SOLUTION{z(w)};
solu1=reshape(solu,[1,(no_of_fam*4)]);
k=1;
    for j=1:2:4
        for i=1:no_of_fam
            A2(k)=A1(i,j);
            k=k+1;
        end
    end
    
    k=1;
    for j=2:2:4
        for i=1:no_of_fam
            A3(k)=A1(i,j);
            k=k+1;
        end
    end
    
    subplot(1,3,w);scatter3(A2,A3,solu1(1:80),'*b')
    axis([-10 10 -10 10 0 20])
    set(gca,'xtick',[-10:5:10])
    set(gca,'ytick',[-10:5:10])
end

figure;
z=[5 20 100];
for w=1:length(z)
x=-10:0.1:10;
y=-10:0.1:10;
for i=1:length(x)
    for j=1:length(y)
        s=[x(i) y(j)];
        solu11(i,j)=xyfunc_d(s);
    end
end
subplot(1,3,w);contour(x,y,solu11);
hold on
scatter(BEST_SOLUTION{1,z(w)}(1,1),BEST_SOLUTION{1,z(w)}(2,1),'*r');
end

figure;
z=[5 20 100];
for w=1:length(z)
    con=functional_values(1:z(w));
    subplot(1,3,w);
    plot(con);
end





























