I am new to sage and programming in general. I want to compute the power series representation of the square root of a polynomial with variable coefficients. Here is my attempt:
i,z=var('i,z') c=function('c') p=sum(c(i)*z^i,i,0,6) sqrt(p)
This returns
sqrt(z^6c(6) + z^5c(5) + z^4c(4) + z^3c(3) + z^2c(2) + zc(1) + c(0)).
How can I get sage to return a power series expression?