
function exact_integral_disk

example=1;

switch example
    case 1
        fpolar=@(r,t) (1*r.*cos(t)+0.5*r.*sin(t)).^10.*r;
        I=integral2(fpolar,0,1,0,2*pi,'AbsTol',10^(-15),'RelTol',10^(-15));
    case 2
        fpolar=@(r,t) (1+1*r.*cos(t)+0.5*r.*sin(t)).^11.*r;
        I=integral2(fpolar,0,1,0,2*pi,'AbsTol',10^(-15),'RelTol',10^(-15));
end

fprintf('\n \t I : %1.15e',I)