magma://Please install CHIMP (https://github.com/edgarcosta/CHIMP) if you want to run this code
chi := DirichletCharacter("626.1");
S:= CuspForms(chi, 4);
N := Newforms(S);
sage:from sage.modular.dirichlet import DirichletCharacter
H = DirichletGroup(626, base_ring=CyclotomicField(2))
chi = DirichletCharacter(H, H._module([0]))
N = Newforms(chi, 4, names="a")
This newform subspace can be constructed as the kernel of the linear operator
\( T_{3}^{25} - 16 T_{3}^{24} - 385 T_{3}^{23} + 7069 T_{3}^{22} + 58602 T_{3}^{21} - 1326517 T_{3}^{20} + \cdots + 80\!\cdots\!96 \)
acting on \(S_{4}^{\mathrm{new}}(\Gamma_0(626))\).