function sM = sSuperQuad(a,b,c)
% sSuperQuad - spline demo applied to superquadrics
%	sM = sSuperQuad(a,b,c)

if( nargin == 0 )
   [x1,y1,z1] = superquad;
else
   [x1,y1,z1] = superquad(a,b,c);
end;



t = 1:prod(size(x1));
  x = x1(t); y = y1(t); z = z1(t);
[az, el, r] = cart2sph(x,y,z);

sM = sFit(az,el, r);
sPlotSurf(sValueM(sM));
figure(gcf+1);
surf(x1,y1,z1);

return
