
%----------------------------------------------------------
% In this driver, we use high order cubature to test the 
% orthogonality of Logan-Shepp polynomials.
%----------------------------------------------------------

[I,err]=cubature_disk('logan_shepp_product', 10^(-10), 4, 3);

fprintf('\n \n');