Numerical approaches to the master equation of open quantum systems