find the volume of the region bounded between the planes x+y+2z=2 and 2x+2y+z=4 in the first octant. in matlab

main.m
syms x y z
f1 = x + y + 2*z - 2;
f2 = 2*x + 2*y + z - 4;

% Solve the equations for the intersection points
sol = solve(f1 == 0, f2 == 0, x, y, z);
sol_x = sol.x; % x values
sol_y = sol.y; % y values
sol_z = sol.z; % z values

% Define the limits of integration
x_limits = [0 sol_x];
y_limits = [0 sol_y];
z_limits = [0 sol_z];

% Calculate the volume using triple integral 
volume = integral3(@(x,y,z)1, x_limits, y_limits, z_limits);

disp("The volume of the region is " + volume)
483 chars
20 lines