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.

The algorithm of modified wavelet analysis is discussed. It is based on the weighted least squares approximation. Contrary to the Gaussian as a weight function, we propose to use a compact weight function. The accuracy estimates using the statistically correct expressions for the least squares approximations with an additional weight function are compared with that obtained using the bootstrap method.

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.

We present a review of the results of long-term photometric monitoring of selected magnetic cataclysmic binary systems, which belong to a class named "Intermediate polars". We found a spin period variability in the V2306 Cygni system. We confirm the strong negative superhump variations in the intermediate polar RX J2133.7+5107 and improved a characteristic time of white dwarf spin-up in this system. We have investigated the periodic modulation of the spin phases with the orbital phase in MU Camelopardalis. We can propose simple explanation as the influence of orbital sidebands in the periodic signal produced by intermediate polar.

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.

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.

Astronomical data are typically irregular in time, e.g. the space (HIPPARCOS/TYCHO, KEPLER, GAIA, WISE etc.) and ground-based CCD (NSVS, ASAS, CRTS, SuperWASP etc.) and photographic (Harvard, Sonneberg, Odessa etc.) photometrical surveys. This leads to cancellation of the conditions, which lead to the orthogonality of the basic functions, and thus the simplified methods give biased parameters of the approximations. We have elaborated a series of algorithms and programs for statistically correct analysis, and have applied them to 2000+ variable stars of different types. The data were obtained from an international collaboration in a frame of the "Inter-Longitude Astronomy" (ILA) campaign. Some highlights of our studies are presented, with an extended list of our original publications. The main improvements were done: 1) for the periodogram analysis - the parameters are determined from a complete set of equations containing the (algebraic polynomial) trend superimposed on the (multi-) harmonic wave, so no "detrending", no "prewhitening" are used; 2) for the approximations - we use additional (multi-) harmonic waves, and also "special shapes" (patterns) for parts of the light curve, which correspond to relatively fast changes (minima of the eclipsing binaries, minima and maxima for the pulsating variables); 3) "auto correlation analysis" (ACF) - taking into account the bias due to a trend removal (previously - only a subtraction of the sample mean was taken into account); ACF for the irregularly spaced data; 4) for the signals with bad coherence, the "scalegram" analysis is proposed, which allows to estimate a characteristic cycle length and the amplitude, as well as to provide a realistic approximation; 5) the extension of the Morlet-type wavelet for more periodic signals and 6) "running" (parabola, sine) approximations for aperiodic and "nearly periodic" variations, respectively.

We introduce the program MAVKA for determination of characteristics of extrema using observations in the adjacent data intervals, with intended applications to variable stars, but it may be used for signals of arbitrary nature. We have used a dozen of basic functions, some of them use the interval near extremum without splitting the interval (algebraic polynomial in general form, "Symmetrical" algebraic polynomial using only even degrees of time (phase) deviation from the position of symmetry argument), others split the interval into 2 subintervals (a Taylor series of the "New Algol Variable", "the function of Prof. Z. Mikul\'a\v{s}ek"), or even 3 parts ("Asymptotic Parabola", "Wall-Supported Parabola", "Wall-Supported Line", "Wall-Supported Asymptotic Parabola", "Parabolic Spline of defect 1"). The variety of methods allows to choose the "best" (statistically optimal) approximation for a given data sample. As the criterion, we use the accuracy of determination of the extremum. For all parameters, the statistical errors are determined. The methods are illustrated by applications to observations of pulsating and eclipsing variable stars, as well as to the exoplanet transits. They are used for the international campaigns "Inter-Longitude Astronomy", "Virtual Observatory" and "AstroInformatics". The program may be used for studies of individual objects, also using ground-based (NSVS, ASAS, WASP, CRTS et al.) and space (GAIA, KEPLER, HIPPARCOS/TYCHO, WISE et al.) surveys.

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.

