The image of Abel-Jacobi embedding of a hyperelliptic Riemann surface of genus less than 5 to its Jacobian is the intersection of several shifted theta divisors with carefully chosen shifts.The necessity for computation of abelian integrals often arises in problems of classical mechanics [11,7], conformal mappings of polygons [4,12,17], solitonic dynamics [1], general relativity [14,8], rational optimization [3,5]-see also [2] and references in the mentioned sources. Function theory on Riemann surfaces allows one to evaluate with computer accuracy such integrals without any quadrature rules. As an example, let us consider Riemann's formula for the 39 abelian integral with two simple poles at the points R, Q of the curve X :where θ [·](·, ·) is Riemann theta function with some odd integer characteristics [ϵ, ϵ ′ ]. It might seem that the usage of this formula still requires the computation of holomorphic abelian integrals involved in the Abel-Jacobi (AJ) map:where g is the genus of the curve X and du s are suitably normalized abelian differentials of the first kind.In this note we show how to avoid the evaluation of AJ map: the image of low genus g < 5 hyperelliptic curve in its Jacobian is given as the solution of a (slightly overdetermined) set of equations containing theta functions. Moving along the curve embedded in its Jacobian we can compute the abelian integral by an explicit formula (e.g. (1)) and simultaneously compute a projection of the curve to the complex projective line (see e.g. [10,4,15]). In this way we get a parametric representation [4,12] of abelian integrals which may be used either for the evaluation of integral, or for its inversion. Alternative function theoretic approaches for the computation of abelian integrals use (hyperelliptic) σ and ℘ functions [8] or Schottky functions [4, Chapter 6].Computation of theta function still requires the knowledge of period matrix Π which should be retrieved from the known moduli of the curve. The latter problem is very important, however it is beyond the scope of this note, see e.g. [9] for the expression of period matrix in terms of theta constants.The techniques developed in this note were used e.g. for the computation of a conformal grid transferred from a rectangle-see Fig. 1 and [13]. The author thanks Oleg Grigoriev for the computations and the anonymous referee for numerous valuable references.