how do I compute the fundmental domain for a congruence subgroup of $SL(2, \mathbb{Z})$
This region is important because the theta function $\theta(z) = q^{n^2}$ is invariant under two transformations $z \mapsto z + 1$ and $z \to - \frac{1}{4z}$, and these two generate the congruence subgroup $\Gamma_0(4)$.
I know there is an algorithm for any Fuchsian group. Here is an erroneous comutation of the fundamental domain. It has too many cusps.

I found a blog with a picture of what could be the fundamental domain of $\Gamma_0(2)$. Therefore the domain for $\Gamma_0(4)$ could be related. I know that $[SL(2,\mathbb{Z}): \Gamma_0(2)]= 3$ and $[SL(2,\mathbb{Z}): \Gamma_0(4)]= 6$ but I can't figure out which two copies to join. And ther's no derivation for this particular case.