The biography of Vladimir Platonovich Tsesevich (11.11.1907 - 28.10.1983), a leader of the astronomy in Odessa from 1944 to 1983, is briefly reviewed, as well as the directions of study, mainly the highlights of the research of variable stars carried out by the members of the scientific school founded by him. The directions of these studies cover a very wide range of variability types - "magnetic" and "non-magnetic" cataclysmic variables, symbiotic, X-Ray and other interacting binaries, classical eclipsers and "extreme direct impactors", pulsating variables from DSct and RR through C and RV to SR and M. Improved algorithms and programs have been elaborated for statistically optimal phenomenological and physical modeling. Initially these studies in Odessa were inspired by ("with a capital letter") Vladimir Platonovich Tsesevich. who was a meticulous Scientist and brilliant Educator, thorough Author and the intelligibly explaining Popularizer, persevering Organizer and cheerful Joker - a true Professor and Teacher. He was "the Poet of the Starry Heavens".

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.

For a group of the Mira-type stars, semi-regular variables and some RV Tau - type stars the limit cycles were computed and plotted using the phase plane diagrams. As generalized coordinates x and $\dot{x},$ we have used $\phi$ - the brightness of the star and its phase derivative. We have used mean phase light curves using observations of various authors from the databases of AAVSO, AFOEV, VSOLJ, ASAS and approximated using a trigonometric polynomial of statistically optimal degree. As generalized coordinates x and $\dot{x}$, we have used m - the brightness of the star and its phase derivative. For a simple sine-like light curve, the limit cycle is a simple ellipse. In a case of more complicated light curve, in which harmonics are statistically significant, the limit cycle has deviations from the ellipse. In an addition to a classical analysis, we use the error estimates of the smoothing function and its derivative to constrain an "error corridor" in the phase plane.

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.

Using the phase plane diagrams, the phase light curves of a group of the Mira-type stars and semi-regular variables are analyzed. As generalized coordinates x and $\dot{x}$, we have used m - the brightness of the star and its phase derivative. We have used mean phase light curves using observations of various authors. The data typically span a large time interval (nearly a century). They were compiled from the databases of AAVSO, AFOEV, VSOLJ, ASAS and approximated using a trigonometric polynomial of statistically optimal degree. As the resulting approximation characterizes the auto-oscillation process, which leads to a photometrical variability, the phase diagram corresponds to a limit cycle. For all stars studied, the limit cycles were computed. For a simple sine-like light curve, in e.g., L$_2$ Pup, the limit cycle is a simple ellipse. In a case of more complicated light curve, in which harmonics are statistically significant, the limit cycle has deviations from the ellipse. In an addition to a classical analysis, we use the error estimates of the smoothing function and its derivative to constrain an "error corridor" in the phase plane.

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.

The algorithm is described, which have been used for photo-polarimetric monitoring of various astronomical objects - cataclysmic variables, asteroids, comets.

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 discuss a method for determination of the size of the emitting region close to the compact star in a binary system with eclipses by a secondary, which fills its Roche lobe. The often used approach is to model the Roche lobe by a sphere with the "effective radius" corresponding to the volume of the Roche lobe. This approach leads to a 4-6% overestimate of the radius, if taking into account the angular dimensions of the Roche lobe seen form the compact star. Andronov (1992) had shown that the projection of the Roche lobe onto the celestial sphere is close to an ellipse and had tabulated these dimensions as a function of the mass ratio. Also he published the coefficients of the approximation similar to that of the Eggleton (1983) for the "sphere corresponding to the same volume". We compare results obtained for the "circle+circle", "ellipse+circle" and "ellipse+point" approximations of the projections of the red dwarf and a white dwarf, respectively. Results are applied to the recently discovered eclipsing polar CSS 081231:071126+440405.

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.

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.

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.

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.

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 dependencies of the phase of eclipse of the white dwarf's centre and the durations of the ascending and descending branches of the light curve on the binary system's parameters were computed using the spherically-symmetric approximation and the more accurate model of the elliptical projection onto the celestial sphere of the companion (red dwarf) that fills its Roche lobe. The parameters of eclipses in the classical eclipsing polar OTJ 071126+440405 = CSS 081231:071126+440405 were estimated.

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 discuss methods for modeling eclipsing binary stars using the "physical", "simplified" and "phenomenological" models. There are few realizations of the "physical" Wilson-Devinney (1971) code and its improvements, e.g. Binary Maker, Phoebe. A parameter search using the Monte-Carlo method was realized by Zola et al. (2010), which is efficient in expense of too many evaluations of the test function. We compare existing algorithms of minimization of multi-parametric functions and propose to use a "combined" algorithm, depending on if the Hessian matrix is positively determined. To study methods, a simply fast-computed function resembling the "complete" test function for the physical model. Also we adopt a simplified model of an eclipsing binary at a circular orbit assuming spherical components with an uniform brightness distribution. This model resembles more advanced models in a sense of correlated parameter estimates due to a similar topology of the test function. Such a model may be applied to detached Algol-type systems, where the tidal distortion of components is negligible.

