You might consider the "sum-of-squares" approach. The idea is to find a set of polynomials so that your quadratic formexpression is the sum of the squaresquares of the elements in the region of interest. For your case, you could replace each $x_i$ with a new variable $z_i^2$, and ask; you are now asking if the quadratic form is a SOS as an (unconstrained)corresponding 4-th degree unconstrained polynomial over $z$is non-negative.
This restatement may not sound like an improvement, but it turns out that SOS problems can be attacked using semidefinite programming techniques (e.g., see this page). You can use a freely available SDP solvers.
This is a sufficient approach, i.e., it may prove that your original quadratic form is positive, but it can't disprove it. Since you're trying to solve a specific problem, though, it may be worth the gamble.