我正在尝试使用有限差分方法和迭代求解器（Gauss-Seidel）求解以下方程组。

连续性方程：

$\frac{\partial V_x}{\partial x} + \frac{\partial V_y}{\partial y} = 0$

简化的斯托克斯方程：

$\eta (\frac{\partial^2 V_x}{\partial x^2} + \frac{\partial^2 V_x}{\partial y^2}) - \frac{\partial P}{\partial x} = 0$

$\eta (\frac{\partial^2 V_y}{\partial x^2} + \frac{\partial^2 V_y}{\partial y^2}) - \frac{\partial P}{\partial y}+ \rho g = 0$

我是否应该使用两个斯托克斯方程来求解速度（）。但是，我该如何解决压力（）？还是我首先使用连续性方程求解速度，最后使用斯托克斯方程求解压力？$v_x,v_y$$P$