var('z')
y=integrate(1/(1+z*z))
p=taylor (y,z,0,6)
show(p)