Posts Tagged ‘patch’

On Calculating Entropy

June 10th, 2014 3 comments

When talking about Pasquali patches that converge and stationary surfaces that they converge to, we discussed that such would stabilize because it wouldn't loop in any shape or form nor would it diverge.  Such Claim we called Entropy.  We didn't provide a mechanism to calculate such. We do that now.

The way we propose we calculate the entropy is by measuring the deviation from the stationary surface.  Thus for the case where we have a Pasquali patch p_t(x,y), t \in \mathbb{Z}^+ we could measure entropy by:

 S(t) = \int_{[0,1]^2} \vert p_t(x,y) - p_\infty(x) \vert \, dA

that is, by taking the absolute deviation from the stationary surface (notice we use the area element dA).  Since in this case the time element is configured at integer times, rather than on continuous time, letting  t \in \mathbb{R} implies using the Pasqualian.  Thus we would have 

 S(t) = \int_{[0,1]^2} \vert p(x,y,t) - p_\infty(x) \vert \, dA

It seems fairly clear that 0 < S(t) < \infty , but here we must agree that 0 means that the system is stable and has the most entropy.

We may however, rather than the absolute difference from the stationary surface, choose to emphasize larger changes and (heavily) discount smaller ones.  Thus we could suggest

 S^{\square}(t) = \int_{[0,1]^2} \left( p(x,y,t) - p_\infty(x) \right)^2 \, dA

for t \in \mathbb{R}.


On Infinite Term Functions with Converging Integrals, Part II

May 5th, 2014 No comments

This post is a continuation of a previous one. I have developed several preparatory (in preparation of other) claims or theorems:  The first and second claim show that a particular collection of finite polynomial functions have area of 1 in the interval [0,1], and are hence also Pasquali patches.  The third and fourth shows that the finite sum of any such functions actually have converging integrals in the same interval.  The corollaries show that this is not the case if the sum is infinite.  A fun way to summarize this information is soon after developed, and these observations, though simple, lead us to classify all Pasquali patches which are functions of x alone, and therefore all stationary/limiting/stable surfaces, eigenfunctions or wavevectors (from the quantum mechanics point-of-view).

Claim 1. Take f_i(x) = (i+1) x^i with i = 0 \ldots n.  Then \int_0^1 f_i(x) \, dx = 1, \forall i.

Proof by Definition of Integration (Inducing).  We show that \int_0^1 f_0(x) \, dx = 1. The expression equals 

\int_0^1 x^0 \, dx = \int_0^1 1 \, dx = x \left. \right\vert_0^1 = 1

 We assume that the kth element \int_0^1 f_k(x) \, dx = 1 although we readily know by the definition of integration that such is true, since 

\int_0^1 (k+1) x^k \, dx = x^{k+1} \left. \right\vert_0^1 = 1^{k+1} = 1

The exact same definition argument applies to the k+1th element and 

\int_0^1 (k+2) x^{k+1} \, dx = x^{k+2} \left. \right\vert_0^1 = 1^{k+2} = 1

Claim 2. The functions f_i(x) = (i+1) x^i with i = 0 \ldots n are Pasquali patches.

Proof.  A Pasquali patch is a function p(x,y) so that \int_0^1 p(x,y) dx = 1.  Let p(x,y) = f_i(x).  Since by Claim 1  \int_0^1 f_i(x) \, dx = 1, \forall i = 1 \ldots n, then applying the definition means  f_i(x) = (i+1) x^i are Pasquali patches \forall i= 0 \ldots n .

Claim 3.  The finite polynomial g(x) = \sum_{i=0}^n (i+1) x^{i} converges in area from [0,1] to n+1.

Proof.  We are looking for

\int_0^1 \sum_{i=0}^n (i+1) x^i \, dx

.  The sum is finite so it converges, and there is no issue exchanging the order of the sum and integral. Thus:

