function a=verhulst(a1, q, n)
    a=a1;
    an=a1;
    for i=2:n
        an=q.*an.*(-an+1.0);
        a=[a; an];
    end
endfunction
