function koule(r)

if(nargin < 1)
   r = 1;
end

theta = 0:0.3:pi;
fi = 0:0.3:2*pi;

[FI,THETA] = meshgrid(fi,theta);

x = r*sin(THETA).*cos(FI);
y = r*sin(THETA).*sin(FI);
z = r*cos(THETA);

surf(x,y,z), axis equal, shading flat
