with(plots):
implicitplot3d(z=x^2-y^2,x=-3..3, y=-3..3, z=-10..10);
plot3d(x^2-y^2,x=-3..3,y=-3..3);
implicitplot3d(y^2-x^2-z^2=16,x=-10..10,y=-20..20, z=-10..10,axes=normal);
plot3d([sqrt(y^2-x^2-16), -sqrt(y^2-x^2-16)], x=-10..10, y=-10..10, grid=[60,60]);