CODE:

function [xx,yy,zz] = mysphere(varargin)

%SPHERE Generate sphere.

%   [X,Y,Z] = SPHERE(N) generates three (N+1)-by-(N+1)

%   matrices so that SURF(X,Y,Z) produces a unit sphere.

%

%   [X,Y,Z] = SPHERE uses N = 20.

%

%   SPHERE(N) and just SPHERE graph the sphere as a SURFACE

%   and do not return anything.

%

%   SPHERE(AX,...) plots into AX instead of GCA.

%

%   See also ELLIPSOID, CYLINDER.

%   Clay M. Thompson 4-24-91, CBM 8-21-92.

%   Copyright 1984-2002 The MathWorks, Inc.

%   $Revision: 5.8.4.1 $  $Date: 2002/09/26 01:55:25 $

% Parse possible Axes input

error(nargchk(0,2,nargin));

[cax,args,nargs] = axescheck(varargin{:});

r = 0.5; %我把半径改为0.5，LZ的半径0.0005

n = 20;

if nargs > 0, n = args{1}; end

% -pi <= theta <= pi is a row vector.

% -pi/2 <= phi <= pi/2 is a column vector.

theta = (-n:2:n)/n*pi;

phi = (-n:2:n)'/n*pi/2;

cosphi = cos(phi); cosphi(1) = 0; cosphi(n+1) = 0;

sintheta = sin(theta); sintheta(1) = 0; sintheta(n+1) = 0;

x = r * cosphi*cos(theta);

y = r * cosphi*sintheta;

z = r * sin(phi)*ones(1,n+1);

if nargout == 0

cax = newplot(cax);

surf(x,y,z,'parent',cax)

else

xx = x; yy = y; zz = z;

end