function vol = volumeTorus2( x) 

% This function computes the volume of the torus from 1 to X.

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