% SR_JACOBI01 Symbolic counterpart of r_jacobi01.m % function ab=sr_jacobi01(dig,N,a,b) if nargin<3, a=0; end; if nargin<4, b=a; end if((N<=0)|(a<=-1)|(b<=-1)) error('parameter(s) out of range'), end sa=vpa(a,dig); sb=vpa(b,dig); cd=sr_jacobi(dig,N,a,b); n=1:N; ab(n,1)=(1+cd(n,1))./2; ab(1,2)=cd(1,2)/2^(sa+sb+1); n=2:N; ab(n,2)=cd(n,2)./4;