I have a lattice, $\Omega$, with basis $\lbrace 2+i, 1+3i\rbrace$ and fundamental region (square) $P$ with vertices $1+2i,\ 2, \ -1+i\ $ and $-i$. I want to find the zeroes and poles of $\wp$ and $\wp^{'}$ in $P$, along with their orders; where $$\wp=\frac{1}{z^2}+\sum_{\omega\in \Omega\backslash\lbrace0\rbrace}\left(\frac{1}{\left(z-\omega\right)^2}-\frac{1}{\omega^2}\right),$$ and $$\wp^{'}=-\frac{2}{z^3}+\sum_{\omega\in \Omega\backslash\lbrace0\rbrace}\left(\frac{-2}{\left(z-\omega\right)^3}\right)=-2\sum_{\omega\in \Omega}\frac{1}{\left(z-\omega\right)^3}.$$
I'd appreciate any guidance you may have to offer.
For any lattice, the poles of $\wp'$ are triple at lattice points, and its zeros are single at half-periods.
For a square lattice, with periods $\omega$ and $i\omega$, the poles of $\wp$ are double at lattice points. Its zeros are also double, at $\frac12(1+i)\omega$ modulo the lattice.