A method for computing accurately the epoch of minimum of an eclipsing variable