magma://Please install CHIMP (https://github.com/edgarcosta/CHIMP) if you want to run this code
chi := DirichletCharacter("5001.1");
S:= CuspForms(chi, 2);
N := Newforms(S);
sage:from sage.modular.dirichlet import DirichletCharacter
H = DirichletGroup(5001, base_ring=CyclotomicField(2))
chi = DirichletCharacter(H, H._module([0, 0]))
N = Newforms(chi, 2, names="a")
This newform subspace can be constructed as the kernel of the linear operator
\( T_{2}^{60} + 13 T_{2}^{59} + 3 T_{2}^{58} - 664 T_{2}^{57} - 2183 T_{2}^{56} + 14182 T_{2}^{55} + \cdots + 8551 \)
acting on \(S_{2}^{\mathrm{new}}(\Gamma_0(5001))\).