function ab=r_mpplusdis(N,a,b,ty)
%R_MPPLUSDIS Recurrence coefficients for the Marchenko-Pastur
%measure plus a discrete measure
%   R_MPPLUSDIS(N,A,B,TY) generates the first N recurrence
%   coefficients for the Marchenko-Pastur measure with 
%   parameters A, B, C and with a discrete measure whose
%   abscissa and jump are given by TY added to it.

global mc mp iq idelta irout DM abmp
mc=1; mp=1; iq=1; idelta=2; irout=1; DM=ty;
eps0=.5e-12; Mmax=N+1;
[ab,Mcap,err,kount]=mcdis(N,eps0,@quad_mpplusdis,Mmax);
Mcap, err, kount
