I am trying to learn how to solve Igusa type local zeta function. Ex. $$\int_{\mathbb{Z}_{p}}||x^3,x^2y,y^2||d\mu(x,y)$$ A nice method I was recently introduced to was to substitute $x=a+px'$ and $y=b+py'$, where $a$ and $b \in \{0,1,\cdots,p-1\}$ . Hence, we should obtain $$\frac{1}{p^2}(\sum_{a=0}^{p-1}\sum_{b=0}^{p-1}\int_{\mathbb{Z}_p}||(a+px)^3,(a+px)^2(b+py),(b+py)^2||d\mu(x,y)) \\ = \frac{1}{p^2}(p^2-1)+p^{2s-2}\int_{\mathbb{Z}_p}||px^3,px^2y,y^2||d\mu(x,y).$$
Here I am trying to find a trick I can use to solve it.
Nice, I solve it using your method and I checked my answer using Zeta SageMath.
Here is my answer, first
$$\int_{\mathbb{Z}_p^2}||x^3,x^2y,y^2||_p^sd\mu(x,y) = \int_{\mathbb{Z}_p^2}||x^3,y^2||_p^sd\mu(x,y).$$
Since, if $|x^2y|_p^s\ge|x^3|_p^s \Rightarrow |y|_p^s\ge|x|_p^s\Rightarrow |y^2|_p^s\ge|xy|_p^s\ge|x^2y|_p^s$, since $|x|_p^s\le 1$. $$ $$ Similarly, if $|x^2y|_p^s\ge|y^2|_p^s\Rightarrow |x^2|_p^s\ge|y|_p^s\Rightarrow|x^3|_p^s\ge|xy|_p^s\ge|x^2y|_p^s$. Hence, $$\int||x^3,y^2||_p^s d\mu(x,y)= 1-p^{-2}+p^{-2s-2}\int||px^3,y^2||_p^s d\mu(x,y)\cdots(1)$$ Now consider, $$\int||px^3,y^2||_p^s d\mu(x,y)=1-p^{-1}+p^{-s-1}\int||x^3,py^2||_p^s d\mu(x,y)\cdots(2)$$ Similarly, $$\int||x^3,py^2||_p^s d\mu(x,y)=1-p^{-1}+p^{-s-1}\int||p^2x^3,y^2||_p^s d\mu(x,y)\cdots(3)$$ Finally, $$\int||p^2x^3,y^2||_p^s d\mu(x,y)=1-p^{-1}+p^{-2s-1}\int||x^3,y^2||_p^s d\mu(x,y)\cdots(4).$$ Now all I have to do is to put (4) in (3) in (2) in (1) and obtain $$L= 1-p^{-2}+p^{-2s-2}[1-p^{-1}+p^{-s-1}[1-p^{-1}+p^{-s-1}[1-p^{-1}+p^{-2s-1}[L]]]],$$ where $L = \int||x^3,y^2||_p^s d\mu(x,y)$. Hence, $$ L = \frac{1-p^{-2}+(1-p^{-1})(p^{-2}T^2+p^{-3}T^3+p^{-4}T^{4})}{1-p^{-5}T^{6}},$$ where $T=p^{-s}$.
Thank you very much for this method it made it much easier. I still need to learn another method which uses vertex geometry, it is more general but pain in the head :)