%SR_FREUD1
%
function ab=sr_freud1(dig,dig0,N,mu,nu)
digits(dig);
mom=smomfreud(dig,N,mu,nu);
ab0=schebyshev(dig,N,mom');
mom=smomfreud(dig+10,N,mu,nu);
ab1=schebyshev(dig+10,N,mom');
err=max(subs(abs((ab1(:,2)-ab0(:,2))./ab1(:,2))))
ab=vpa(ab0,dig0)
