Search

The Prelle-Singer method allows determining an elementary first integral admitted by a polynomial vector field in the plane. It is a semi-algorithm whose nonlinear step consists of determining the Darboux polynomials of the vector field. In this article we construct a linear procedure to determine the Darboux polynomials present in the integrating factor of a polynomial vector field in the plane. Next, we extend the procedure to deal with rational 2ODEs that admit an elementary first integral

Here we present an efficient method for finding and using a nonlocal symmetry admitted by a rational second order ordinary differential equation (rational 2ODE) in order to find a Liouvillian first integral (belonging to a vast class of Liouvillian functions). In a first stage, we construct an algorithm (improving the methodde veloped in [1]) that computes a nonlocal symmetry of a rational 2ODE. In ase cond stage, based on the knowledge of this symmetry, it is possible to construct three polynomial vector fields (in R2), which "share" the Liouvillian first integral with the rational 2ODE. These "plane" polynomial vector fields can be used to construct a procedure (based on an idea developed in [2]) to determine an integrating factor for the rational 2ODE with a fast probabilistic algorithm. The main advantages of the proposed method are: the obtaining of the nonlocal symmetry is algorithmic and very efficient and, furthermore, its use to find an integrating factor is a sequence of linear or quasilinear processes.

We have already dealt with the problem of solving First Order Differential Equations (1ODEs) presenting elementary functions before in [1, 2]. In this present paper, we have established solid theoretical basis through a relation between the 1ODE we are dealing with and a rational second order ordinary differential equation, presenting a Liouvillian first Integral. Here, we have expanded the results in [3], where we have establish a theoretical background to deal with rational second order ordinary differential equations (2ODEs) via the S-function method. Using this generalisation and other results hereby introduced, we have produced a method to integrate the 1ODE under scrutiny. Our methods and algorithm are capable to deal efficiently with chaotic systems, determining regions of integrability.

Here we present a new approach to compute symmetries of rational second order ordinary differential equations (rational 2ODEs). This method can compute Lie symmetries (point symmetries, dynamical symmetries and non-local symmetries) algorithmically. The procedure is based on an idea arising from the formal equivalence between the total derivative operator and the vector field associated with the 2ODE over its solutions (Cartan vector field). Basically, from the formal representation of a Lie symmetry it is possible to extract information that allows to use this symmetry practically (in the 2ODE integration process) even in cases where the formal operation cannot be performed, i.e., in cases where the symmetry is nonlocal. Furthermore, when the 2ODE in question depends on parameters, the procedure allows an analysis that determines the regions of the parameter space in which the integrable cases are located.

Here we present a very efficient method to search for Liouvillian first integrals of second order rational ordinary differential equations (rational 2ODEs). This new algorithm can be seen as an improvement to the S-function method we have developed [24]. Here, we show how to further use the knowledge of the S-function to find an integrating factor of a set of first order rational ordinary differential equations (rational 1ODEs) which is shared by the original 2ODE, without having to actually solving these 1ODEs. This new use of the S-function, that is the theoretical basis of our new method to compute the integrating factor, proved to be a linear process of computation for a vast class of non-linear rational 2ODEs, making it much more efficient.

Here we present an efficient method to compute Darboux polynomials for polynomial vector fields in the plane. This approach is restricetd to polynomial vector fields presenting a Liouvillian first integral (or, equivalently, to rational first order differential equations (rational 1ODEs) presenting a Liouvillian general solution). The key to obtaining this method was to separate the procedure of solving the (nonlinear) algebraic systems resulting from the equation that translates the condition of existence of a Darboux polynomial into feasible steos (procedures that requires less memory consumption). We also present a brief performance analysis of the algorithms developed.

Dynamic systems have a fundamental relevance in the description of physical phenomena. The search for more accurate and faster numerical integration methods for the resolution of such systems is, therefore, an important topic of research. The present work introduces a new approach for the numerical integration of dynamic systems. We propose an association of numerical integration methods (integrators) in order to optimize the performance. The standard we apply is the balance of the duo : precision obtained x running time. The numerical integration methods we have chosen, for this particular instance of association, were the Runge-Kutta of fourth order and seventheighth order. The algorithm was implemented in C++ language. The results showed an improvement in accuracy over the lower grade numerical integrator (actually, we have achieved, basically, the precision of the top integrator) with a processing time performance closer to the one of the lower grade integrator. Similar results can be obtained for other pairs of numerical integration methods.

