function vol = volumeTorus3( t )

% The volumeTorus3 function computes the volume of the torus from 
% y = -t to y = t.

vol = 2*quad( 'areaTorus3', 0, t ) ;