I came across this question while preparing a talk on [this paper](http://imrn.oxfordjournals.org/content/early/2015/09/09/imrn.rnv251.full.pdf). There we gave the explicit formula for $p_n$ (in the GOE) that the other answers anticipated could be found by random matrix methods. The formula is given as follows: First, let $n\geq 1$ be any integer, and define $n':=2 \lceil n/2 \rceil$. When real symmetric matrices are chosen according to the $n$-dimensional Gaussian Orthogonal Ensemble, the probability of positive definiteness is given by: \begin{equation}p_n=\frac{{\mathrm{ Pf}}(A)}{2^{n(n+3)/4}\prod_{m=1}^n\Gamma(\frac m 2)},\label{theorem3} \end{equation} where $A$ is the $n'\times n'$ skew-symmetric matrix whose $(i,j)$-entry $a_{ij}$ is given for $i<j$ by: \begin{align}\label{pfaffianformula} a_{ij} = \begin{cases} \,2^{i+j-2}\Gamma(\frac{i+j}{2})\left(\beta_{\frac{1}{2}}(\frac{i}{2},\frac{j}{2})- \beta_{\frac{1}{2}}(\frac{j}{2},\frac{i}{2})\right) & \mbox{if }~ i<j\leq n, \\[.075in] \,2^{i-1}\Gamma(\frac{i}{2}) & \mbox{if }~ i<j=n+1.\\[.075in] \end{cases} \end{align} $($Note that the second case arises only when $n$ is odd.$)$ Here $\mathrm{Pf}$ is the Pfaffian, $\Gamma$ denotes the usual gamma function $\Gamma(s)=\int_{0}^\infty x^{s-1}e^{-x} \mathrm{d} x$, and $\beta_{t}$ is the usual incomplete beta function $\beta_{t}(i,j) = \int_0^{t} x^{i-1}(1-x)^{j-1} \mathrm{d} x.$ This formula follows from de Bruijn's identity for evaluating determinantal integrals. As for the asymptotics, the leading term of $\log p_n$ is indeed $-\frac{\log(3)}{4} n^2$ as noted in Mikael's answer. If you are interested in the leading order asymptotic for $p_n$ it can be found from the results of [Borot, Eynard, Majumdar and Nadal](https://arxiv.org/pdf/1009.1945v4.pdf). Specialising to the GOE they give: \begin{align} p_n \sim \exp\left(-\frac{\log3}{4} n^2 - \frac{1}{2}\log\left(1+\frac{2}{\sqrt{3}}\right)n - \frac{1}{24}\log n + \alpha_1 \right)(1+o(1)).\label{asympexp} \end{align} The constant $\alpha_1$ is given by: \begin{align} \alpha_1= \frac{\zeta'(-1)}{2} - \frac{\log(2)}{12} -\frac{\log(3)}{16}+ \frac{1}{4}\log\left(1+\frac{2}{\sqrt{3}}\right)\simeq-0.0172, \end{align} where $\zeta(z)$ is the Riemann zeta function.