
clear 
format long


    %how many parameter values?
    numPar=500;

    %final time
    tf=300;
    
    %Time axis
    t=linspace(0, tf, tf+1);

    %parameters
    r=linspace(1.9, 2.9, numPar);
    
    %an initial condition
    x_initial=0.6;
    
    for i=1:numPar
        X(i,1:(tf+1))=populationSeries(x_initial,tf,r(i));
    end %end of numPar loop    
        
    %plot
    figure
    hold on
    for j=1:numPar
        thisParameter=r(j)*ones(1,tf+1-200);
        plot(thisParameter(:),X(j,200:tf),'.')
    end %end of plot loop    
        