(2**19937-1)%2*F*52
solve(52*F, F)
diff(52*F, F)
series(52*F, F, 0, 10)
Want a full Python shell? Use SymPy Live.