In [1,2] we have developed a method (we call it the S-function method) that is successful in treating certain classes of rational second order ordinary differential equations (rational 2ODEs) that are particularly `resistant' to canonical Lie methods and to Darbouxian approaches. In this present paper, we generalize the S-function method making it capable of dealing with a class of elementary 2ODEs presenting elementary functions. Then, we apply this method to a Duffing-Van der Pol forced oscillator, obtaining an entire class of first integrals.

Here we present a new approach to deal with first order ordinary differential equations (1ODEs), presenting functions. This method is an alternative to the one we have presented in [1]. In [2], we have establish the theoretical background to deal, in the extended Prelle-Singer approach context, with systems of 1ODEs. In this present paper, we will apply these results in order to produce a method that is more efficient in a great number of cases. Directly, the solving of 1ODEs is applicable to any problem presenting parameters to which the rate of change is related to the parameter itself. Apart from that, the solving of 1ODEs can be a part of larger mathematical processes vital to dealing with many problems., Comment: we have expanded the theoretical results presented on the piece and we have also introduced the presentation of a computer package allowing th euser to materialize the calculations in an easier way (the paper had 31 pages now it has 55)

In [1], we have presented the theoretical background for finding the Elementary Invariants for a 3D system of first order rational differential equations (1ODEs). We have also provided an algorithm to find such Invariants. Here we introduce new theoretical results that will lead to a novel, more efficient, approach to determine these invariants. Furthermore, one important aspect of such dynamical systems is that the integrability can be an issue. We will show that the present theoretical development allows for an Integrability analysis in the case where the system has free parameters., Comment: 10 pages

Here we present a new approach to search for first order invariants (first integrals) of rational second order ordinary differential equations. This method is an alternative to the Darbouxian and symmetry approaches. Our procedure can succeed in many cases where these two approaches fail. We also present here a Maple implementation of the theoretical results and methods, hereby introduced, in a computational package -- {\it InSyDE}. The package is designed, apart from materializing the algorithms presented, to provide a set of tools to allow the user to analyse the intermediary steps of the process., Comment: 42 pages

We have been working in many aspects of the problem of analyzing, understanding and solving ordinary differential equations (first and second order). As we have extensively mentioned, while working in the Darboux type methods, the most costly step of our methods and algorithms of solution is the determination of Darboux polynomials for the associated differential operators. Here, we are going to present some algorithms to greatly reduce the time expenditure in determining these needed Darboux polynomials. Some of them are based on a detailed analysis of the general structure of second order differential equations regarding the associated differential invariants. In order to perform this analysis, we produce a theorem concerning the general form for the differential invariants in terms of the Darboux polynomials.

In [Solving second order ordinary differential equations by extending the Prelle-Singer method, J. Phys. A: Math.Gen., 34, 3015-3024 (2001)] we defined a function (we called S) associated to a rational second order ordinary differential equation (rational 2ODE) that is linked to the search of an integrating factor. In this work we investigate the relation between these $S$-functions and the Lie symmetries of a rational 2ODE. Based on this relation we can construct a semi-algorithmic method to find the Lie symmetries of a 2ODE even in the case where it presents no Lie point symmetries.

The vast use of computers on scientific numerical computation makes the awareness of the limited precision that these machines are able to provide us an essential matter. A limited and insufficient precision allied to the truncation and rounding errors may induce the user to incorrect interpretation of his/hers answer. In this work, we have developed a computational package to minimize this kind of error by offering arbitrary precision numbers and calculation. This is very important in Physics where we can work with numbers too small and too big simultaneously., Comment: contribution to X hadron Physics - 2007 to appear on International Journal of Modern Physics E

Here we present/implement an algorithm to find Liouvillian first integrals of dynamical systems in the plane. In \cite{JCAM}, we have introduced the basis for the present implementation. The particular form of such systems allows reducing it to a single rational first order ordinary differential equation (rational first order ODE). We present a set of software routines in Maple 10 for solving rational first order ODEs. The package present commands permitting research incursions of some algebraic properties of the system that is being studied.

In this paper, we are concerned with improving the forecast capabilities of the Global approach to Time Series. We assume that the normal techniques of Global mapping are applied, the noise reduction is performed, etc. Then, using the mathematical foundations behind such approaches, we propose a method that, without a great computational cost, greatly increase the accuracy of the corresponding forecasting.

Here we present an algorithm to find elementary first integrals of rational second order ordinary differential equations (SOODEs). In \cite{PS2}, we have presented the first algorithmic way to deal with SOODEs, introducing the basis for the present work. In \cite{royal}, the authors used these results and developed a method to deal with SOODEs and a classification of those. Our present algorithm is based on a much more solid theoretical basis (many theorems are presented) and covers a much broader family of SOODEs than before since we do not work with restricted ansatz. Furthermore, our present approach allows for an easy integrability analysis of SOODEs and much faster actual calculations.

Here we present a method to find elementary first integrals of rational second order ordinary differential equations (SOODEs) based on a Darboux type procedure \cite{ManMac,firsTHEOps1,secondTHEOps1}. Apart from practical computational considerations, the method will be capable of telling us (up to a certain polynomial degree) if the SOODE has an elementary first integral and, in positive case, finds it via quadratures.

We present a semi-decision procedure to tackle first order differential equations, with Liouvillian functions in the solution (LFOODEs). As in the case of the Prelle-Singer procedure, this method is based on the knowledge of the integrating factor structure., Comment: 11 pages, latex

We present an algorithm to solve First Order Ordinary Differential Equations (FOODEs) extending the Prelle-Singer (PS) Method. The usual PS-approach miss many FOODEs presenting Liouvillian functions in the solution (LFOODEs). We point out why and propose an algorithm to solve a large class of these previously unsolved LFOODEs. Although our algorithm does not cover all the LFOODEs, it is an elegant extension mantaining the semi-decision nature of the usual PS-Method., Comment: 11 pages, latex

In math-ph/0107007, we present a method to tackle first order ordinary differential equations whose solutions contain Liouvillian functions (LFOODEs), many of them missed by the usual PS-approach. Here, we demonstrate an important result concerning the generality of that method., Comment: 8 pages, latex

A set of MapleV R5 software routines for solving first order ordinary differential equations (1ODEs) is presented. The package implements the Prelle-Singer Method in its original form plus its extension to include elementary functions (ELFs) on the integrating factor . The package also presents a theoretical extension to deal with all 1ODEs presenting liouvillian solutions (LIS)., Comment: 19 pages (no figures), latex

An extension of the ideas of the Prelle-Singer procedure to second order differential equations is proposed. As in the original PS procedure, this version of our method deals with differential equations of the form y''={M(x,y,y')}/{N(x,y,y')}, where M and N are polynomials with coefficients in the field of complex numbers C. The key to our approach is to focus not on the final solution but on the first-order invariants of the equation. Our method is an attempt to address algorithmically the solution of SOODEs whose first integrals are elementary functions of x, y and y'., Comment: 10 pages

The possibility of interaction between Maple and numeric compiled languages in performing extensive numeric calculations is exemplified by the Ndynamics package, a tool for studying the (chaotic) behavior of dynamical systems. Programming hints concerning the construction of Ndynamics are presented. The system command, together with the application of the black-box concept, is used to implement a powerful cooperation between Maple code and some other numeric language code., Comment: Latex2e, 9 pages, 4 figures. All the material (including the package, paper, figures) soon to be available on http://www.dft.if.uerj.br/symbcomp.htm

A set of MapleV R.4/5 software routines for calculating the numerical evolution of dynamical systems and flexibly plotting the results is presented. The package consists of an initial condition generator (on which the user can impose quite general constraints), a numerical solving manager, plotting commands that allow the user to locate and focus in on regions of possible interest and, finally, a set of routines that calculate the fractal dimension of the boundaries of those regions. A special feature of the software routines presented here is an optional interface in C, permitting fast numerical integration using standard Runge-Kutta methods, or variations, for high precision numerical integration, Comment: Latex, 19 pages, 7 figures. All the material (including the package, paper, figures) soon to be available on http://www.dft.if.uerj.br/symbcomp.htm

An update of the ODEtools Maple package, for the analytical solving of 1st and 2nd order ODEs using Lie group symmetry methods, is presented. The set of routines includes an ODE-solver and user-level commands realizing most of the relevant steps of the symmetry scheme. The package also includes commands for testing the returned results, and for classifying 1st and 2nd order ODEs., Comment: 24 pages, LaTeX, Soft-package (On-Line help) and sample MapleV sessions available at: http://dft.if.uerj.br/odetools.htm or http://lie.uwaterloo.ca/odetools.htm

A set of Maple V R.3/4 computer algebra routines for the analytical solving of 1st. order ODEs, using Lie group symmetry methods, is presented. The set of commands includes a 1st. order ODE-solver and routines for, among other things: the explicit determination of the coefficients of the infinitesimal symmetry generator; the construction of the most general invariant 1st. order ODE under given symmetries; the determination of the canonical coordinates of the underlying invariant group; and the testing of the returned results., Comment: 14 pages, LaTeX, submitted to Computer Physics Communications. Soft-package (On-Line Help) and sample MapleV session available at: http://dft.if.uerj.br/symbcomp.htm or ftp://dft.if.uerj.br/pdetools

We study a class of superstring models compactified in the 3-generation Calabi-Yau manifold of Tian and Yau. Our analysis includes the complete $E_6$-singlet sector, which has been recently evaluated using techniques of spectral and exact sequences. We use the discrete symmetries of the models to find flat directions of symmetry breaking that leave unbroken a low energy matter parity and make all leptoquarks heavy while preserving light Higgs fields. Then we classify the patterns of ordinary quark mass matrices and show that (without invoking effects due to nonrenormalizable terms) only one structure can accommodate the observed value of fermion masses and mixing angles, with preference for a heavy {\it top} quark ( $m_t\ge 170$ GeV for $V_{13}\le 0.013$ ). The model, which unifies perturbatively and predicts a realistic structure of quark mass matrices with texture zeroes, is one of the many possible string vacua. However, in contrast with what is often assumed in the search for realistic unified scenarios, it is highly nonminimal near the unification scale and the predicted mass matrices have no simple symmetry properties., Comment: 30 (including Tables and Figures), UG-FT-38/94