A mathematical model for simulating th