%PLOT_RAM
%
IR=zeros(101,4);
j=0;
for N=10:10:40
  i=0; j=j+1
  for t=0:100
    i=i+1
    IR(i,j)=Ramanujan(dig,N,t);
  end
  for k=1:3
    err=max(abs(IR(:,k+1)-IR(:,k)))
  end
end
t=(0:100)';
plot(t,IR(:,4),'LineWidth',2)
axis([-10 110 0 1.1])
set(gca,'FontSize',15)
xticks([0 20 40 60 80 100])
yticks([0 .2 .4 .6 .8 1])
