sol = NDSolve[{
\!\(
\*SubscriptBox[\(\[PartialD]\), \(\[Theta]\)]\(v[\[Theta],
r]\)\) == -\!\(
\*SubscriptBox[\(\[PartialD]\), \({r, 2}\)]\(u[\[Theta], r]\)\) - 2,
\!\(
\*SubscriptBox[\(\[PartialD]\), \(\[Theta]\)]\(u[\[Theta], r]\)\) ==
v[\[Theta], r],
u[0, r] == 1,
u[\[Theta], 0] == 1 - \[Theta]^2, u[\[Theta], 1] == 1 - \[Theta]^2,
v[0, r] == 0, v[\[Theta], 0] == -2 \[Theta],
v[\[Theta], 1] == -2 \[Theta]}
{u, v},
{r, 0, 1},
{\[Theta], 0, 1}, MaxStepSize -> 0.001]