\sum_{i=0}^n \int_0^1 (1+i) x^i \, dx =\sum_{i=0}^n \left( x^{i+1} \left. \right\vert_0^1 \right) = \sum_{i=0}^n 1^{i+1} =\sum_{i=0}^n 1 = n+1

Claim 4. Pick n functions from the pool of f_i(x) = (i+1) x^i.  For example, pick f_3(x), f_5(x), and f_7(x).  Create the function h(x) = \sum_i f_i(x).  Then \int_0^1 h(x) \, dx = n.

Proof by induction.  Since by Claim 2 all f_i(x) are Pasquali patches, it follows their integral is 1 in the interval (Claim 1).  Picking 1 function from the pool thus gives an integral of 1 in the interval.  Suppose that picking k functions gives k units at the integral in the interval. Now pick k+1 functions.  The first k functions give k units at the integral in the interval, and the 1 additional function contributes 1 unit at the integral in the interval.  Thus k+1 functions contribute k+1 units at the integral in the interval.

Corollary 1. The infinite polynomial a(x) = \sum_{i=0}^\infty (i+1) x^i diverges in area in the interval from [0,1].

Proof.  Take

\int_0^1 \left( \lim_{n \to \infty} \sum_{i=0}^n (1+i) x^i \right) \, dx =\lim_{n \to \infty} \int_0^1\sum_{i=0}^n (1+i) x^i \, dx

Here exchanging the order of limit and integral is justified by the fact that, term-wise, the integral converges. Next  

\lim_{n \to \infty} n+1 = \infty

Here the second to last step is justified by Claim 3.

Corollary 2.  The infinite polynomial a(x) - h(x) diverges in area in the interval from [0,1].

Proof.  Take the limit

 \lim_{n \to \infty} \left[ a(x) - h(x) \right]

Taking n to infinity applies to a(x) only which we know diverges by Corollary 1.  The same limit  has no effect on h(x) as the sum it is composed of is finite and adds up to an integer constant, say m.  We conclude that any infinite collection of terms of f_i(x) diverges, even when a finite number of them may be absent from the sum.

And now sushi.

Corollary 3.  The infinite polynomial  a(x) - b(x) diverges in area in the interval from [0,1] with a(x), b(x) are infinite polynomials constructed by sums of functions picked from the pool f_i(x) = (i+1) x^i and with no repetitions. (Note that the difference of these two infinite polynomials must also be infinite).

Proof. Since the a(x) - b(x) is an infinite polynomial, the integral of such will be an infinite string of ones since the functions it contains are f_i(x) and these are Pasquali patches (Claim 2) and there are no repetitions.  Such infinite sum of ones clearly diverges.

Remark 1.  We can view what we have learned in the claims from a slightly different vantage point.  Create the infinite identity matrix

 I = \left[ \begin{array}{cccc} 1 & 0 & 0 & \ldots \\ 0 & 1 & 0 & \ldots \\ \vdots & \vdots & \vdots & \ddots \end{array} \right]

Next create the following polynomial differential vector

 D =\left[ \begin{array}{c} 1 \\ 2x \\ 3x^2 \\ \vdots \end{array} \right]

It is clear that

 \int_0^1 I_i \cdot D \, dx =1

for all rows  i of  I .  We can omit the little  i because this definition applies to all rows and:

 \int_0^1 I \cdot D \, dx = \int_0^1 D \, dx= \left[ \begin{array}{c} 1 \\ 1 \\ \vdots \end{array} \right] = \bf{1}

This of course summarizes Claims 1 and 2.  Next, define the matrix J consisting of rows which are finite sums of rows of I (so that each row of J consists of a finite number of ones at any position, namely n such coming from n picked rows of I).  Claims 3 and 4 are summarized in the statement  

 \int_0^1 J\cdot D \, dx = S

where S is the vector consisting of the sum of the rows of J, which, since it is made up of a finite number of ones at each row, adds up to a constant integer at each row:

