<s>[INST] Determinant of big matrix Let:
$$M_n:=\left(\begin{array}{ccccc} 1 & 2 & 3 & \dots & n \\ 1^3 & 2^3 & 3^3 & \dots & n^3 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 1^{2n-1} & 2^{2n-1} & 3^{2n-1} & \dots & n^{2n-1} \end{array}\right),$$
for $n\in\mathbb N$. How could I compute $\det M_n$? I tried computing the first few and seeing if I saw an apparent pattern to try inducting, but the third one was already too much of a mess and no pattern was evident. This was an exercise in an exam, so there must be a smart way of doing it. Is there? Which is it?
 [/INST] $$M_n:=\left(\begin{array}{ccccc} 1 & 2 & 3 & \dots & n \\ 1^3 & 2^3 & 3^3 & \dots & n^3 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 1^{2n-1} & 2^{2n-1} & 3^{2n-1} & \dots & n^{2n-1} \end{array}\right) = \left(\begin{array}{ccccc} 1 \times 1^0 & 2 \times 2^0 & 3 \times 3^0 & \dots & n \times n^0 \\ 1 \times 1^2 & 2 \times 2^2 & 3 \times 3^2 & \dots & n \times n^2 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 1 \times 1^{2n-2} & 2 \times 2^{2n-2} & 3 \times 3^{2n-2} & \dots & n \times n^{2n-2} \end{array}\right)$$
Factor out $n!$ (in the form of $j$ from the $j$th column):
$$\text{det} = n! \left \vert \left(\begin{array}{ccccc} 1^0 & 2^0 & 3^0 & \dots & n^0 \\ 1^2 & 2^2 & 3^2 & \dots & n^2 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 1^{2n-2} & 2^{2n-2} & 3^{2n-2} & \dots & n^{2n-2} \end{array}\right) \right \vert$$
This is a Vandermonde matrix with $a_i = i^2$; the determinant of that matrix is $$\prod_{1 \leq i < j \leq n} (i^2 - j^2)$$
Therefore $$\text{det}(M_n) = n! \prod_{1 \leq i < j \leq n} (i+j)(i-j)$$
If you want to "simplify" it a bit more, you can do the following to get a feel for it, but it doesn't yield a nice closed form that I can see.
First, let $n$ be odd.
In the $i+j$ terms, $3$ and $4$ turn up once; $5$ and $6$ turn up twice; $7$ and $8$ three times; and so on up to $n$ and $n+1$, which occurs $\frac{n-1}{2}$ times. Then $n+2$ turns up $\frac{n-1}{2}$ times, $n+3$ and $n+4$ turn up $\frac{n-3}{2}$ times, and so on down to $2n-2$ and $2n-1$ which happen once. Therefore $$\prod_{1 \leq i < j \leq n} (i+j) = 3 \times 4 \times 5^2 \times 6^2 \times \dots \times n^{(n-1)/2} \times (n+1)^{(n-1)/2} \times (n+2)^{(n-1)/2} \times (n+3)^{(n-1)/2-1} \times \dots \times (2n-2) \times (2n-1)$$
You can do similar things if $n$ is even, and then analyse the $i-j$ terms in a similar way. (The $i-j$ terms lead to a product of factorials.)
</s>