function vol = volumeTorus2( t ) 

% The volumeTorus2 function computes the volume of the torus from 
% x = 1 to x = t.

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