for N in [11..5000] do 
   F := Factorization(N);
   if #F eq 2 and F[1][2] eq 1 and F[2][2] eq 1 then
      M := ModularSymbols(N);
      print N,"\t", RealTamagawaNumber(CS(M));
      DisownChildren(M);
   end if;
end for;
