plot3(Unstablemanifold(:,1),Unstablemanifold(:,3),Unstablemanifold(:,2),'b.','markersize',5)

psection=zeros(length(Unstablemanifold(:,1)),3);
tend=100;
for i=1:length(Unstablemanifold(:,1))
    init=Unstablemanifold(i,:);
    [t,solX,a, b, c]=ode15s('chemostatmodelus',[0,tend],init,options,parameters);
                psection(i,:)=b;
                strt=1;
                t=t(strt:end,1);
                xx=solX(strt:end,1);
                yy=solX(strt:end,2);
                zz=solX(strt:end,3);
                if mod(i,3)==0
                plot3(xx,zz,yy,'b-','linewidth',1,'markersize',5)
                hold on
                end
end
plk=1;
psection1=psection(plk:end,:);

psectionfamily=psection1;

