magma://Please install CHIMP (https://github.com/edgarcosta/CHIMP) if you want to run this code
chi := DirichletCharacter("209.45");
S:= CuspForms(chi, 4);
N := Newforms(S);
sage:from sage.modular.dirichlet import DirichletCharacter
H = DirichletGroup(209, base_ring=CyclotomicField(6))
chi = DirichletCharacter(H, H._module([0, 2]))
N = Newforms(chi, 4, names="a")
This newform subspace can be constructed as the kernel of the linear operator
\( T_{2}^{50} - 2 T_{2}^{49} + 155 T_{2}^{48} - 272 T_{2}^{47} + 13477 T_{2}^{46} - 21509 T_{2}^{45} + \cdots + 10\!\cdots\!04 \)
acting on \(S_{4}^{\mathrm{new}}(209, [\chi])\).