Recent works (Steeghs & Casares 2002, Casares et al. 2003, Hynes et al. 2003) have demonstrated that Bowen flourescence is a very efficient tracer of the companion star in LMXBs. We present a numerical code to simulate time-delayed transfer functions in LMXBs, specific to the case of reprocessing in emission lines. The code is also able to obtain geometrical and binary parameters by fitting observed (X-ray + optical) light curves using simulated annealing methods. In this work we present the geometrical model for the companion star and the analytical model for the disc and show synthetic time-delay transfer functions for different orbital phases and system parameters.