%LAGPLUSDIS  Laguerre measure plus a discrete measure
%
function [ab,Mcap,kount]=r_lagplusdis(N,alpha,ty)
global mc mp iq idelta irout DM AB 
global a
a=alpha;
mc=1; mp=size(ty,1); iq=1; idelta=2;
Mmax=N+1;
DM=ty; AB=[0,Inf]; eps0=1e3*eps;
[ab,Mcap,kount]=mcdis(N,eps0,@quadlp,Mmax);
