A systematic path to non-Markovian dynamics: new response probability density function evolution equations under Gaussian coloured noise excitation