magma://Please install CHIMP (https://github.com/edgarcosta/CHIMP) if you want to run this code
chi := DirichletCharacter("3936.3361");
S:= CuspForms(chi, 2);
N := Newforms(S);
sage:from sage.modular.dirichlet import DirichletCharacter
H = DirichletGroup(3936, base_ring=CyclotomicField(2))
chi = DirichletCharacter(H, H._module([0, 0, 0, 1]))
N = Newforms(chi, 2, names="a")
sage:traces = [22,0,0,0,-4,0,0,0,-22,0,0,0,0,0,0,0,0,0,0,0,4,0,-8,0,30,0,0,
0,0,0,-16]
f = next(g for g in N if [g.coefficient(i+1).trace() for i in range(31)] == traces)
gp:f = lf[1] \\ Warning: the index may be different
This newform subspace can be constructed as the intersection of the kernels of the following linear operators acting on \(S_{2}^{\mathrm{new}}(3936, [\chi])\):