A program WWZ is introduced, which realizes the wavelet analysis using an improved modification of the algorithm of the Morlet wavelet for a general case of irregularly spaced data, which is typical for the databases available in virtual observatories. Contrary to the well-known analogs, working with regularly spaced (equidistant in time) arguments, we have implemented an improved algorithm presented by Andronov, (1998KFNT...14..490A, 1999sss..conf...57A), which significantly increases the signal-to-noise ratio. The program has been used to study semi-regular pulsating variable stars (U Del et al.), but can be used for the analysis of signals of any nature.

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.

We made a comparative analysis of numerical methods for multidimensional optimization. The main parameter is a number of computations of the test function to reach necessary accuracy, as it is computationally "slow". For complex functions, analytic differentiation by many parameters can cause problems associated with a significant complication of the program and thus slowing its operation. For comparison, we used the methods: "brute force" (or minimization on a regular grid), Monte Carlo, steepest descent, conjugate gradients, Brent's method (golden section search), parabolic interpolation etc. The Monte-Carlo method was applied to the eclipsing binary system AM Leo.

Variability of the spin period of the white dwarf in the V405 Aur (RX J0558.0+5353) system using our observations and previously published maxima timings is analyzed. As the phase light curve contains two nearly equal photometric waves, one maximum was set as a "primary" one. The ephemeris for the maxima of the "spin variability" (due to rotation of the magnetized white dwarf) for recent seasons 2010-2012 is Tmax = HJD 2455882.470614(25)+0.00631314602(46)E. This corresponds to a significant negative trend at the "O-C" diagram. Due to significant gaps in the observational data and statistical error of timings, there may be some suggestions on the spin period variability - a fast period "jump" in 2007y; secular period variations; a cubic ephemeris (which may be interpreted by a precession of the magnetic white dwarf at a time-scale of decades) or a periodic change with a period of 6.2 years and semi-amplitude of 17.2\pm1.8 sec. For the present observations, more reliable are two latter models. To distinguish between them, a continuation of monitoring is needed. The periodic variations may be interpreted by a light time effect caused by a third body of mass (M3 > 0.09M\odot), which corresponds to a low-mass star, but not to an extra-solar planet. In this case, the system belongs to a rare class of cataclysmic variables with a third body.

We present the results of photometric CCD observations of the magnetic cataclysmic variable V405 Aurigae (RX J0558.0+5353 = 1RXS J055800.7+535358) obtained using different instruments. We analysed variability of the spin period of the white dwarf in the V405 Aur (RX J0558.0+5353) system using our observations and previously published maxima timings. The spin period of the system in 2010-2012 is $P=545.4558163(94)$s. As we have gaps in observational data, we present 2 hypotheses of the spin period variability of this system - a cubic ephemeris which may be interpreted by a precession of the magnetic white dwarf or a periodic change with a period of 6.2 years and semi-amplitude of $17.2\pm1.8$ sec. The periodic variations may be interpreted by a light-time effect caused by a low-mass star ($M_3\ge0.09M_\odot$). In this case, the system belongs to a rare class of cataclysmic variables with a third body.