S = \left[ \begin{array}{c} n_1 \\ n_2 \\ \vdots \end{array} \right]

 Finally, the corollaries can be summarized in the statement in which we create a matrix  K consisting of rows with a finite number of zeroes (and an infinite number of ones) or an infinite number of zeroes but an infinite number of ones as well.  It is clear then that

 \int_0^1 K\cdot D \, dx = \infty

Remark 2. The cool thing about this notation is that it gives us power to conclude several interesting things.  For example, scaling of matrices  I and  J as by a constant  t shows convergence at the integral in the interval  \left[ 0,1 \right] of every one of the scaled sums  represented by the rows of such matrices.  Thus:

Corollary 4. Let  I^* = t \cdot I and J^* = t \cdot J with  t is a scaling factor.  Then the area of each of the infinitely many polynomials represented by the matrices I^*, J^* dot D in the interval from 0 to 1 converge.

Proof.  On the one hand, we have 

 \int_0^1 I^* \cdot D \, dx = \int_0^1 t \cdot I \cdot D \, dx =t \left( \int_0^1 I \cdot D \, dx \right) = \bf{t}

 On the other hand,

 \int_0^1 J^* \cdot D \, dx =\int_0^1 t \cdot J\cdot D \, dx = t \left( \int_0^1 J\cdot D \, dx \right) = t \cdot S =\left[ \begin{array}{c} t \cdot n_1 \\ t \cdot n_2 \\ \vdots \end{array} \right]

Remark 3. Next consider the infinite-matrix formed by convergent sequences (at the sum) at each row,

A = \left[ \begin{array}{cccc} \vdots & \vdots & \vdots & \vdots \\ 1 & \frac{1}{2^2} & \frac{1}{3^2} & \ldots \\ \vdots & \vdots & \vdots & \vdots \end{array} \right]

Depicted is the reciprocals of squares which we know converges at the sum (Basel problem), simply for illustration, but all convergent sequences would be in the ith row of A.  We have 

 \int_0^1 A_i\cdot D \, dx = \sum_j a_{i,j}

is convergent by definition.  The cool thing is we can easily prove in one swoop that all sequences that are scaled will also converge at the sum (and the infinite polynomials with coefficients A \cdot D have converging area in the interval from 0 to 1).

Corollary 5. Let  A^* = t \cdot A with  t is a scaling factor.  Then the area of each of the infinitely many polynomials represented by the matrix entries of A^* \cdot D in the interval from 0 to 1 converge.

Proof.  We have 

 \int_0^1 A^*_i\cdot D \, dx = \sum_j a_{i,j}

for all i, so this equals 

 \int_0^1 t \cdot A_i\cdot D \, dx = t \left(\int_0^1 A_i\cdot D \, dx \right) = t \cdot\sum_j a_{i,j}

for all i.

All of these small and obvious observations lead to this:

Claim 5. The Grand Classification Theorem of Limiting Surfaces (A General and Absolutely Complete Classification of Pasquali patches which are functions of x alone).  All Pasquali patches which are functions of x alone (and therefore possible limiting surfaces) take the form

 p(x) = \frac{A_i \cdot D}{\sum_j a_{i,j}}

Proof. We have that, since such  p(x) is a Pasquali patch, it must conform to the definition.  Thus 

 \int_0^1 p(x) \, dx = \int_0^1 \frac{A_i \cdot D}{\sum_j a_{i,j}} \, dx = \frac{\int_0^1 A_i \cdot D \, dx}{\sum_j a_{i,j}} = \frac{\sum_j a_{i,j}}{\sum_j a_{i,j}} = 1

