function ab=sr_meixner_ops(dig,N,b,c)
%SR_MEIXNER_OPS Recurrence coefficients for the Meixner 
%orthogonal polynomials.
%   AB=SR_MEIXNER_OPS(DIG,N,B,C) generates in DIG-digit
%   arithmetic the Nx2 array AB of the first N recurrence 
%   coefficients for the weight function 
%   w(k)=C^k*gamma(k+B)/(gamma(B)*k!), k=0:Inf, B>0, 0<C<1.

digits(dig); 
ab=vpa(zeros(N,2)); sb=vpa(b); sc=vpa(c); s1=vpa(1);
k=vpa(1:N);
ab(:,1)=((s1+sc)*(k-1)+sb*sc)/(s1-sc);
k=vpa(1:N-1);
ab(:,2)=[(s1-sc)^(-sb) (sc/((s1-sc)^2))*k.*(k+sb-s1)]';