The observations of the two-periodic magnetic cataclysmic system EX Hya have been carried out, using the telescopes RC16 and TOA-150 of the Tzec Maun observatory. 6 nights of observations were obtained in 2010-2011 (alternatively changing filters VR). Also the databases of WASP, ASAS and AAVSO have been analyzed. Processing time series was carried out using the program MCV. We analyzed changes in the rotation period of the white dwarf, and based on our own and previously published moments of maximum. The ephemeris was determined for the maxima of the radiation flux associated with the rotation of the magnetic white dwarf: Tmax=2437699.89079(59) +0.0465464808(69).E-6.3(2)*10^{-13}E^2, which corresponds to the characteristic timescale of the rotation spin-up of 4.67(14)*10^6 years. This contradicts the estimated value of the mass of the white dwarf of 0.42M_\odot, based on X-ray observations made by Yuasa et al (2010), however, is consistent with estimates of the masses of 0.79 M_\odot (white dwarf) and 0.108 M_\odot (red dwarf) previously published Beuermann and Reinsch (2008), and the assumption that the capture of accreted plasma by magnetic field of the white dwarf is near the border of the Roche lobe. Analyzed moments do not support the assumption of Mauche et al (2009) for a statistically significant cubic term in the ephemeris. Despite the presence of outbursts in EX Hya, there are significant differences from the DO Dra, which supports the introduction to a detailed classification of the intermediate polars the groups of "outbursting intermediate polars" and "magnetic dwarf novae."

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.

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).

The variation of average brightness during the time and the variation of the amplitude of the Mira-type stars T Cep, X CrB, U UMi, U Cyg, BG Cyg is studied. For the research, the observations of the members of the French association of observers of variable stars (AFOEV) covering almost 100 years are used. All stars show cyclic variations of the specified parameters. For T Cep, U UMi and U Cyg, the values of quasi-period of the variations of average brightness are about 1400, 1025 and 1680 days, respectively. For the stars T Cep, U UMi and BG Cyg, the period of the brightness variations changes, but we used a mean one. The dependences {\guillemotleft}average brightness - JD{\guillemotright} and {\guillemotleft}amplitude - JD{\guillemotright} are plotted. Results are discussed.

We present results of CCD photometric study of the intermediate polars FO Aqr, EX Hya and RXS J180340.0+401214 = USNO-A2.0 1275-09738647 and corresponding corrected ephemerids.

We discuss different methods for phenomenological approximations of signals with (generally) irregularly spaced arguments. Such signals may be classified as periodic, multi-periodic, quasi-periodic (cyclic), burst-type and flicker-type (see e.g. a review by Andronov (2003)). In this talk, we concentrate on (nearly) periodic light curves of the Algol-type eclipsing variables. Having no possibility to determine all physical parameters for all stars and to make physical modeling, the phenomenological fits are used. This is especially important for discovery papers, when there are often available only mono-chromatic photometric data, without estimates of temperature and orbital velocities of the components.

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.

During processing the observations of the intermediate polar 1RXS J180340.0+401214, obtained 26.05.2012 at the 60-cm telescope of the Mt. Suhora observatory (Krakow, Poland), variability of 2MASS J18024395+4003309 was discovered. As this object was not listed in the "General Catalogue of Variable Stars" or "Variable Stars Index", we registered it as VSX J180243.9+400331. Additionally we used 189 separate observations from the Catalina Sky Survey spread over 7 years. The periodogram analysis yields the period of 0d.3348837{\pm}0d.0000002.The object was classified as the Algol-type eclipsing binary with a strong effect of ellipticity. The depths of the primary and secondary minima are nearly identical, which corresponds to a brightness (and maybe) mass ratio close to 1. The statistically optimal degree of the trigonometric polynomial n=4. The most recent minimum occurred at HJD 2456074.4904. The brightness range from our data is 16.56-17.52 (V), 16.18-17.08 (R). The NAV ("New Algol Variable") algorithm was applied for statistically optimal phenomenological modeling and determination of corresponding parameters.

