Abstract: | This paper is a continuation of work by Forest and Lee [1,2]. In [1,2] it was proved that the function theory of periodic soliton solutions occurs on the Riemann surfaces ? of genusN, where the integrals over paths on ? play the most fundamental role. In this paper a numerical method is developed to evaluate these integrals. Predisely, the aim is to develop a computational code for integrals of the form $$intlimits_gamma {f(z)frac{{dz}}{{R(z)}}, or} intlimits_gamma {f(z)R(z)dz,} $$ wheref(z) is any single-valued analytic function on the complex planeC, andR(z) is a two-valued function onC of the form $$R^2 (z) = prodlimits_{k = 1}^{2N + delta } {(z - z_0 (k)), delta = 0 or 1,} $$ where {z 0(k),1≤k≤2N+δ} are distinct complex numbers which play the role of the branch points of the Riemann surface ? = {(z, R(z))} of genusN?1+δ. The integral path γ is continuous on ?. The numerical code is developed in “Mathematica” [3]. |