function vol = volumeTorus3( s )
% This function computes the volume of the torus from -s to s.

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