We introduce a special class of functions for mathematical modeling of periodic signals of special shape with irregularly spaced arguments. This method was developed for determination of phenomenological characteristics of the light curves, which are necessary for registration in the "General Catalogue of Variable Stars" and other databases. For eclipsing binary stars with smooth light curves - of types EB and EW - it is recommended a trigonometric polynomial of optimal degree in a complete or symmetric form. For eclipsing binary systems with narrow minima (EA-type), statistically optimal is an approximation in a class of non-polynomial spline functions. It is used a combination of the second-order trigonometric polynomial (TP2, what describes effects of "reflection", "ellipsoidality" and "spotness") and localized contributions of minima (parametrized in depth and profile separately for primary and secondary minima). Such an approach is characterized by a statistical accuracy of the smoothing curve, which is up to ~1.5-2 times better than the trigonometric polynomial of statistically optimal degree, and the absence of false "waves" in the light curve associated with the effect of Gibbs. In addition to the minimum width, which can not be determined by a trigonometric polynomial approximation, the method allows to determine with better accuracy its depth, as well as to separate the effects of the eclipse and out-of-eclipse parts. For multi-color observations, improving the accuracy of the smoothing of the curve in each filter will allow to obtain with better accuracy the curves of the color index variations. Effectivity of the proposed method increases with decreasing eclipse depth. The method called NAV ("New Algol Variable"), was applied to eclipsing binary systems VSX J022427.8-104034=USNO-B1.0 0793-0023471 and BM UMa. For VSX0224, an alternative model of "double period" is discussed.

Two recently discovered variable stars (CzeV134 = GSC 3682 0018 = USNO-A2.0 1425-1870026 and CzeV135 = GSC 3682 2051 = USNO-A2.0 1425-1825909 = V1094 Cas), which have been identified in the field of the W UMa variable star BS Cas, are studied in the present paper. The phase curves and finding charts for these stars are presented. The ephemeris and other photometric parameters were computed. The phenomenological features indicate that the first star (CzeV134) is probably a low-amplitude RRc Lyrae - type variable star with the period P = 0.419794\pm0.000029 d and the initial epoch T0 = HJD2453236.50412\pm0.00056. The amplitude and the shape of the light curve are variable possibly indicating the Blazhko phenomenon. The second star (CzeV135) was classified as an EW-type binary system of subtype A. However, a {\beta} Lyrae type may not be excluded, as various classification parameters lie in a range of overlapping values for both classes. The period P=0.51429090\pm0.00000012 d and the initial epoch T0 = HJD2454543.7920\pm0.0006. The O'Connell effect is clearly visible. There are slight changes of this effect, noticeable while comparing different seasons of observations. O-C diagrams for these stars were analyzed.

We present the results of two color VR observation of the Algol-type binary system WZ Crv (12h44m15.19s, -21d25m35.4s) which were obtained using the remotely controlled telescope TOA-150 of Tzec Maun Observatory. We determined the moments of individual minima, the orbital period and its derivative, the initial epoch, color indices V-R and temperature estimates of the components. Also we noticed that the phase curve is asymmetric: the second maximum is higher than the first one. It indicates that there is a spot in the photosphere of one of the stars in this system.

The results of three-color BVR photometric study of a recently discovered unusual binary system VSX J052807.9+725606 = USNO-B1.0 1629-0064825 are presented. This system is very similar to V361 Lyr, which was assumed to be unique. We confirmed a strong asymmetry of the phase curve and found the wavelength dependence of the amplitude. This is interpreted by a "direct impact" model with a "hot spot" in the atmosphere of the accreting component. Color temperatures are determined. Characteristics of the "hot spot" are estimated. We also calculated the new ephemeris for VSX J052807.9+725606.

We study the variability of the nova-like cataclysmic variable TT Ari, on time-scales of between minutes and months. The observations in the filter R were obtained at the 40-cm telescope of the Chungbuk National University (Korea). TT Ari was in the "negative superhump" state after its return from the "positive superhump" state, which lasted 8 years. The ephemeris for 12 best pronounced minima is $T_{min}=BJD 2453747.0700(47)+0.132322(53)E.$ where numbers in digits are errors in units of the last digit. The phases of minima may reach 0.2, indicating non-eclipse nature of these minima. The quasi-periodic oscillations (QPO) are present with a mean "period" of 21.6 min and mean semi-amplitude of 36 mmag. This value is consistent with the range 15-25 minutes reported for previous "negative superhump" states and does not support the hypothesis of secular decrease of the QPO period. Either the period, or the semi-amplitude show significant night-to-night variations. According to the position at the two-parameter diagrams, the interval of observations was splitted into 5 parts, showing different characteristics which are discussed in this paper. The system is an excellent laboratory to study processes resulting in variations at time-scales from seconds to decades and needs further monitoring at various states of activity.

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

