%
% Computes the volume of the solid obtained by slicing the torus
%      (x - 2)^2 + y^2 =1
% parallel to the x-axis at t and -t.
%
function vol = volumeDisk(t)

% take the volume from 0 to t and double it to get the volume from 
% -t to t
vol = 2 * quad('areaSlice', 0, t);