Consider the Heston model expressed as \begin{align} dS_t &= \mu S_t dt + S_t \sqrt{V_t} \big(\rho dW_t^{(1)}+\sqrt{1-\rho^2}dW_t^{(2)} \big); \tag*{(1)} \\ dV_t &= \kappa(\theta - V_t)dt + \sigma \sqrt{V_t}dW_t^{(1)}, \tag*{(2)} \end{align} where $(W^{(1)},W^{(2)})$ is a two-dimensional standard Brownian motion (under the probability measure $P$) and $\mu, \rho, \kappa, \theta$ and $\sigma$ are constants. We assume that the Feller condition is satisfied, i.e. $$2 \kappa \theta > \sigma^2,$$ which ensures that $V_t >0.$
In Shreve's book, I read that the solution $(S_t,V_t)_{0 \leq t \leq T}$ to the two-dimensional SDE above is a Markov process but he doesn't prove it. I have already checked a couple of books and I only have found a sufficient condition, which requires the coefficients (drift and diffusion functions) to satisfy the Lipschitz and linear growth conditions. This is not the case for this SDE, so I don't know how to proceed. Any ideas?
Edit: I see in the comments asking for the definition of a Markov process. Any definition is fine as long as I can get a rigorous proof. For example:
The solution $(X_t,V_t)_{0 \leq t \leq T}$ of the above SDE is a Markov process if for any bounded Borel measurable function $f:\mathbb{R}^2 \rightarrow \mathbb{R}$ and for all $0 \leq s \leq t \leq \infty,$ we have $$E[f(X_t,V_t) | \mathscr{F}_s]=[E[f(X_t,V_t) |(X_s,V_s) ],$$ or we could also use the transition probability function of the Markov process.