%
% Computes the volume of the solid obtained by cutting the torus along
% the x-axis at t.
%
function vol = volumeTorus(t)

vol = quad('areaShell', 1, t);