function ab=sr_OPbimod(nofdig,N)
%SR_OPBIMOD Recurrence coefficients for the bimodal weight
%function.
%   AB=SR_OPBIMODAL(NOFDIG,N) generates in dig-digit
%   arithmetic the first N recurrence coefficients for the 
%   weight function w(x)=exp[-(x^2-1)^2/(4*epsi)] on [-Inf,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 value of mu  
%   if epsi<1/10.

global eps0 epsi dig h
eps0=.5*10^(-nofdig);
dig=input('enter the value of dig')
epsi=input('enter the value of epsi')
if epsi<1/10
h0=exp(-1/(4*epsi));
%
% input mu larger than, or equal to 1
%
  mu=input('enter the value of mu')
  h=exp(mu)*h0;
end
ab=sOPbimod_gp(dig,N,epsi); ab=vpa(ab,nofdig);