Results of a CCD study of the variability of the cataclysmic variable BG CMi obtained at the Korean 1.8m telescope in 2002-2005 are presented. The "multi-comparison star" method had been applied for better accuracy estimates. The linear ephemeris based on 19 mean maxima for 2002--2005 is HJD 2453105.31448(6)+0.01057257716(198)(E-764707). The period differs from that predicted by the quadratic ephemeris by Pych et al. (1996) leading to a possible cycle miscount. The statistically optimal ephemeris is a fourth-order polynomial, as a quadratic or even a cubic ephemeris leads to unaceptably large residuals: Min.HJD=$ 2445020.28095(28)+0.0105729609(57)E -1.58(32)\cdot10^{-13}E^2-5.81(64)\cdot10^{-19}E^3+4.92(41)\cdot10^{-25}E^4.$ Thus the rate of the spin-up of the white dwarf is decreasing. An alternative explanation is that the spin-up has been stopped during recent years. The deviations between the amplutudes of the spin variability in V and R, as well as between phases are not statistically significant. However, the orbital light curves exhibit distinct difference; the corresponding color index shows a nearly sinusoidal shape with a maximum at orbital phase $\sim0.2.$ The variations of the amplitude of spin waves shows a short maximum at the phase of the orbital dip. The corrected ephemeris for orbital minima is Min.HJD=$2448368.7225(12)+0\fd13474841(6)\cdot(E-24849)$ with a narrow dip occuring 0.07P later. The rate of the spin period variation seems to be changed, justifying the necessity of regular observations of intermediate polars.

We report the discovery of a bright X-ray transient object, XMMU J174554.4-285456, observed in outburst with XMM-Newton on October 3, 2002,and located at 6.3' from SgrA*, the supermassive black hole at the Galactic center.This object exhibits a very large X-ray luminosity variability of a factor of about 1300 between two X-ray observations separated by four months. The X-ray spectrum is best fitted by a power-law with a photon index of 1.6+/-0.2 and absorption column density of 14.1 (+1.6,-1.4) x 10^22 cm^-2. This large absorption suggests this source is located at the distance of the Galactic center, i.e., 8 kpc. The 2-10 keV luminosity is about 1.0 x 10^35(d/8kpc)^2 erg/s. A pulsation period of about 172 s is hinted by the timing analysis. The X-ray properties strongly suggest a binary system with either a black hole or a neutron star for the compact object.

The processes of accretion of the gaseous blobs with different masses and densities onto strongly magnetized white dwarfs in the systems of polars have been modeled. We have proved that shot noise in blue wavelengths represents accretion of the smaller and denser blobs than in redder wavelengths. Using combined "smooth particle hydrodynamics - drag force" model, we have predicted a shape of the accretion stream and active regions on the white dwarf surface.

A study of AM Her low states in September 1990 and 1991 and June-July 1997 is reported from a coordinated campaign with observations obtained at the Haute-Provence observatory, at the 6-m telescope of the Special Astrophysical Observatory and at the 2.6m and 1.25m telescopes of the Crimean observatory. Spectra obtained at different dates when the source was in low states at a comparable V magnitude, show the presence of strong Zeeman absorption features and marked changes in emission lines with a day-to-day reappearance of the HeII (4686\AA) emission lines in 1991. Despite this variability, the magnetic field inferred from the fitting of the absorption spectrum with Zeeman hydrogen splitting, is remarkably constant with a best value of (12.5$\pm$0.5)MG. Detailed analysis of the UBVRI light curves shows the presence of repetitive moderate amplitude ($\sim$ 0.3-0.5 mag) flares predominantly red in colour. These flares are attributed to small accretion events and are compared to the large ($\sim$ 2 mag.) blue flare reported by Shakhovskoy et al. (1993). We suggest that the general flaring activity observed during the low states is generated by accretion events. The different characteristics of the flares (colour and polarization) are the results of different shock geometries depending on the net mass accretion flux.