We report an analysis of two poorly studied systems GSC 04396-00605 and GSC 04395-00485, which were recently named as V0455 Dra and V0454 Dra, respectively. For two eclipsing stars, the periods and epochs were significantly corrected using our extensive data. The phenomenological characteristics of the mean light curves were determined using the "New Algol Variable" (NAV) algorithm. The individual times of Maxima/Minima (ToM) were determined using the newly developed software MAVKA, which outputs the accurate parameters using the "asymptotic parabola" approximations. The light curves were approximated using the Wilson-Devinney model, the best fit parameters are listed. For both systems, the presence of the dark spot is justified, the parameters are estimated. Both systems classified to be overcontact binaries of EW type. The observations were obtained using the 1-m Korean telescope of the LOAO while monitoring the cataclysmic variable DO Dra within the "Inter-Longitude Astronomy" (ILA) project.

Advanced MAVKA software for the approximation of extrema observations is used to analyze the variability of the brightness of pulsating and eclipsing stars, but may be useful in analyzing signals of any nature. A new algorithm using a parabolic (quadratic) spline is proposed. In contrast to the traditional definition of a spline as a piecewise-defined function at fixed intervals, a spline is proposed to be divided into three intervals, but the positions of the boundaries between the intervals are additional parameters. The spline defect is 1, that is, the function and its first derivative are continuous and the second derivative can be discontinuous at the boundaries. Such a function is an enhancement of the "asymptotic parabola" (Marsakova and Andronov 1996). The dependence of the fixed signal approximation accuracy on the location of the boundaries of the interval is considered. The parameter accuracy estimates using the least squares method and bootstrap are compared. The variability of the semi-regular pulsating star Z UMa is analyzed. The presence of multicomponent variability of an object, including, four periodic oscillations and significant variability of the amplitudes and phases of individual oscillations is shown.

Multiple algorithms of time series analysis are briefly reviewed and partially illustrated by application to the visual observations of the semi-regular variable DY Per from the AFOEV database. These algorithms were implemented in the software MCV (Andronov and Baklanov, 2004), MAVKA (Andrych and Andronov, 2019; Andrych et al., 2019). Contrary to the methods of physical modeling, which need to use too many parameters, many of which may not be determined from pure photometry (like temperature/spectral class, radial velocities, mass ratio), phenomenological algorithms use smaller number of parameters. Beyond the classical algebraic polynomials, in the software MAVKA are implemented other algorithms, totally 21 approximations from 11 classes. Photometric observations of DY Per from the AFOEV international database were analyzed. The photometric period has switched from P=851.1d(4.1) to P=780.5d(2.7) after JD 2454187(9)d. A parameter of sinusoidality is introduced, which is equal to the ratio of effective semi-amplitudes of the signal determined from a sine fit and the running parabola scalegram.

