function area = areaTorus2( t )

% The areaTorus2 function computes the cross-sectional area of the torus
% at x = t using the method of shells.

y = sqrt( 1 - (t-2).^2 ) ;
height = 2.*y ;
area = 2*pi*height.*t ;