I would like a nod in the right direction with the following problem.
Let $0<t_1<t_2<\cdots<t_n<\infty$, and $\{a_i\}_1^n$ be real numbers, find a function $f(x_1,\ldots,x_n)$ so that
$$P\{B(t_1)\leq a_1, B(t_2)\leq a_2,\ldots, B(t_n)\leq a_n\} > =\int_{-\infty}^{a_n} \cdots \int_{-\infty}^{a_1} f(x_1,\ldots,x_n)\,dx_1\cdots dx_n$$
My assumption is that you try to re-write the expression so that you can use the fact that Brownian motion has independent increments, however I am getting nowhere.
$\newcommand{\var}{\operatorname{var}}\newcommand{\cov}{\operatorname{cov}}$ Notice that for $i\le j$ we have \begin{align} & \cov(B(t_i),B(t_j)) \\[10pt] = {} & \cov(B(t_i), (B(t_j)-B(t_i)) + B(t_i)) \\[10pt] = {} & \cov(B(t_i), B(t_j)- B(t_i)) + \cov(B(t_i),B(t_i)) \\[10pt] = {} & 0 + \var(B(t_i)) = t_i = t_{\min\{i,j\}} = \min\{t_i,t_j\}. \end{align}
Hence $$ \var\begin{bmatrix} B(t_1) \\ \vdots \\ B(t_n) \end{bmatrix} = \begin{bmatrix} t_1 & t_1 & t_1 & t_1 & \cdots & t_1 \\ t_1 & t_2 & t_2 & t_2 & \cdots & t_2 \\ t_1 & t_2 & t_3 & t_3 & \cdots & t_3 \\ t_1 & t_2 & t_3 & t_4 & \cdots & t_4 \\ \vdots & \vdots & \vdots & \vdots & & \vdots \\ t_1 & t_2 & t_3 & t_4 & \cdots & t_n \end{bmatrix}. $$ Call this matrix $V$. then the joint density of this random vector whose variance is $V$ is \begin{align} f(x_1,\ldots,x_n) = {} & \frac 1 {\sqrt{2\pi}^n}\cdot\frac 1 {(\det V)^{1/2}} \exp \left( \frac{-1} 2 x^\top V^{-1} x \right) \tag 1 \\[10pt] & \text{ where } x = \begin{bmatrix} x_1 \\ \vdots \\ x_n \end{bmatrix}. \end{align}
PS: If I'm not mistaken, then when $t_i=i$ for $i=1,2,3,\ldots$ then
PPS: Elementary row operations find the determinant: $t_1(t_2-t_1)(t_3-t_2) \cdots (t_n-t_{n-1})$. I think you can also give a geometric argument and a probabilistic argument for that result.