a, q, k, n = var('a, q, k, n') show(sum(a*q^k, k, 0, n))
sage.version.version
%display latex
x, b, c = var('x b c') solve([x^2 + b*x + c == 0],x)