Section 13.5 #37

Picture and Double Integral

>	with(plots):

Warning, the name changecoords has been redefined

>	cylinder1:=implicitplot3d(x^2+z^2=1,x=-1.5..1.5,y=-1.5..1.5,z=-1.5..1.5,axes=boxed,scaling=constrained):

>	cylinder2:=implicitplot3d(y^2+z^2=1,x=-1.5..1.5,y=-1.5..1.5,z=-1.5..1.5,axes=boxed,scaling=constrained):

>	display(cylinder1,cylinder2);

Here we compute the desired surface area.

>	16*Int(Int(1/sqrt(1-x^2),y=0..x),x=0..1);

>

value(%);

>