shows this is indeed the case.  To show that "all" Pasquali patches that are functions of x alone are of the form of p(x), we argue by contradiction.  Suppose that there is a Pasquali patch that is a function of x alone which does not take the form of p(x).  It couldn't possibly be one such that is a finite polynomial, since  A_i was defined to be that matrix formed by all convergent sequences at the sum at each row and it can be scaled any which way we like, and this includes sequences with a finite number of nonzero coefficients.  But now it couldn't be any infinite polynomial either, by the same definition of  A_i which includes infinite sequences so that \sum_j a_{i,j} is convergent.  Thus it must be a polynomial formed by dotting divergent sequences (at the sum), but all such have been happily excluded from the definition of A.

Remark 4.  Thus, EVERY convergent series has an associated Pasquali patch (which is solely a function of  x), and vice versa, covering the totality of the Pasquali patch functions of x universe and the convergent series universe bijectively.

Remark 5.  Notice how the definition takes into account Taylor polynomial coefficients (thus all analytic functions are included) and those that are not (even those that are as yet unclassified), and all sequences which may be scaled by a factor as well.

Claim 6. Let f(x) is Maclaurin-expandable so that

 f(x) = \sum_{n=0}^\infty \frac{f^n(0) x^n}{n!}


 \sum_{n=0}^\infty\frac{f^n(0)}{(n+1)!} = \int_0^1 f(x) \, dx


\int_0^1 f(x) \, dx = \int_0^1 A_i \cdot D \, dx

for some i row of A.  Such a row would have to be of form

 A_i = \left[ \begin{array}{cccc} f(0) & \ldots & \frac{f^n(0)}{n! (n+1)} & \ldots \end{array} \right]

 Then the integral

\int_0^1 A_i \cdot D \, dx = \sum_j a_{i,j} =\sum_{n=0}^\infty \frac{f^n(0)}{n! (n+1)} = \sum_{n=0}^\infty \frac{f^n(0)}{(n+1)!}

Remark 6. Notice that all Maclaurin-expandable functions converge in area (have stable area) in the interval from 0 to 1, a remarkable fact.

Example 1.  Take

f(x) = e^x = \sum_{n=0}^\infty \frac{x^n}{n!}

 By applying Claim 6, it follows that

 \sum_{n=0}^\infty \frac{1}{(n+1)!} = \int_0^1 e^x \, dx = e - 1

Remark 7. Now we have a happy way to construct (any and all) Pasquali patches which are functions of x alone, merely by taking a sequence which is convergent at the sum.

Remark 8. Quantum mechanically, we now know all possible shapes that a stationary (limiting) eigen wavevector can take.

Remark 9. This gives us extraordinary power to calculate convergent sums via integration, as the next examples show.  It also gives us extraordinary power to express any number as an infinite sum, for example.


Compendium of Claims and Proofs, Including New Ones, Part I

December 3rd, 2012 No comments

I've condensed this exceptional mathematical wisdom here, which is still transforming as I organize and jot down ideas.

Part I v16 (latest, but very unorganized after @Dynamics)

Part I v15

Part I v14

Part I v13

Part I v12

Part I v11

Part I v10

Part I v9

Part I v8

Part I v7

Part I v6

Part I v5

Part I v4

Part I v3

On Patch Stationariness II (RWLA,MCT,GT,AM Part X)

December 2nd, 2012 No comments

Claim:  Suppose p(x,y) = f_1(x) \cdot g_1(y) + f_2(x) \cdot g_2(y) and \int p(x,y) dx = 1 .  This last restriction causes freedom of choice of three functions, say the first three, and then g_2(y) = \frac{1 - g_1 (y) \int_0^1 f_1(x) dx}{\int_0^1 f_2(x) dx} = \frac{1 - g_1(y) F_1}{F_2} .  This basic construction of p(x,y) is the premise of all claims regarding patches I've proven here on my blog.  Take an a(x) with \int_0^1 a(x) dx = 1, with form that will be defined in the course of the proof.   Then B = \int a(1-y) g_1(y) dy = \sum_{i=0}^{\infty} a^i(1-y) G_1^{i+1} \vert_0^1 provided the sum converges.  If it does, this gives rise to the "stationary" p_\infty (x,y) = \frac{f_2(x)}{F_2} - \left( \frac{f_2(x) F_1}{F_2} - f_1(x) \right) B (and it is also a patch).