The expert system for time series analysis of irregularly spaced signals is reviewed. It consists of a number of complementary algorithms and programs, which may be effective for different types of variability. Obviously, for a pure sine signal, all the methods should produce the same results. However, for irregularly spaced signals with a complicated structure, e.g. a sum of different components, different methods may produce significantly different results. The basic approach is based on classical method of the least squares (1994OAP.....7...49A). However, contrary to common "step-by-step" methods of removal important components (e.g. mean, trend ("detrending"), sine wave ("prewhitening"), where covariations between different components are ignored, i.e. erroneously assumed to be zero, we use complete mathematical models. Some of the methods are illustrated on the observations of the semi-regular pulsating variable RY UMa. The star shows a drastic cyclic change of semi-amplitude of pulsations between 0.01 to 0.37mag, which is interpreted as a bias between the waves with close periods and a beat period of 4000d (11yr). The dominating period has changed from 307.35(8)d before 1993 to 285.26(6)d after 1993. The initial epoch of the maximum brightness for the recent interval is 2454008.8(5). It is suggested that the apparent period switch is due to variability of amplitudes of these two waves and an occasional swap of the dominating wave.

The methods for determination of the characteristics of the extrema are discussed with an application to irregularly spaced data, which are characteristic for photometrical observations of variable stars. We introduce new special functions, which were named as the "Wall-Supported Polynomial" (WSP) of different orders. It is a parabola (WSP), constant line (WSL) or an "asymptotic" parabola (WSAP) with "walls" corresponding to more inclined descending and ascending branches of the light curve. As the interval is split generally into 3 parts, the approximations may be classified as a "non-polynomial spline". These approximations extend a parabolic/linear fit by adding the "walls" with a shape, which asymptotically corresponds to the brightness variations near phases of the inner contact. The fits are compared to that proposed by Andronov (2010, 2012) and Mikulasek (2015) and modified for the case of data near the bottom of eclipses instead of wider intervals of the light curve. The WSL method is preferred for total eclipses showing a brightness standstill. The WSP and WSAP may be generally recommended in a case of transit eclipses, especially by exoplanets. Other two methods, as well as the symmetrical polynomials of statistically optimal order, may be recommended in a general case of non-total eclipses. The method was illustrated by application to observations of a newly discovered eclipsing binary GSC 3692-00624 = 2MASS J01560160+5744488, for which the WSL method provides 12 times better accuracy.

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.

The phenomenological parameters of eclipsing binary stars, which are the prototypes of the EA, EB and EW systems are determined using the expert complex of computer programs, which realizes the NAV ("New Algol Variable") algorithm (Andronov 2010, 2012) and its possible modifications are discussed, as well as constrains for estimates of some physical parameters of the systems in a case of photometric observations only, such as the degree of eclipse, ratio of the mean surface brightnesses of the components. The half-duration of the eclipse is 0.0617(7), 0.1092(18) and 0.1015(7) for Algol, $\beta$ Lyrae and W UMa, respectively. The brightness ratio is 6.8$\pm$1.0, 4.9$\pm$1.0 and 1.15$\pm$0.13. These results show that the eclipses have distinct begin and end not only in EA (as generally assumed), but also in EB and EW - type systems as well. The algorithm may be applied to classification and study of the newly discovered (or poorly studied) eclipsing variables based on own observations or that obtained using photometric surveys.

A comparative analysis of the special shapes (patterns, profiles) of the eclipses applied for the phenomenological modeling of the light curves of eclipsing binary stars is conducted. Families of functions are considered, generalizing local approximations (Andronov, 2010, 2012) and the functions theoretically unlimited in a width, based on a Gaussian (Mikulasek, 2015). For an analysis, the light curve of the star V0882 Car = 2MASS J11080308 - 6145589 of the classic Algol - subtype (\beta Persei) is used. By analyzing dozens of modified functions with additional parameters, it was chosen the 14 best ones according to the criterion of the least sum of squares of deviations. The best are the functions with an additional parameter, describing profiles, which are limited in phase.

We review the method NAV (New Algol Variable) first introduced in 2012Ap.....55..536A, which uses the locally-dependent shapes of eclipses in an addition to the trigonometric polynomial of the second order (which typically describes the "out-of-eclipse" part of the light curve with effects of reflection, ellipticity and O'Connell). Eclipsing binary stars are believed to show distinct eclipses only if belonging to the EA type. With a decreasing eclipse width, the statistically optimal value of the trigonometric polynomial s (2003ASPC..292..391A) drastically increases from ~2 for elliptic (EL) variables without eclipses, ~6-8 for EW and up to ~30-50 for some EA with narrow eclipses. In this case of large number of parameters, the smoothing curve becomes very noisy and apparent waves (the Gibbs phenomenon) may be seen. The NAV set of the parameters may be used for classification in the GCVS, VSX and similar catalogs. The maximal number of parameters is m=12, which corresponds to s=5, if correcting both the period and the initial epoch. We have applied the method to few stars, also in a case of multi-color photometry (2015JASS...32..127A), when it is possible to use the phenomenological parameters from the NAV fit to estimate physical parameters using statistical dependencies. We conclude that the NAV approximation is better than the TP one even for the case of EW-type stars with much wider eclipses. It may also be used to determine timings (see 2005ASPC..335...37A for a review of methods) or to determine parameters in the case of variable period, using a complete light curve modeling the phase variations. The method is illustrated on 2MASS J11080447-6143290 (EA-type), USNO-B1.0 1265-0306001 and USNO-B1.0 1266-0313413 (EW-type) and compared to various other methods from the literature.

We report on analysis of the two-color VR CCD observations of the newly discovered variable 2MASS J18024395+4003309=VSX J180243.9+400331 obtained using the 1-m telescope of the Mt. Lemmon Observatory (LOAO) in the field of the intermediate polar V1323 Her. The extended version of this conference talk we published in 2015JASS...32..127A. The variability was reported in 2012OAP....25..150A, and the object was monitored. The two-color observations covered all phase interval. The object is classified as an Algol-type variable with tidally distorted components, and shows an asymmetry of the maxima (the O\'Connell effect). For phenomenological modeling, we used the trigonometric polynomial approximation of statistically optimal degree, and a recent method "NAV" (New Algol Variable) using local specific shapes for the eclipse. Methodological aspects are described, especially for the case of few color observations. Estimates of the physical parameters based on analysis of phenomenological parameters, are presented.

Phenomenological characteristics of the sample of the Algol-type stars are revised using a recently developed NAV ("New Algol Variable") algorithm (2012Ap.....55..536A, 2012arXiv 1212.6707A) and compared to that obtained using common methods of Trigonometric Polynomial Fit (TP) or local Algebraic Polynomial (A) fit of a fixed or (alternately) statistically optimal degree (1994OAP.....7...49A, 2003ASPC..292..391A). The computer program NAV is introduced, which allows to determine the best fit with 7 "linear" and 5 "non-linear" parameters and their error estimates. The number of parameters is much smaller than for the TP fit (typically 20-40, depending on the width of the eclipse, and is much smaller (5-20) for the W UMa and beta Lyrae - type stars. This causes more smooth approximation taking into account the reflection and ellipsoidal effects (TP2) and generally different shapes of the primary and secondary eclipses. An application of the method to two-color CCD photometry to the recently discovered eclipsing variable 2MASS J18024395 + 4003309 = VSX J180243.9 +400331 (2015JASS...32..101A) allowed to make estimates of the physical parameters of the binary system based on the phenomenological parameters of the light curve. The phenomenological parameters of the light curves were determined for the sample of newly discovered EA and EW - type stars (VSX J223429.3+552903, VSX J223421.4+553013, VSX J223416.2+553424, US-NO-B1.0 1347-0483658, UCAC3-191-085589, VSX J180755.6+074711= UCAC3 196-166827). Despite we have used original observations published by the discoverers, the accuracy estimates of the period using the NAV method are typically better than the original ones.

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.

Phenomenological modeling of variable stars allows determination of a set of the parameters, which are needed for classification in the "General Catalogue of Variable Stars" and similar catalogs. We apply a recent method NAV ("New Algol Variable") to eclipsing binary stars of different types. Although all periodic functions may be represented as Fourier series with an infinite number of coefficients, this is impossible for a finite number of the observations. Thus one may use a restricted Fourier series, i.e. a trigonometric polynomial (TP) of order s either for fitting the light curve, or to make a periodogram analysis. However, the number of parameters needed drastically increases with decreasing width of minimum. In the NAV algorithm, the special shape of minimum is used, so the number of parameters is limited to 10 (if the period and initial epoch are fixed) or 12 (not fixed). We illustrate the NAV method by application to a recently discovered Algol-type eclipsing variable 2MASS J11080308-6145589 (in the field of previously known variable star RS Car) and compare results to that obtained using the TP fits. For this system, the statistically optimal number of parameters is 44, but the fit is still worse than that of the NAV fit. Application to the system GSC 3692-00624 argues that the NAV fit is better than the TP one even for the case of EW-type stars with much wider eclipses. Model parameters are listed.

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.

We review one of complementary methods for time series analysis - the method of "Running Sines". "Crash tests" of the method include signals with a large period variation and with a large trend. The method is most effective for "nearly periodic" signals, which exhibit "wavy shape" with a "cycle length" varying within few dozen per cent (i.e. oscillations of low coherence). This is a typical case for brightness variations of long-period pulsating variables and resembles QPO (Quasi-Periodic Oscillations) and TPO (Transient Periodic Oscillations) in interacting binary stars - cataclysmic variables, symbiotic variables, low-mass X-Ray binaries etc. General theory of "running approximations" was described by Andronov (1997A &AS..125..207A), one of realizations of which is the method of "running sines". The method is related to Morlet-type wavelet analysis improved for irregularly spaced data (Andronov, 1998KFNT...14..490A, 1999sss..conf...57A), as well as to a classical "running mean" (="moving average"). The method is illustrated by an application to a model signal with strongly variable period, as well as to a semi-regular variable AF Cyg. Some other stars studied with this method are discussed, e.g. RU And (switching between "Mira-type" large amplitude oscillations and time intervals of "constancy"), intermediate polars MU Cam (1RXS J062518.2+733433) and BG CMi, magnetic dwarf nova DO Dra, symbiotic stars UV Aur and V1329 Cyg.

We review methods for determination of "quasi-periods" (or "cycle length") of signals of low coherence. Such type of variability was called "cyclic" for semi-regular red variables, or "quasi-periodic oscillations" (QPO) for fast variability in cataclysmic variables and related objects. Methods are illustrated by application to AF Cygni.

Photometric analysis of photometric variability of the semi-regular pulsating variable U Del is analyzed. From the international AFOEV database, 6231 brightness values in the time interval JD 2451602-55378 were chosen. For the periodogram analysis, we have used a trigonometric polynomial fit. Using the criterion of minimal variance of the approximation at arguments of observations, the optimal degree is s=1. Initial epoch for maximum is T0=JD2453340(3), the period P=1198(4)d. Mean brightness at maximum is 6.624(5), at minimum 7.124(5), i.e. the amplitude is 0.499(5)m. Besides this slow variability, there is a faster oscillation of a period of 119.45(6)d, amplitude 0.303(5)m and an initial epoch for maximum 2453215.1(5). These results are mean during the time interval after that analyzed in the catalogue of Chinarova and Andronov (2000). Also the method of "running sines" with a filter half-width Delta t=0.5P. The local mean (averaged over a short period) brightness varies in a range 6.58-7.41m, the semi-amplitude exhibits very strong variations from 0.01 to 0.46m. The phase is also variable - typically of a full amplitude of 0.5. Close to JD 2452589, occurred a phase jump by a complete period during a descending branch of a slow wave. This effect was not observed during other cycles. No significant correlation between mean brightness and amplitude of short-period oscillation was found. Despite significant variability of amplitude, the periodic contribution is not statistically significant. Also characteristics of individual brightness extrema were found. The study was made in a course of the international project "Inter-Longitude Astronomy" (Andronov et al., 2010) and national project "Ukrainian Virtual Observatory (Vavilova et al., 2012).

We study the variability of the cataclysmic variable DO Dra, on time-scales of between minutes and decades. The characteristic decay time dt/dm=0.902(3) days/mag was estimated from our 3 nights of CCD R observations. The quiescent data show a photometric wave with a cycle about 303(15)d. We analyzed the profile of the composite (or mean) outburst. We discovered however, that a variety of different outburst heights and durations had occurred, contrary to theoretical predictions. With increasing maximum brightness, we find that the decay time also increases; this is in contrast to the model predictions, which indicate that outbursts should have a constant shape. This is interpreted as representing the presence of outburst-to-outburst variability of the magnetospheric radius. A presence of a number of missed weak narrow outbursts is predicted from this statistical relationship. A new type of variability is detected, during 3 subsequent nights in 2007: periodic (during one nightly run) oscillations with rapidly-decreasing frequency from 86 to 47 cycles/day and a semi-amplitude increasing from 0.06 to 0. 10, during a monotonic brightness increase from 14. 27 to 14. 13. This phenomenon was observed only during an unusually prolonged event of about 1 mag brightening in 2007 (lasting till autumn), during which no (expected) outburst was detected. We refer to this behaviour as to the transient periodic oscillations (TPO). To study this new and interesting phenomenon, new regular photometric and spectral (in a target of opportunity mode) observations are required.

Results of a multi-color study of the variability of the magnetic cataclysmic variable BY Cam are presented. The observations were obtained at the Korean 1.8m and Ukrainian 2.6m, 1.2m and 38-cm telescopes in 2003-2005, 56 observational runs cover 189 hours. The variations of the mean brightness in different colors are correlated with a slope dR/dV=1.29(4), where the number in brackets denotes the error estimates in the last digits. For individual runs, this slope is much smaller ranging from 0.98(3) to 1.24(3), with a mean value of 1.11(1). Near the maximum, the slope becomes smaller for some nights, indicating more blue spectral energy distribution, whereas the night-to-night variability has an infrared character. For the simultaneous UBVRI photometry, the slopes increase with wavelength from dU/dR=0.23(1) to dI/dR=1.18(1). Such wavelength dependence is opposite to that observed in non-magnetic cataclysmic variables, in an agreement to the model of cyclotron emission. The principal component analysis shows two (with a third at the limit of detection) components of variablitity with different spectral energy distribution, which possibly correspond to different regions of emission. The scalegram analysis shows a highest peak corresponding to the 200-min spin variability, its quarter and to the 30-min and 8-min QPOs. The amplitudes of all these components are dependent on wavelength and luminosity state. The light curves were fitted by a statistically optimal trigonometrical polynomial (up to 4-th order) to take into account a 4-hump structure. The dependences of these parameters on the phase of the beat period and on mean brightness are discussed. The amplitude of spin variations increases with an increasing wavelength and with decreasing brightness