What is the sum of the coefficients of the equation?
There is a minimal volume four-dimensional ellipsoid with these points on its surface: {2,7,2,1}, {2,2,5,1}, {3,5,3,1}, {5,5,5,2} and {7,2,3,5}.
Find the equation of that ellipsoid with at least four digits to the right of the decimal point accuracy. The answer was developed with ten digits to the right of the decimal point.Then, sum the coefficients of the terms (include the constant term as a coefficient). The answer is entered with five digits to the right of the decimal point. Of course, the usual Brilliant allowance for error in a real number answer apples.
The methodology used to solve this problem is now documented in the Brilliant Wiki pages.
This section requires Javascript.
You are seeing this because something didn't load right. We suggest you, (a) try
refreshing the page, (b) enabling javascript if it is disabled on your browser and,
finally, (c)
loading the
non-javascript version of this page
. We're sorry about the hassle.
See Steiner ellipse, minimal area through three points and Affine transformations . What works in two dimensions works in one dimensional (but then the problem is trivial) and in higher dimensions than two-dimensional. The methodology remains the same.
The equation with rounded coefficients is 1 . x 1 2 − 0 . 7 7 7 2 5 4 x 2 x 1 − 1 . 1 8 0 9 1 x 3 x 1 − 2 . 9 4 9 9 1 x 4 x 1 + 5 . 8 1 5 5 6 x 1 + 0 . 2 5 5 1 5 6 x 2 2 + 0 . 5 4 5 6 6 9 x 3 2 + 2 . 3 0 2 3 x 4 2 − 4 . 1 0 3 1 2 x 2 + 0 . 6 4 2 3 1 x 2 x 3 − 5 . 9 3 2 8 2 x 3 + 1 . 3 0 0 5 3 x 2 x 4 + 1 . 8 9 6 8 8 x 3 x 4 − 1 0 . 2 9 0 5 x 4 + 1 7 . 8 7 4 5 = 0
The ends of the ellipsoid's axes with rounded coordinates: { 2 . 4 2 0 4 4 , 3 . 0 3 4 2 9 , 5 . 5 4 3 6 5 , 0 . 4 1 8 1 1 7 } { 5 . 1 3 5 4 7 , 5 . 8 9 9 3 6 , 3 . 3 7 6 4 7 , 2 . 1 2 3 0 3 } { 5 . 2 7 8 4 3 , 1 . 6 8 4 6 5 , 3 . 2 0 8 6 2 , 3 . 5 3 2 6 2 } { 0 . 7 7 2 1 9 9 , 6 . 3 1 3 3 7 , 2 . 2 6 6 1 8 , − 0 . 1 7 5 6 8 4 } { 5 . 1 7 9 5 6 , 5 . 3 6 5 7 1 , 1 . 6 5 6 3 5 , 3 . 5 8 1 8 8 } { 2 . 4 6 4 5 3 , 2 . 5 0 0 6 4 , 3 . 8 2 3 5 3 , 1 . 8 7 6 9 7 } { 2 . 3 2 1 5 7 , 6 . 7 1 5 3 5 , 3 . 9 9 1 3 8 , 0 . 4 6 7 3 8 4 } { 6 . 8 2 7 8 , 2 . 0 8 6 6 3 , 4 . 9 3 3 8 2 , 4 . 1 7 5 6 8 }
OK, OK OK, how did I get there:
The n array, the regular 5-cell vertex homogeneous coordinates in radial form and which locations are on a unit four-dimension sphere it is easy to prove that the unit regular 5-cell has the maximum volume within that sphere: ⎝ ⎜ ⎜ ⎜ ⎜ ⎜ ⎛ 1 0 0 0 − 4 1 4 1 5 0 0 − 4 1 − 4 3 5 6 5 0 − 4 1 − 4 3 5 − 2 6 5 − 2 2 5 − 4 1 − 4 3 5 − 2 6 5 2 2 5 ⎠ ⎟ ⎟ ⎟ ⎟ ⎟ ⎞
From here on, I am providing only the numeric forms of the matrices as the radical forms take minutes to compute and take up pages to present.
The same matrix in numeric form: \left( \begin{array}{|l|l|l|l|l|} \hline 1. & -0.25 & -0.25 & -0.25 & -0.25 \ \hline 0. & 0.968246 & -0.322749 & -0.322749 & -0.322749 \ \hline 0. & 0. & 0.912871 & -0.456435 & -0.456435 \ \hline 0. & 0. & 0. & -0.790569 & 0.790569 \ 1. & 1. & 1. & 1. & 1. \ \hline \end{array} \right)]
The inverse of n in numeric form: ⎝ ⎜ ⎜ ⎜ ⎜ ⎛ 0 . 8 − 0 . 2 − 0 . 2 − 0 . 2 − 0 . 2 0 0 . 7 7 4 5 9 7 − 0 . 2 5 8 1 9 9 − 0 . 2 5 8 1 9 9 − 0 . 2 5 8 1 9 9 0 0 0 . 7 3 0 2 9 7 − 0 . 3 6 5 1 4 8 − 0 . 3 6 5 1 4 8 0 0 0 − 0 . 6 3 2 4 5 6 0 . 6 3 2 4 5 6 0 . 2 0 . 2 0 . 2 0 . 2 0 . 2 ⎠ ⎟ ⎟ ⎟ ⎟ ⎞
The equation for a four-dimensional sphere: x 1 2 + x 2 2 + x 3 2 + x 4 2 = 1 .
The o , which contains the five four-dimensional points in homogeneous coordinate form: ⎝ ⎜ ⎜ ⎜ ⎜ ⎛ 2 . 7 . 2 . 1 . 1 . 2 . 3 . 3 . 1 . 1 . 2 . 5 . 2 . 1 . 1 . 5 . 3 . 3 . 1 . 1 . 7 . 2 . 3 . 5 . 1 . ⎠ ⎟ ⎟ ⎟ ⎟ ⎞
At this point, the determinant of o is checked to be non-zero so that it is invertible, i. e., the points are not contained in a subspace.
At this point both o and n are known and therefore a can be computed from the a . o = n by inverting o and right multiplying n giving a : \[\left( \begin{array}{|l|l|l|l|l|} \hline -0.707547 & 0.377358 & 0.212264 & 0.990566 & -1.64151 \\ \hline -0.645497 & 0. & 0.322749 & 0.645497 & 0. \\ \hline 0.766467 & -0.439211 & -0.732019 & -1.4382 & 4.44379 \\ \hline -0.074582 & -0.223746 & -0.37291 & 0.104415 & 2.35679 \\ \hline 0. & 0. & 0. & 0. & 1. \\ \hline
\end{array} \right)\]
The inverse of a can be computed by inverting a : ⎝ ⎜ ⎜ ⎜ ⎜ ⎛ − 1 . 8 2 . 8 − 1 . 6 − 1 . 0 . − 2 . 3 2 3 7 9 − 1 . 5 4 9 1 9 1 . 0 3 2 8 − 1 . 2 9 0 9 9 0 . − 2 . 1 9 0 8 9 1 . 0 9 5 4 5 − 0 . 7 3 0 2 9 7 − 1 . 8 2 5 7 4 0 . 1 . 2 6 4 9 1 − 1 . 8 9 7 3 7 − 1 . 2 6 4 9 1 1 . 8 9 7 3 7 0 . 3 . 8 4 . 2 3 . 6 2 . 1 . ⎠ ⎟ ⎟ ⎟ ⎟ ⎞
At this point, the Singular Value Decomposition of the affine matrix portion of the inverse of a is computed to get the scale factors of the ellipsoid's axes (the diagonal of the σ and the rotation to the directions of those axes (\bm{v}). By the way, this is where the eigenvalues and eigenvectors mentioned by Mark Henning come in to the solution: u = ⎝ ⎜ ⎜ ⎛ 0 . 6 4 4 9 3 9 − 0 . 5 4 0 8 4 6 0 . 0 9 6 9 2 0 8 0 . 5 3 1 1 7 4 0 . 4 1 8 5 9 7 0 . 6 9 2 4 9 6 − 0 . 5 1 0 9 5 6 0 . 2 9 0 0 8 7 − 0 . 4 0 1 7 5 5 − 0 . 4 2 2 3 4 − 0 . 7 8 7 1 8 4 0 . 2 0 1 4 0 5 0 . 4 9 7 4 1 7 − 0 . 2 2 2 6 3 1 − 0 . 3 3 1 4 6 9 − 0 . 7 7 0 1 5 6 ⎠ ⎟ ⎟ ⎞ diagonal ( σ ) = { 5 . 4 3 7 7 7 , 3 . 4 4 8 5 7 , 1 . 9 5 0 9 1 , 0 . 4 1 4 6 4 2 } v = ⎝ ⎜ ⎜ ⎛ − 0 . 6 1 8 1 7 7 − 0 . 2 2 9 2 2 5 − 0 . 5 6 0 1 6 1 0 . 5 0 1 5 3 1 0 . 4 9 6 7 1 5 − 0 . 8 5 4 7 7 5 − 0 . 0 9 1 3 3 6 8 0 . 1 1 9 5 5 2 0 . 3 0 6 8 8 1 0 . 2 6 3 9 1 0 . 3 2 0 2 1 7 0 . 8 5 6 5 2 6 − 0 . 5 2 6 2 6 6 − 0 . 3 8 3 6 2 5 0 . 7 5 8 5 1 1 0 . 0 2 3 1 8 1 3 ⎠ ⎟ ⎟ ⎞
Yes, I checked the a and a − 1 were, in fact, inverses of each other.
Using a , { x 1 , x 2 , x 3 , x 4 } were projected into the unit sphere environment and substituted into the Unit sphere equation ( x 1 2 + x 2 2 + x 3 2 + x 4 2 = 1 ) and then the resultant equation was simplified. The resultant equation was given above and not repeated here.
Forward mapping the positive and negative unit vectors using the v matrix and using the inverse of the a back the original environment of the specified points gives the ends of the axes of the ellipsoid. These also were provided above and not repeated here