... We present a computer program for modelling energy flux distribution and light curves of eclipsing binaries from far ultraviolet to infrared regions. The Roche model is assumed. Proximity effects of components (reflection of radiation and surface gravity darkening) are taken into account. ... Figure 1: The surface temperature distribution for Algol A-B system before (upper plot) and after (lower plot) correction for temperature due to the irradiation effect. ...