clear
clc
clf

aph(1,1)=0.7;
aph(2,1)=0.7;

for ite=2:1:500; 
ee=50000000*2;
    uz=ite^2/(ee+ite);
    aph(ite,1)=aph(ite-1,1)/(1+uz);
    % 1/ite scheme
    %    aph(ite,1)=aph(ite-1,1)/ite;
end

plot(aph)

%50500
%ee=5000000000000*1.2;

%1100
%ee=50000000*1.2;

%1500
%ee=500000000*7; slow
%ee=50000000*7; medium
%ee=5000000*7; fast alpha died at 800

%500
%ee=50000000*2;slow
%ee=5000000*2; medium
%ee=500000*2; fast alpha died at 250

