figure(1); clf;
h = [];
i = 1;
for a = 1.5:-.01:0
    p_a = 1/a;
    h(i) = log(a);
    rectangle('Position', [0 0 a p_a]);
    hold on;
    i = i+1;
end
%%
h_2 = [];
i = 1;
for a = 1:-.0001:0
    p_a = 1/a;
    h_2(i) = log(a);
    i = i+1;
end
