function [pso G] = pso_3D()
%region=zeros(3,2);
d=[-3,3;-3,3;-3,3];
a=[1,2,3;2,3,4;3,4,5;1,3,5;-2,-1,1;-3,-2,-1;-2,-1,2];
for j=1:7;
if Region_in(a(j,:),d)
    pso=1;
end
end


function flag=Region_in(b,c)
flag=1
for j=1:3
flag=flag&(b(1:j)>=c(j,1))&(b(1:j)<=c(j,2));
end