function ab=r_blockhahn(n,N)
%R_BLOCKHAHN
%
global mc mp iq idelta 
mc=2*N+1; mp=0; iq=1; idelta=2; eps0=1; Mmax=200;
[ab,Mcap,kount]=mcdis(n,eps0,@quad_blockhahn,Mmax);
