
function exact_integral_disk_chebfun

example=1;

switch example
    case 1
        f=@(x,y) (1*x+0.5*y).^10;
    case 2
        f=@(x,y) (1+1*x+0.5*y).^11;
end

fc=diskfun(f);
Ic=sum2(fc);
fprintf('\n \t Ic: %1.15e \n',Ic)
