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.

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.

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.

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

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.