On numerical study of periodic solutions of a delay equation in biological models