Simulation of 2D 1H homo- and 1H–13C heteronuclear NMR spectra of organic molecules by DFT calculations of spin–spin coupling constants and 1H and 13C-chemical shifts