clear;clc;
a=2;b=3;c=4;n=5;%自己改
[x,y,z]=meshgrid(linspace(-a,a),linspace(-b,b),linspace(-c,c));
v=x.^2/a^2+y.^2/b^2+z.^2/c^2-1;
isosurface(x,y,z,v,0);%只要骨架的这句去掉
h=contourslice(x,y,z,v,linspace(-a,a,n),linspace(-b,b,n),linspace(-c,c,n),[0 0]);
grid on;xlabel('x');ylabel('y');zlabel('z');axis equal;view(3);
title(sprintf('$$\\frac{x^2}{%d^2}+\\frac{y^2}{%d^2}+\\frac{z^2}{%d^2}=1$$',a,b,c),'interpreter','latex')