Modular equations of order $p$ and Theta functions

An algorithm to obtain equations between theta functions with integral characteristics evaluated at $\tau$ and $p\tau$ for $g>1$ is presented.