sage:from sage.modular.dirichlet import DirichletCharacter
H = DirichletGroup(2025, base_ring=CyclotomicField(12))
chi = DirichletCharacter(H, H._module([4, 3]))
N = Newforms(chi, 1, names="a")
magma://Please install CHIMP (https://github.com/edgarcosta/CHIMP) if you want to run this code
chi := DirichletCharacter("2025.757");
S:= CuspForms(chi, 1);
N := Newforms(S);
sage:traces = [4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
4,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]
f = next(g for g in N if [g.coefficient(i+1).trace() for i in range(46)] == traces)
gp:f = lf[1] \\ Warning: the index may be different