clear;
numer = [0 1];
denom = [1 1];
for i = 1:10
    for j = 1:length(numer)-1;
        new_numer(j) = numer(j+1) + numer(j);
        new_denom(j) = denom(j+1) + denom(j);
    end
    for k = 1:length(numer)
        temp_numer(2*k-1) = numer(k);
        temp_denom(2*k-1) = denom(k);
    end
    for l = 1:length(new_numer)
        temp_numer(2*l) = new_numer(l);
        temp_denom(2*l) = new_denom(l);
    end
    numer = temp_numer;
    denom = temp_denom;
end
rat = numer./denom;
plot(rat,'.k')