%R_MARCHENKO_PASTUR_DIS

function [ab,Mcap,err,kount]=r_marchenko_pastur_dis(N,a,b,c)
global mc mp iq idelta irout DM abmp
%if c<=1
%  error('the parameter c is out of range');
%end
[mom,mu]=smom_marchenko_pastur(20,1,a,b,c);
abmp=r_marchenko_pastur_mod(110,100,a,b,c); abmp(1,2)=1/c;
ty=[0,1-1/c];
mc=1; mp=1; iq=1; idelta=2; irout=1; DM=ty;
eps0=.5e-13; Mmax=100; %Mmax=N+1;
[ab,Mcap,err,kount]=mcdis(N,eps0,@quad_marchenko_pastur_dis,Mmax);
