I want to calculate the square root of a polynomial with variable coefficients in sage. For example:
i,z=var('i,z') c=function('c') p=sum(c(i)*z^i,i,0,6) sqrt(p)
returns
sqrt(z^6c(6) + z^5c(5) + z^4c(4) + z^3c(3) + z^2c(2) + zc(1) + c(0))
I want a Taylor or Laurent expansion in z. How do I get this?