function [ab,Mcap,kount]=sr_OPhrbimod(dig,nofdig,N,b10,b20)
%SR_OPHRBIMOD Recurrence coefficients for the half-range
%bimodal weight function.
%   [AB,MCAP,KOUNT]=SR_OPHRBIMOD(DIG,NOFDIG,N,B10,B20) 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. The input variables
%   B10 and B20, B10<B20, are close estimates of the two x-values
%   for which w(x)=1/2. The output variables MCAP and KOUNT are
%   inherited from the rotuine mcdis.m.

global mc mp iq idelta irout DM AB epsi eps0 
digits(dig);
mp=0; mc=3; iq=0; idelta=1; irout=1; DM=[0 0];
eps0=.5e-12; Mmax=2000;
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);
