Yes missed parentheses from my part.
The context is easy :
Let be the parametrization of a unit sphere ##\vec{r}(\theta,\phi)=(cos\theta cos\phi,cos\theta sin\phi,sin\theta),\theta\in [-\pi/2,\pi/2],\phi\in [0,2\pi] ## and the associated basis vectors ##\vec{e}_\theta=\frac{\partial\vec{r}}{\partial\theta}## and ##\vec{e}_\phi##
The metric is ##g=diag(1,cos(\theta)^2)##
I wanted to orthonormalize the metric via the associated quadratic form : ##\underbrace{(du\vec{e}_\theta+dv\vec{e}_\phi)}_{\vec{x}}g\vec{x}^T##
The ##du,dv## are infinitesimals on the tangent plane and I express them as a function of ##\theta,\phi##: ##du=tan(d\theta)##
##dv'=cos(\theta)^2tan(d\phi)##
Instead of linearizing I did what I wrote in the OP with successive integrations (##\int(d\theta)^n=\theta^n/n!## by fixing all the successive constants of integration to 0) and found
##u=\sum_{n=0}^\infty(-1)^n2^{2n+2}(2^{2n+2}-1)B_{2n+2}\frac{\theta^{2n+1}}{(2n+2)!(2n+1)!}##
I wanted to compute the shape of u,v which if you want should be a map of the sphere.
Is it periodic and what's the asymptotic behaviour ?
Could anyone help me or give me references in webpages, it should be found in common textbooks about differential geometry or map projections I suppose.
What I found weird is that if I consider a path on the sphere ##\theta(t),\phi(t)## then the endpoint coordinate on the map could depend on the path itself ?!

Thanks.