D((x,y,z)) <= c0 * ( (x^2+y^2)^2 + z^2 )^{1/4}