
function exact_integral

d=3;

switch d
    case 2
f=@(x,y) (0.3*x+0.9*y).^10;
I=integral2(f,-1,1,-1,1,'AbsTol',10^(-15),'RelTol',10^(-15));
    case 3
        f=@(x,y,z) (0.3*x+0.9*y+0.8*z).^10;
        I=integral3(f,-1,1,-1,1,-1,1,'AbsTol',10^(-15),'RelTol',10^(-15));
end


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