function ab=sr_hermite(dig,N,mu)
% SR_HERMITE Symbolic counterpart of r_hermite
% 
if nargin<3, mu=0; end
if N<=0|mu<=-1/2, error('parameter(s) out of range'), end
m0=gamma(vpa(mu+1/2));
if N==1, ab=[0 m0]; return, end
N=N-1; n=1:N; nh=vpa(n/2); nh(1:2:N)=nh(1:2:N)+mu;
A=zeros(1,N+1); B=[m0 nh];
ab=[A' B'];
