We present a code which can be used to simulate the light curve of a binary system harbouring a non-circular accretion disk, and show that such a system geometry can produce light curves with asymmetric minima and different heights of the maxima. We have applied this code to check for the possible existence of a non-circular accretion disk, as suggested by Smak & Plavec (1997), in the Algol-type binary KU Cyg. We compare the best solutions for two sets of multicolour light curves, obtained within circular and non-circular disk models.