26 research outputs found
On the Radius of Convergence of the Simplest Power Series Solution of Painleve-I equation II(A Singular Solution)
In our previous paper [1] we considered the simplest power series solution of the Painleve-I equation which is regular at the origin. This note is a sequel to it. Here we consider another simplest Laurent series solution which is singular at the origin. Important feature of this solution is the location of the singularities. The location of the nearest singularity from the origin is given by the radius S of convergence of this Laurent series. The value of S is calculated numerically by the same method as in [1]. We obtained S = 2.56.... Various theoretical bounds for S are also obtained. The mathematical part of this work was done by Kametaka and the numerical part by Noda
A hybrid approach for the integration of a rational function
AbstractA hybrid integration algorithm obtaining an indefinite integral of a rational function (say q/r, q and r are polynomials) with floating-point but real coefficients is proposed. The algorithm consists of four steps and is based on combinations of symbolic and numeric computations (hybrid computation). The first step is a hybrid preprocessing stage. An integrand is decomposed into rational and logarithmic parts by using an approximate Horowitz' method which allows floating-point coefficients. Here, we replace the Euclidean GCD algorithm with an approximate-GCD algorithm which was proposed by Sasaki and Noda recently. It is easy to integrate the rational part. The logarithmic part is integrated numerically in the second step. Zeros of a denominator of it are computed by the numerical Durand-Kerner method which computes all zeros of a polynomial equation simultaneously. The integrand is then decomposed into partial fractions in the third step. Coefficients of partial fractions are determined by residue theory. Finally, in the fourth step, partial fractions are transformed into the resulting indefinite integral by using well-known rules of integrals. The hybrid algorithm proposed here gives both indefinite integrals and accurate values of definite integrals. Numerical errors in the hybrid algorithm depend only on errors in the second step. The algorithm evaluates some problems where numerical methods are inefficient or incapable, or a pure symbolic method is theoretically insufficient