function area = areaSphere( x) 

% This function computes the area of a half a cross section of a sphere.
rad = sqrt( 9 - x.^2) ;
area = 0.5*pi*rad.^2 ;