write an array in latex

Question

I am unable to write this matrix in latex in a full width page. The matrix is given in the below link. Please help me to write this matrix in latex.

Edit: Thanks for your feedback. I did the following syntax, but unfortunately I am not able to get the output within a page. Help me.

\[\begin{array}{*{9}{c}} 
{1+\Delta t\theta } & {\rho _{,x} \Delta t} & {\rho_{,y} \Delta t} & {0} & {\rho \Delta t} & {0} & {0} & {0} & {0} \\ 
{\left(uu_{,x} +vu_{,y} \right)\Delta t} & {\rho +\rho u_{,x} \Delta t} & {\rho u_{,y} \Delta t} & {0} & {0} & {0} & {0} & {0} & {0} \\ 
{\left(uv_{,x} +vv_{,y} +\frac{1}{2\varepsilon {\kern 1pt} {\rm Fr}} \right)\Delta t} & {\rho v_{,x} \Delta t} & {\rho +\rho v_{,y} \Delta t} & {0} & {0} & {0} & {0} & {0} & {0} \\ 
{\left(uT_{,x} +vT_{,y} +\frac{\left(\gamma -1\right){\rm M}^{2} v}{2\varepsilon {\kern 1pt} {\rm Fr}} \right)\Delta t} & {\left(\rho T_{,x} -\left(\gamma -1\right){\rm M}^{2} p'_{,x} \right)\Delta t} & {\left(\rho T_{,y} -\left(\gamma -1\right){\rm M}^{2} p'_{,y} +\frac{\left(\gamma -1\right){\rm M}^{2} \rho }{2\varepsilon {\kern 1pt} {\rm Fr}} \right)\Delta t} & {\rho } & {-\frac{8}{3} \frac{\left(\gamma -1\right){\rm M}^{2} }{Re} \theta \Delta t} & {-2\frac{\left(\gamma -1\right){\rm M}^{2} }{Re} \omega \Delta t} & {0} & {0} & {-\left(\gamma -1\right){\rm M}^{2} } \\ 
{0} & {0} & {0} & {0} & {1} & {0} & {0} & {0} & {0} \\ 
{0} & {0} & {0} & {0} & {0} & {1} & {0} & {0} & {0} \\ 
{0} & {0} & {0} & {0} & {0} & {0} & {1} & {0} & {0} \\ 
{0} & {0} & {0} & {0} & {0} & {0} & {0} & {1} & {0} \\ 
{0} & {0} & {0} & {0} & {0} & {0} & {0} & {0} & {0} 
\end{array}\] 

Answer 1

\documentclass{amsart}
\begin{document}
\[
\begin{pmatrix}
{1+\Delta t\theta } & 
{\rho _{,x} \Delta t} & 
{\rho {,y} \Delta t} & 
{0} & 
{\rho \Delta t} & 
{0} & 
{0} & 
{0} & 
{0} \\ 
{\left(uu{,x} +vu_{,y} \right)\Delta t} & 
{\rho +\rho u_{,x} \Delta t} & 
{\rho u_{,y} \Delta t} & {0} & {0} & {0} & {0} & {0} & {0} \\ 
{\left(uv_{,x} +vv_{,y} +\frac{1}{2\varepsilon {\kern 1pt} {\rm Fr}} \right)\Delta t} & 
{\rho v_{,x} \Delta t} & 
{\rho +\rho v_{,y} \Delta t} & 
{0} & {0} & {0} & {0} & {0} & {0} \\ 
X &
Y & 
Z  & 
{\rho } & 
W & 
{-2\frac{\left(\gamma -1\right){\rm M}^{2} }{Re} \omega \Delta t} & 
{0} & {0} & 
{-\left(\gamma -1\right){\rm M}^{2} } \\ 
{0} & {0} & {0} & {0} & {1} & {0} & {0} & {0} & {0} \\ 
{0} & {0} & {0} & {0} & {0} & {1} & {0} & {0} & {0} \\ 
{0} & {0} & {0} & {0} & {0} & {0} & {1} & {0} & {0} \\ 
{0} & {0} & {0} & {0} & {0} & {0} & {0} & {1} & {0} \\ 
{0} & {0} & {0} & {0} & {0} & {0} & {0} & {0} & {0} 
\end{pmatrix}
\]
where
\begin{align*}
X
&=
{\left(uT_{,x} +vT_{,y} +\frac{\left(\gamma -1\right){\rm M}^{2} v}{2\varepsilon {\kern 1pt} {\rm Fr}} \right)\Delta t},
\\
Y
&=
{\left(\rho T_{,x} -\left(\gamma -1\right){\rm M}^{2} p'{,x} \right)\Delta t},
\\
Z
&=
{\left(\rho T{,y} -\left(\gamma -1\right){\rm M}^{2} p'_{,y} +\frac{\left(\gamma -1\right){\rm M}^{2} \rho }{2\varepsilon {\kern 1pt} {\rm Fr}} \right)\Delta t},
\\
W
&=
{-\frac{8}{3} \frac{\left(\gamma -1\right){\rm M}^{2} }{Re} \theta \Delta t},
\end{align*}
\end{document}