Finding Volumes of Revolution

Since this involves straightforward integration, one may simply issue commands such as:

  > f:= x-> sqrt(x) +1;

  > vol:= int(Pi*f(x)^2, x=0..9);

  > evalf(vol);