% ACC_ISTARAL  Accuracy of Istaral
%
function y=acc_Istaral(alpha,beta)
global N abj abj1 abj2 abj3 abj4 abl1 abl2 abl3 abl4 
N=50;
abj=r_jacobi01(N,alpha,beta);
abj1=r_jacobi01(N,alpha,0); abj2=r_jacobi01(N,beta,0);
abj3=r_jacobi01(N,alpha,2*beta+1); 
abj4=r_jacobi01(N,beta,2*alpha+1);
abl1=r_jacobilog(N,alpha,0); abl2=r_jacobilog(N,beta,0);
abl3=r_jacobilog(N,alpha,2*beta+1);
abl4=r_jacobilog(N,beta,2*alpha+1);
i=0;
for n=N:-1:N-1
  i=i+1;
  [Istaralg,Istarlog]=Istaral(n);
  Ist(i)=Istarlog;
end
y=-log10(abs((Ist(2)-Ist(1))/Ist(2)));
