magma://Please install CHIMP (https://github.com/edgarcosta/CHIMP) if you want to run this code
chi := DirichletCharacter("800.63");
S:= CuspForms(chi, 2);
N := Newforms(S);
sage:from sage.modular.dirichlet import DirichletCharacter
H = DirichletGroup(800, base_ring=CyclotomicField(20))
chi = DirichletCharacter(H, H._module([10, 0, 19]))
N = Newforms(chi, 2, names="a")
This newform subspace can be constructed as the kernel of the linear operator
\( T_{3}^{64} + 6 T_{3}^{61} - 111 T_{3}^{60} + 148 T_{3}^{59} - 42 T_{3}^{58} + 358 T_{3}^{57} + \cdots + 110166016 \)
acting on \(S_{2}^{\mathrm{new}}(800, [\chi])\).