In[]:=
RegionPlot3D[0<=z<=4&&
x^2+y^2
<=(4-z)^316,{x,-5,5},{y,-5,5},{z,-5,5}]