A famous conjecture (due I think to Hardy and Littlewood) states that if $P(x)$ denotes the number of primes of the form $n^2+1$ less than or equal to $x$, then $$P(x)\sim \frac{C\sqrt x}{\log x}$$ for some constant $C$ (their Conjecture F gives an expression for $C$, which I think is predicted to be somewhere in the area of $2$, although I'm not confident of that.
However, having calculated the first $1.35$ million such primes, it is not at all obvious that this is accurate. My question is the following: is there substantial numerical evidence which supports the conjecture? If so, is it simply not apparent until much higher values, or are my calculations wrong?
Up to $10^{24}$ there are $25814570672$ primes of the form $n^2+1,$ and the formula with $C=1.3728134628182\ldots$ predicts $24841887983$ (3.9% error). You can get better by replacing $C\sqrt x/\log x$ with $C\operatorname{li}(\sqrt{x})/2$ which gives $25814350227$ (0.00085% error). Note that $C\sqrt x/\log x \sim C\operatorname{li}(\sqrt x)/2.$
You mentioned the first 1.35 million primes of the form $n^2+1.$ The 1.35 millionth prime of that form is $1016179851635777$, and the two functions predict 1.266450 million and 1.350213 million, respectively -- not bad at all.
So yes, there is substantial numerical evidence in favor of this conjecture.