Division by 2 on odd degree hyperelliptic curves and their jacobians

Let $K$ be an algebraically closed field of characteristic different from 2, $g$ a positive integer, $f(x)$ a degree $(2g+1)$ polynomial with coefficients in $K$ and without multiple roots, $C:y^2=f(x)$ the corresponding genus $g$ hyperelliptic curve over K, and $J$ the jacobian of $C$. We identify $C$ with the image of its canonical embedding into $J$ (the infinite point of $C$ goes to the identity element of $J$). It is well known that for each $mathfrak{b} in J(K)$ there are exactly $2^{2g}$ elements $mathfrak{a} in J(K)$ such that $2mathfrak{a}=mathfrak{b}$. M. Stoll constructed an algorithm that provides Mumford representations of all such $mathfrak{a}$, in terms of the Mumford representation of $mathfrak{b}$. The aim of this paper is to give explicit formulas for Mumford representations of all such $mathfrak{a}$, when $mathfrak{b}in J(K)$ is given by $P=(a,b) in C(K)subset J(K)$ in terms of coordinates $a,b$. We also prove that if $g>1$ then $C(K)$ does not contain torsion points with order between $3$ and $2g$.