Proof: Since we are looking for stationary p_\infty (x,y), we are looking for  a(x) \star p(x,y) = \lambda a(x) with  \lambda = 1 in direct analogy with discrete Markov chains.  The star operator is explicitly as I've defined it previously implies \int_0^1 a(1 - y) p(x,y) dy = a(x) .  By the definition of  p(x,y) , we have that \int_0^1 a(1 - y) \left( f_1(x) g_1(y) + f_2(x) \left( \frac{1 - g_1(y) F_1}{F_2} \right) \right) dy = a(x) .  Expansion results in a(x) = f_1(x) \int_0^1 a(1-y) g_1(y) dy + \frac{f_2(x)}{F_2} \cdot 1- \frac{F_1}{F_2} f_2(x) \int_0^1 a(1-y) g_1(y) dy , where we have simplified  \int_0^1 a(1-y) dy to 1 because the transformation to the y-axis does not change the integral result.

Next we have the expression

a(x) = \frac{f_2(x)}{F_2} - \left( \frac{F_1 f_2(x)}{F_2} - f_1(x) \right) \int_0^1 a(1-y) g_1(y) dy


a(x) = \frac{f_2(x)}{F_2} - \left( \frac{F_1 f_2(x)}{F_2} - f_1(x) \right) B and derivatives

a^i(x) = \frac{f_2^i(x)}{F_2} - \left( \frac{F_1 f_2^i(x)}{F_2} - f_1^i(x) \right) B

We want to obtain B, so that the expression  \int_0^1 a(1-y) g_1(y) dy can be clearly defined.  We use the tabular method to simplify the integration by parts.

 \begin{array}{ccccc} \vert & Derivatives & \vert & Integrals & \vert \\ \vert & a(1-y) & \vert & g_1(y) & \vert \\ \vert & -a'(1-y) & \vert & G_1^1(y) & \vert \\ \vert & a''(1-y) & \vert & G_1^2(y) & \vert \\ \vert & \vdots & \vert & \vdots & \vert \end{array}

