\[\wp\left(u\right)=\left(\frac{\pi}{2\omega_{1}}\right)^{2}\left(\frac{\theta_{%
1}'''\left(0,q\right)}{3\theta_{1}'\left(0,q\right)}-\frac{{\mathrm{d}}^{2}}{{%
\mathrm{d}z}^{2}}\ln\theta_{1}\left(z,q\right)\right),\]