(z(r)=((c(r^2))/(1+(Sqrt[1-(1+k)(c^2)(r^2)])))+(Subscript[a, 2])(r^2)+(Subscript[a, 4])(r^4)+(Subscript[a, 6])(r^6)+(Subscript[a, 8])(r^8)+…   .(2))