\[\wp\left(z\right)=\frac{\pi^{2}}{4\omega_{1}^{2}}\left(-\frac{1}{3}+{\csc}^{2}%
\left(\frac{\pi z}{2\omega_{1}}\right)+8\left(1-\cos\left(\frac{\pi z}{\omega_%
{1}}\right)\right)q^{2}+O\left(q^{4}\right)\right),\]