We present new analysis of O-C diagrams variations of three Algol-type eclipsing binary stars AD And, TW Cas and IV Cas. We have used all published minima times (including visual and photographic) as well as new determined ones from our and SuperWasp observations. We determined orbital parameters of the 3rd bodies in the systems with statistically significant errors, using our code based on genetic algorithms and Markov Chain Monte Carlo simulation. We confirmed multiple nature of AD And and triple-star model of TW Cas and we proposed quadruple-star model of IV Cas.

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 present our observations of the low-mass eclipsing binary GJ~3236. We have analyzed a phased $R_C$ light-curve and confirmed previously determined fundamental parameters of the components. We detected evolution of the spot(s) and found that there exists a large spot near a polar region of the primary component and another spot either on the primary or the secondary component. We also observed 7 flare events and determined a flare rate of about 0.1 flares per hour. We observed two high energy, long-term flares with a complex light curve and possibly four weak short-term flaring events. A majority of the flares was detected in the $R_C$ filter, which indicate their high energy.

We present new observations of the transiting system TrES-3 obtained from 2009 to 2011 at several observatories. The orbital parameters of the system were re- determined and a new linear ephemeris was calculated. The best quality light curve was used for light curve analysis, an other datasets were used to determine mid-transit times, Tc, and study transit time variation (TTV). For planet parameter determination we used two independent codes and finally, we concluded that our parameters are in agreement with previous studies. Based on our observations, we determined 14 mid-transit times. Together with published Tc we found that the timing residuals showed no significant deviation from the linear ephemeris. We concluded that a periodic TTV signal with an amplitude greater than 1 minute over a 4-year time span seems to be unlikely. Our analysis of an upper mass limit allows us to exclude an additional Earth-mass planet close to inner 3:1, 2:1, and 5:3 and outer 3:5, 1:2, and 1:3 mean-motion resonances. Using the long-term integration and applying the method of maximum eccentricity, the region from 0.015 au to 0.05 au was found unstable and the region beyond the 0.05 au was found to have a chaotic behaviour and its depletion increases with increasing values of the initial eccentricity as well as inclination.

We present basic informations about package SPHOTOM for an automatic multicolour photometry. This package is in development for a creation of photometric pipe-line, which we plan to use in near future with our new instruments. It could operate in two independent modes, (i) GUI mode, in which user can select images and control functions of package through interface and (ii) command line mode, in which all processes are controlled using a main parameter file. SPHOTOM is developed as an universal package for Linux based systems with easy implementation for different observatories. Photometric part of package is based on Sextrator code, what allow us to detect all objects on the images and perform their photometry with different apertures. We can also perform astrometric solution for all images for a correct cross-identification of the stars on the images. The result is a catalogue of all objects with their instrumental photometric measurements which are consequently used for a differential magnitudes calculations with one ore more comparison stars, transformations to international system and colour indices determinations.

Tightest known quadruple systems VW LMi consists of contact eclipsing binary with P_12 = 0.477551 days and detached binary with P_34 = 7.93063 days revolving in rather tight, 355.0-days orbit. This paper presents new photometric and spectroscopic observations yielding 69 times of minima and 36 disentangled radial velocities for the component stars. All available radial velocities and minima times are combined to better characterize the orbits and to derive absolute parameters of components. The total mass of the quadruple system was estimated at 4.56 M_sun. The detached, non-eclipsing binary with orbital period P = 7.93 days is found to show apsidal motion with U approximately 80 years. Precession period in this binary, caused by the gravitational perturbation of the contact binary, is estimated to be about 120 years. The wide mutual orbit and orbit of the non-eclipsing pair are found to be close to coplanarity, preventing any changes of the inclination angle of the non-eclipsing orbit and excluding occurrence of the second system of eclipses in future. Possibilities of astrometric solution and direct resolving of the wide, mutual orbit are discussed. Nearby star, HD95606, was found to form loose binary with quadruple system VW LMi.

Photometric analysis of $BVR_C$ light curves of newly discovered eclipsing binary GSC 0008-00901 is presented. The orbital period is improved to 0.28948(11) days. Photometric parameters are determined, as well. The analysis yielded to conclusion that system is an over-contact binary of W UMa type with components not in thermal contact. The light curves from 2005 show the presence of a spot on the surface of one of the components, while light curves from 2006 are not affected by maculation.

The ~15.1 years period found in the long-term UBV photoelectric and photographic photometry of the symbiotic nova V1016 Cyg is detected also in the (J-K) colour index and in the UV continuum and emission line fluxes from IUE and HUT spectra. It could be interpreted either as the effect of recurrent enhanced mass loss episodes from the Mira type variable companion to a hot component along its ultra-wide orbit (proposed from recent HST observations) or the true orbital period of the inner, unresolved binary of a triple system. A 410-day delay of the maximum of UV emission lines fluxes with respect to the maximum of continuum was found. The pulsation period of the Mira type variable was improved to 474+/-6 days.