In our research, we studied 9 eclipsing binary stars: AR Lac, U CrB, S Equ, SU Boo, VV UMa, WW Gem, YY Eri, V0404 Lyr, HP Aur. We collected large sets of moments of minima from BRNO and observational data from AAVSO databases. Then, we obtained moments of minima for AAVSO observations (totally - 397 minima) using method of approximation with symmetric polynomial, realized in software MAVKA. This software was provided by Kateryna D. Andrych and Ivan L. Andronov. After we combined obtained moments of minima with data taken from BRNO and plotted O-C diagrams. For all stars these diagrams represented sinusoidal-like oscillations with superposition of linear (for SU Boo, VV UMa, WW Gem, V0404 Lyr and HP Aur) or parabolic trend (for AR Lac, U CrB, S Equ and YY Eri). The oscillations could be described as presence of the third component, which does not take part in eclipses, but causes the well-known light-time effect (LTE). That effect could be easily detected using long data series of observations. However, the oscillations have clear asymmetry, which we interpreted as elliptical shape of the third components orbit. Parabolic trend we explained as mass transfer between components of binary system. For all stars we computed minimal possible mass of the third component. In addition, we developed our own code in computed language Python and using it we computed orbital elements of the third component. Moreover, for stars with parabolic trend we calculated rate of the mass transfer. Finally, for all computed values we estimated errors.

We investigated two variable stars AH Tau and ZZ Cas. They are eclipsing binaries of W UMa and $\beta$ Lyr types. The period changes with time. The reason for steady changes of the period could be the mass transfer (the flow of matter) between components of these stellar systems. For ZZ Cas the changes of the period are cyclic caused by the presence of the hypothetical third component (either a small star or a large planet). The cyclic changes of the period for AH Tau superimpose on steady ones (the period decrease) indicating for a third component and mass transfer. The third components do not take part in the eclipses. However, due to their gravity, they make the visible close binary systems rotate and become closer or further to an observer, causing the Light-Time Effect (LTE). Generally, an orbit of a third component is not a circle, but an ellipse and it is inclined relatively to the observer's line of sight. Using the O-C curve, we estimated the parameters of a third component's orbit, such as a semi-major axis, an eccentricity, angles of orientation and a period of a third component's rotation. We elaborated a modeling program using the computer language Python. We calculated errors of third component's orbit parameters and even estimated its mass. The values of masses of the third components within errors of calculations show that the third components are probably stars. All these calculations were made using all available data from international database BRNO (Brno Regional Network for Observers). Moreover, we used moments of minima which we calculated as the result of observation processing from the AAVSO database.

In our research, we have studied 6 close binary stellar systems which are eclipsing variables of $\beta$ Lyrae and W Ursae Majoris types of variability. We have studied their O-C curves. To build them, we used moments of minima, listed in the BRNO database. Also, we used ones, which we obtained as the result of processing of observational data taken from AAVSO database. As the result of the O-C analysis, we detected that all of these stars have parabolic O-C curves, which might be caused by mass transfer from one component to another. In an addition, 3 from researched stars (U Peg, V0523 Cas and WZ Cyg) have superimposed parabolic and cyclic O-C changes that could be caused by presence of the third components in the systems. Also, we calculated minimal possible masses of third components, rates of the mass transfer for these systems and corrected the ephemerides.

