Hacker News Reader

Wouldn’t it be even faster to solve that problem conditionally? P(x >= a)P(y >= a | x >= a), since the conditional distribution is known and normal?

You still need to integrate over x. y|x=c is Gaussian for each c, but there are many such c's in [a,inf].

If you're curious:

Ah, duh. Not enough coffee.