On the non-linear spectroscopy including saturated absorption and four-wave mixing in two and multi-level atoms: a computational study