
function demo_disk_zernike


thetamin=0; thetamax=2*pi; rmin=0; rmax=1; degree_GL=32; cub_type=4;

[x,y,rep_r,r_weights,theta_weights]=integration_annulus(...
    thetamin, thetamax, rmin, rmax, degree_GL, cub_type);


% [theta,r] = cart2pol(X,Y);
%       idx = r<=1;
%       z = nan(size(X));
%       z(idx) = zernfun(5,1,r(idx),theta(idx));

[thz,rz]=cart2pol(x,y);
f_pts=zernfun(5,1,rz,thz,'norm').*zernfun(5,1,rz,thz,'norm');

I=integral_evaluation(x,y,f_pts,rep_r,r_weights,theta_weights,degree_GL);

fprintf('\n \n  \t >> demo_disk_zernike \n')
fprintf('\n \t INTEGRAL: %1.15e \n \n',I);
