function area = areaTorus3( t )

% The areaTorus3 function computes the cross-sectional area of the torus
% at y = t using the method of washers.

rad1 = 2 - sqrt( 1 - t.^2 ) ;
rad2 = 2 + sqrt( 1 - t.^2 ) ;

area = pi*( rad2.^2 - rad1.^2 ) ;