function iterate_example(ics,N,b,bnds)
    nics = size(ics,1);
    for j = 1:nics,
        x = ics(j,1);
        y = ics(j,2);
        plot(x,y,'rd')
        for k = 1:N,
            x1 = (1-x)*x+b*y;
            y1 = y/2 + x;
            if (x1 < bnds(1) | x1 > bnds(2) | y1 < bnds(3) | y1 > bnds(4)) then
                break;
            end
            plot([x x1],[y y1],'r-d')
            x = x1;
            y = y1;
        end
    end
    a = gca();
    set(a,'data_bounds',bnds);
endfunction
    
