for N in [11..5000] do 
   M := ModularSymbols(N);
   print N,"\t", RealTamagawaNumber(CS(M));
   DisownChildren(M);
end for;