e = linspace(0,.3,10);
e = e(2:length(e));
th_max = zeros(1,length(e));

for k = 1:length(e)
    th_max(k) = orbit(e(k));
end

plot(e,th_max,'o');