The types of instability in the interacting binary stars are reviewed. The project "Inter-Longitude Astronomy" is a series of smaller projects on concrete stars or groups of stars. It has no special funds, and is supported from resources and grants of participating organizations, when informal working groups are created. Totally we studied 1900+ variable stars of different types. The characteristic timescale is from seconds to decades and (extrapolating) even more. The monitoring of the first star of our sample AM Her was initiated by Prof. V.P. Tsesevich (1907-1983). Since more than 358 ADS papers were published. Some highlights of our photometric and photo-polarimetric monitoring and mathematical modelling of interacting binary stars of different types are presented: classical, asynchronous, intermediate polars and magnetic dwarf novae (DO Dra) with 25 timescales corresponding to different physical mechanisms and their combinations (part "Polar"); negative and positive superhumpers in nova-like and many dwarf novae stars ("Superhumper"); eclipsing "non-magnetic" cataclysmic variables; symbiotic systems ("Symbiosis"); super-soft sources (SSS, QR And); spotted (and not spotted) eclipsing variables with (and without) evidence for a current mass transfer ("Eclipser") with a special emphasis on systems with a direct impact of the stream into the gainer star's atmosphere, or V361 Lyr-type stars. Other parts of the ILA project are "Stellar Bell" (interesting pulsating variables of different types and periods - M, SR, RV Tau, RR Lyr, Delta Sct) and "Novice"(="New Variable") discoveries and classification with a subsequent monitoring for searching and studying possible multiple components of variability. Special mathematical methods have been developed to create a set of complementary software for statistically optimal modelling of variable stars of different types.

A computer program is introduced, which allows to determine statistically optimal approxi-mation using the "Asymptotic Parabola" fit, or, in other words, the spline consisting of polynomials of order 1,2,1, or two lines ("asymptotes") connected with a parabola. The function itself and its derivative is continuous. There are 5 parameters: two points, where a line switches to a parabola and vice versa, the slopes of the line and the curvature of the parabola. Extreme cases are either the parabola without lines (i.e.the parabola of width of the whole interval), or lines without a parabola (zero width of the parabola), or "line+parabola" without a second line. Such an approximation is especially effective for pulsating variables, for which the slopes of the ascending and descending branches are generally different, so the maxima and minima have asymmetric shapes. The method was initially introduced by Marsakova and Andronov (1996OAP.....9..127M) and realized as a computer program written in QBasic under DOS. It was used for dozens of variable stars, particularly, for the catalogs of the individual characteristics of pulsations of the Mira (1998OAP....11...79M) and semi-regular (200OAP....13..116C) pulsating variables. For the eclipsing variables with nearly symmetric shapes of the minima, we use a "symmetric" version of the "Asymptotic parabola". Here we introduce a Windows-based program, which does not have DOS limitation for the memory (number of observations) and screen resolution. The program has an user-friendly interface and is illustrated by an application to the test signal and to the pulsating variable AC Her.

We used long-term visual amateur observations of several symbiotic variables for detection of periods that may be caused by pulsation. The examples of multiple periodicities are discussed individually in each case.

The studies of pulsating variable stars are traditional subjects of astronomers in Odessa. In the last half of the 20th century, the studies of the physical variable stars were the topics of 15 PhD theses of the collaborators of the Odessa State (now National) University. Continuing the tradition of studies of long-period variable stars in the Astronomical Observatory of Odessa University, we try to conduct a detailed classification of the Mira-type stars, semi-regular variable stars and also find out the location of symbiotic stars among other long-period variables in their evolution to the planetary nebulae, using the compiled material from own observations and that of other authors. For the research, we have used observations from the databases of the French Association of Variable Stars Observers (AFOEV) and Variable Stars Observers League of Japan (VSOLJ), which allow study of the variability of these stars during the time interval of about 100 years. Some stars were studied using the observations of the American Association of Variable Star Observers (AAVSO). Modern research of long-period variable stars in Odessa is conducted in several directions: 1. The analysis of correlations between the photometric parameters. 2. The study of changes on time of individual characteristics of the light curves. 3. The study of variability of the objects which are transitional between the Mira-type and the semi-regular variables. The conclusions are discussed.

The group of Mira-type and semi-regular variables with similar periodicity (multiperiodicty) is analyzed. They have periods of 230-260 days and 140-150 days and show intervals of periodical (Mira-type) variability with relatively high amplitude and "semi-regular" (SR-type) small-amplitude oscillations. Results of periodogram analysis are represented.