and obtain

 B = a(1-y) G_1^1(y) + a^{'}(1-y) G_1^2(y) + \ldots \vert_0^1 = \sum_{i=0}^\infty a^i (1-y) G_1^{i+1} \vert_0^1

Lastly the steady-state "stationary" patch is exactly this a(x), so that p_\infty(x,y) = a(x) as we desired to show.

Notice that \int_0^1 p_\infty(x,y) dx = 1 by hypothesis, since it is equal to a(x).

Existence.  Lastly, notice that p_\infty (x,y) exists provided the sum converges.  In previous exercises, I have shown that such converges in two specific examples: when derivatives of a(x) vanish (when a(x) is a finite polynomial) and when derivatives of a(x) are periodic but vanish for the periodic terms.

On Eigen(patch(ix))values, II - (RWLA,MCT,GT,AM Part IX)

March 22nd, 2011 No comments

So remember my little conjecture from last time, that the number of patch(ix) (kernel) eigenvalues would depend on the number of x terms that composed it?  I started working it out by writing all expressions and trying to substitute them and I got sums of sums of sums and it became nightmarish, and since math is supposed to be elegant, I opted for a different track.  A little proposition did result, but I'm not sure yet if it means what I want it to mean. Haha.

If you recall, last time we figured that

 B_1 = \frac{B_2 \sum_{i=0}^\infty f_2^i(1-y)G_1^{i+1}(y)\vert_0^1}{\lambda - \sum_{i=0}^\infty f_1^i(1-y)G_1^{i+1}(y)\vert_0^1}


 B_2 = \frac{B_1 \sum_{i=0}^\infty f_1^i(1-y) G_2^{i+1}(y)\vert_0^1}{\lambda - \sum_{i=0}^\infty f_2^i(1-y) G_2^{i+1}(y)\vert_0^1}

Let's rename the sums by indexing over the subscripts, so that

 \begin{array}{ccc} C_{1,1} & = &\sum_{i=0}^\infty f_1^i(1-y)G_1^{i+1}(y)\vert_0^1 \\ C_{1,2} & = &\sum_{i=0}^\infty f_1^i(1-y) G_2^{i+1}(y)\vert_0^1 \\ C_{2,1} & = &\sum_{i=0}^\infty f_2^i(1-y)G_1^{i+1}(y)\vert_0^1 \\ C_{2,2} & = &\sum_{i=0}^\infty f_2^i(1-y) G_2^{i+1}(y)\vert_0^1 \end{array}

Renaming therefore the constants we get:

 B_1 = \frac{B_2 C_{2,1}}{\lambda - C_{1,1}}


 B_2 = \frac{B_1 C_{1,2}}{\lambda - C_{2,2}}

Last time we substituted one equation into the other to figure additional restrictions on  \lambda .  A faster way to do this is to  notice:

 \left( \lambda - C_{1,1} \right)B_1 = B_2 C_{2,1}


 \left( \lambda - C_{2,2} \right) B_2 = B_1 C_{1,2}

If we multiply these two expressions we get

 \left( \lambda - C_{1,1} \right)\left( \lambda - C_{2,2} \right) B_1 B_2 = B_1 B_2 C_{1,2} C_{2,1}

Finally, dividing out both  B_1, B_2 we arrive at the quadratic expression on  \lambda of before:

 \left( \lambda - C_{1,1} \right)\left( \lambda - C_{2,2} \right) = C_{1,2} C_{2,1}

Now.  Let's posit that, instead of  a(x) = B_1 f_1(x) + B_2 f_2(x) we have  a^*(x) = B_1 f_1(x) + B_3 f_3(x) .  Then by all the same arguments we should have an expression of  B_1 that is the same, and an expression of  B_3 that is:

 \left( \lambda - C_{3,3} \right) B_3 = B_1 C_{1,3}

with the similar implication that

 \left( \lambda - C_{1,1} \right)\left( \lambda - C_{3,3} \right) = C_{1,3} C_{3,1}

An  a^{**}(x) = B_2 f_2(x) + B_3 f_3(x) would give the implication

 \left( \lambda - C_{2,2} \right)\left( \lambda - C_{3,3} \right) = C_{2,3} C_{3,2}

If we are to multiply all similar expressions, we get

 \left( \lambda - C_{1,1} \right)^2\left( \lambda - C_{2,2} \right)^2 \left( \lambda - C_{3,3} \right)^2 = C_{1,2} C_{2,1}C_{1,3} C_{3,1}C_{2,3} C_{3,2}


 \left( \lambda - C_{1,1} \right) \left( \lambda - C_{2,2} \right) \left( \lambda - C_{3,3} \right) = \sqrt{C_{1,2} C_{2,1}C_{1,3} C_{3,1}C_{2,3} C_{3,2}}

In other words, we want to make a pairwise argument to obtain the product of the  \lambda -expressions and a polynomial in  \lambda .  Next I'd like to show the proposition:

 \left( \lambda - C_{1,1}\right) \cdot \left( \lambda - C_{2,2} \right) \cdot \ldots \cdot \left( \lambda - C_{n,n} \right) = \sqrt[n-1]{\prod_{\forall i, \forall j, i \neq j}^n C_{i,j}}

and for this I want to begin with a combinatorial argument.  On the left hand side, the number of pairwise comparisons we can make depends on the number of  \lambda factors of the  \lambda polynomial (or, the highest degree of the  \lambda polynomial).  That is to say, we can make  \binom{n}{2} pairwise comparisons, or  \frac{n!}{(n-2)!2!} = \frac{n (n-1)}{2} comparisons.  Now, I don't know whether anyone has ever noticed this, but this last simplified part looks exceptionally like Gauss's sum of consecutive integers (pyramidal series), so in other words, this last part is in effect  \sum_{i=1}^{n-1} i which I find very cool, because we have just shown, quite accidentally, the equivalence:

 \binom{n}{2} = \binom{n}{n-2} = \sum_{i=1}^{n-1} i

The way I actually figured this out is by noticing that, in our pairwise comparisons, say for the 3rd-degree-polynomial-in- \lambda case, by writing the pairwise comparisons first of the  (\lambda - C_{1,1}) products, then of the  (\lambda - C_{2,2}) (in other words, ordering logically all  \binom{3}{2} products), there were 2 of the first and 1 of the second (and none of the  (\lambda - C_{3,3}) ).  If we do the same for the 4th-degree, there are 3 of the  (\lambda - C_{1,1}) , 2 of the  (\lambda - C_{2,2}) , and 1 of the  (\lambda - C_{3,3}) , with none of the  (\lambda - C_{4,4}) .  In other words, the  \binom{4}{2} pair-products could be written as the sum of the cardinality of the groupings:  3 + 2 + 1 .

Now Gauss's sum of integers formula is already known to work in the general case (just use an inductive proof, e.g.), so the substitution of it into the binomial equivalence needs no further elaboration: it generalizes automatically for all  n .

So if we are to multiply all pairwise comparisons, notice there will be  n - 1 products of each  \lambda -factor: there are  n - 1 products belonging to the  (\lambda - C_{1,1}) grouping (because this first grouping has n-1 entries, from the Gauss formula equivalence), there are  n - 2 products belonging to the  (\lambda - C_{2,2}) PLUS the one already counted in the  (\lambda - C_{1,1}) grouping, for a total of, again,  n - 1 .  The  kth grouping  (\lambda - C_{k,k}) has  n - k products listed for itself PLUS one for each of the previous k - 1 groupings, for a total of  n - k + k - 1 = n - 1, and the k+1th grouping  (\lambda - C_{k+1, k+1}) has  n - (k+1) products listed for itself PLUS one for each of the previous  k groupings, for a total of n - (k+1) + k = n - 1.  We are left in effect with:

 (\lambda - C_{1,1})^{n-1} \cdot(\lambda - C_{2,2})^{n-1} \cdot \ldots \cdot (\lambda - C_{n,n})^{n-1}

The right hand side of each pairwise comparison was nothing more than the simple product on the cross indexes of  C , so it's not difficult to argue then that, if we multiply  \binom{n}{2} such pairs, we get  \prod_{\forall i, \forall j, i \neq j}^n C_{i,j} .   We then take the  n-1 th root on both sides of the equation.

Since the  n + 1 case follows the same basic structure of the argument, we are done with proving our proposition.

What I want this proposition to mean may be very different than what it actually is, I'm hopeful nevertheless but I agree that it requires a bit of further investigation.  As I hinted before, I would like that

 \left( \lambda - C_{1,1}\right) \cdot \left( \lambda - C_{2,2} \right) \cdot \ldots \cdot \left( \lambda - C_{n,n} \right) = \sqrt[n-1]{\prod_{\forall i, \forall j, i \neq j}^n C_{i,j}}

with, for example,  n = 3 represent the constraint on the eigen(patch(ix))values of  a^\circ = B_1 f_1(x) + B_2 f_2(x) + B_3 f_3(x) or, if not that, maybe  a^\circ_\star = a(x) + a^*(x) + a^{**}(x) = 2B_1 f_1(x) + 2 B_2 f_2(x) + 2 B_3 f_3(x) , which brings into question the superposition of functions and their effect on the eigenvalues.  I may be wildly speculating, but hey!  I don't really know better!  I'll do a few experiments and see what shows.