
function demo_disk_notadaptive

example=3;

switch example
    case 1
        intfcn=inline('y.*sin(x)+x.*cos(y)','x','y');
        intfcn_dblquad=inline('(y.*sin(x)+x.*cos(y)).*(x.^2+y.^2<=1)','x','y');
    case 2
        intfcn=inline('3*x','x','y');
        intfcn_dblquad=inline('(3*x).*(x.^2+y.^2<=1)','x','y');
    case 3
        intfcn=inline('exp(x)+exp(y)','x','y');
        intfcn_dblquad=inline('(exp(x)+exp(y)).*(x.^2+y.^2<=1)','x','y');
end



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);

f_pts=feval(intfcn,x,y);

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

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