function [ab,Mcap,kount]=sr_morse(nofdig,N)
%SR_MORSE Recurrence coefficients for the Morse weight
%function.
%   [AB,MCAP,KOUNT]=SR_MORSE(NOFDIG,N) generates in dig-digit
%   arithmetic the first N recurrence coefficients for the 
%   weight function w(x)=exp(-5*x+6*(exp(-x)-1)) on [0,Inf]
%   and places them to NOFDIG digits into the Nx2 array AB.
%   (Be sure to specify the correct weight function wf.m at
%   the end of the routine squadgp.) The value of dig has to
%   be entered when prompted. The output variables MCAP and 
%   KOUNT are inherited from the routine smcdis.

global eps0 dig
eps0=.5*10^(-nofdig);
dig=input('enter the value of dig')
ab=smorse(dig,N); ab=vpa(ab,nofdig);
