function [ab,Mcap,kount]=sr_OPhrbimod(nofdig,N)
%SR_OPHRBIMOD Recurrence coefficients for the half-range
%bimodal weight function.
%   [AB,MCAP,KOUNT]=SR_OPHRBIMOD(NOFDIG,N) generates in
%   dig-digit arithmetic the first N recurrence coefficients 
%   for the weight function w(x)=(x^2)*exp[-(x^2-1)^2/(4*epsi)]
%   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.m.) The values of
%   dig and epsi have to be entered when prompted, and also 
%   the values of b10 and b20 (b10<b20), where x1=sqrt(1+b10) and
%   x2=sqrt(1+b20) are assumed to be close estimates of the two 
%   x-values for which w(x)=1/2. The output variables MCAP and
%   KOUNT are inherited from the routine smcdis.m.

global mc mp iq idelta irout DM AB dig epsi eps0 
mc=3; mp=0; iq=0; idelta=1; irout=1; DM=[0 0];
eps0=.5*10^(-nofdig); Mmax=2000;
dig=input('enter the value of dig')
digits(dig);
epsi=input('enter the value of epsi')
b10=input('enter the value of b10')
b20=input('enter the value of b20')
b1=xhrbimod(epsi,b10); b2=xhrbimod(epsi,b20);
b1=vpa(b1); b2=vpa(b2);
AB=[[0 b1];[b1 b2];[b2 Inf]];
[ab,Mcap,kount]=smcdis(dig,N,eps0,@squadgp,Mmax); ab=vpa(ab,nofdig);
