Investigation of the properties of MHD waves in the presence of non-uniform equilibria and flow using a unique numerical approach