> f1 := 1/( 1 - x^2 -y^2 );
> with(plots);
> plot3d(f1, x=-2.0..2.0, y=-2.0..2.0);
> y1 := (1-x^2)^(1/2);
> plot(y1, x =-2.0..2.0);