How to creat 3D mesh about superellipse hole in a cube

The equation of the superellipse is (x/a)^p+(y/b)^p+(z/c)^p=1. p<2.a,b,c are constant. I want to draw a 3D mesh of the cube include this superellipse. I download iso2mesh which can generate 3D mesh but not succeed.

Answer 1

[KSSV]

th = linspace(-pi/2,pi/2) ;
phi = linspace(-pi,pi) ;

[T,P ] =meshgrid(th,phi) ;

a = 2; b = 2; c = 1 ;
X = a*cos(T).*cos(P) ;
Y = b*cos(T).*sin(P) ;
Z = c*sin(T) ;

surf(X,Y,Z)

Read here: https://en.wikipedia.org/wiki/Ellipsoid

[KOU DU]

thank you for you answer. But I don't want to mesh a ellipse.You can see in my question, I give the equation of the hole: (x/a)^p+(y/b)^p+(z/c)^p=1. I want to draw a cube with a hole in the center.