Skip to main content

Latest Post

Power Essence Coefficients and Bernoulli Numbers

Previously, we have explored methods to compute the essence of power functions $Ȣx^n$ which involves solving a large system of linear equations. This method is equivalent to solving for the inverse of a large $n\times n$ matrix where entries are values of pascal's triangle. Though the matrix method allow us to solve for large number of essences at once, it does not extend easily to solve for next iterations of essence coefficients. Rather than reusing the values we have already solved for, we will have to solve for inverse of a separate larger matrix again. Here we will introduce an iterative method for solving for these coefficients. Chapter 0: Recap Let us first remind ourselves of the definition of essence. For a function $f(x)$, we want to find the transformation $Ȣf(x)$ such that we are able to 'smooth out' its series: $$\sum_{i=a}^b f(i) = \int_{a-1}^b Ȣf(x) dx$$ For example, we can solve for the following functions: $$\begin{align*}Ȣ1 &= 1 \\ Ȣx &= x +

Sequence Curve: Essence of Power Function with Induction Formula and Sequence Flatlining


Last time, we introduced the concept of Sequence Curve. This time we will explore how the Essence Transformation will actually work. 
To start this, we will consider the simplest function: linear function $f(x) = x$.
$$\sum_{i=n}^m i = \int_{n-1}^m Ȣ\{x\} dx$$

How can we find $Ȣ\{x\}$?

For this, we can utilize the well known Arithmetic Series formula to simplify this problem.
$$\sum_{i=n}^m i = \frac{m+1 - n}{2}(m+n) $$
This lends to that
$$\frac{m+1-n}{2}(m+n) = \int_{n-1}^m Ȣ\{x\} dx $$
From this relation we will be able to solve for $Ȣ\{x\}$.

Let’s consider the case when $n$ is a constant $c$. 
Then it must be case that
$$ \int_{c-1}^m Ȣ\{x\} dx = \frac{m+1-c}{2}(m+c) $$
From First Fundamental Theorem of Calculus, we find that
$$\begin{align*} \frac{d}{dm}\int_{c-1}^m Ȣ\{x\} dx & = \frac{d}{dm}\frac{m+1-c}{2}(m+c) \\
Ȣ\{m\} & = \frac{m+1-c}{2} + \frac{m+c}{2} \\
& = m + \frac{1}{2}
\end{align*}$$

Let's see this in action to demonstrate that this works:
$$\begin{align*} \sum_{i=3}^7 i & = \frac{7+1-3}{2}(7+3)  & = \int_{3-2}^7 x + \frac{1}{2} dx\\
3 + 4 + 5 + 6 + 7 & = \frac{5}{2}(10) & = \left[ \frac{x^2}{2} + \frac{x}{2} \right]^7_2 \\
& & = \frac{49}{2} + \frac{7}{2} - \frac{4}{2} - \frac{2}{2}\\
25 &= 25 &= 25 \end{align*}$$

We can continue pairing Induction with Essences, but that will limit us to get only as far as Inductions we already know. What we are trying to achieve through Essence Transformation is to find a systematic approach of rather than relying on Induction Theory. 

Sequence Flatlining


Let's try to solve for $Ȣ\{x\}$ in a different method that does not involve using induction formulas.
Consider the case when $m=n$:
$$\sum_{i=n}^m i = \sum_{i=n}^n i = n$$
Our Essence Transformation should still hold true even for these cases, leading us to
$$\sum_{i=n}^n i = n = \int_{n-1}^n Ȣ\{x\} dx$$
For the sake of ease, let's define $Ȣ\{x\} = f(x)$.
Like with the previous method, we apply FTC and find that
$$\begin{align*}
n = & \int_{n-1}^n f(x) dx\\
1 = & f(n)-f(n-1) \\
0 = & f'(n) - f'(n-1)
\end{align*}$$
At this point I want to introduce a lemma:

Naught Curve Difference Postulate

If the difference of function $f$ between $x$ and $x+c$ is 0, then it must be that $f$ is constant or that $f$ oscillates with period of $c$.
Given that postulate, and since there is no reason to believe that $f'$ would be oscillating, we can assume that $f'(x) = c$. (We will discuss how to solve for $Ȣ\{x\}$ without making such assumptions later on, but for now know that this assumption still gives us consistent results)
Since, $f'(x)$ is a constant, it must follow that
$$f'(x) = c_1\\
f(x) = c_1x + c_2\\
\int f(x) dx = \frac{c_1}{2}x^2 + c_2x + c_3$$
Let us now substitute these into the identities above, we find that
$$\begin{align*} 1 = & f(n) - f(n-1)\\
= & c_1n + c_2 - c_1(n-1) - c_2\\
= & c_1n - c_1n + c_1\\
= & c_1 \end{align*}$$ And then $$\begin{align*} n = & \int_{n-1}^nf(x)dx\\
= & \int_{n-1}^n x+c_2 dx\\
= & \frac{1}{2}n^2 + c_2n - \frac{1}{2}(n-1)^2 - c_2(n-1)\\
= & \frac{1}{2}(n^2 - (n-1)^2) +  c_2n - c_2(n-1)\\
= & \frac{1}{2}(2n - 1) + c_2\\
n = & n - \frac{1}{2} + c_2\\
\frac{1}{2} = & c_2
 \end{align*}$$
Combing them together we find that
$$f(x) = Ȣ\{x\} = x+\frac{1}{2}$$ which is the same result we got from deriving $Ȣ\{x\}$ from Arithmetic Series formula. $\blacksquare$

Notice that this method relies differentiating the sequence formula until it eventually reaches 0, at which point we apply Naught Difference Postulate and substitute values backwards to solve for Essence. This will work for all Integer Power Series.
We can think of this process as 'flatlining' the sequence; as morbid as it sounds, it captures the spirit of the technique pretty well.

Let's use this method on next simplest Sequence $x^2$.
$$\begin{align*} \sum_{i=n}^ni^2 =  n^2 = & \int_{n-1}^nf(x)dx\\
 2n = & f(n) - f(n-1)\\
 2 = & f'(n) - f'(n-1)\\
 0 = & f''(n) - f''(n-1)\\
& \therefore f''(x) = c_1 \text{, by Naught Difference Postulate}\\ \\
& \therefore f'(x) = c_1x + c_2\\
\therefore 2 = & c_1n + c_2 - c_1(n-1) - c_2\\
2 =& c_1\\ \\
& \therefore f(x) = x^2 + c_2x + c_3\\
\therefore 2n = & n^2 + c_2n + c_3 - (n-1)^2 - c_2(n-1) - c_3\\
= & 2n - 1 + c_2\\
1 =& c_2 \\ \\
& \therefore F(x) = \frac{1}{3}x^3 + \frac{1}{2}x^2 + c_3x + c_4\\
\therefore n^2 =& \frac{1}{3}n^3 + \frac{1}{2}n^2 + c_3n - \frac{1}{3}(n-1)^3 - \frac{1}{2}(n-1)^2 - c_3(n-1) \\
= & \frac{1}{3}(n^3 - (n-1)^3) + \frac{1}{2}(n^2 - (n-1)^2) + c_3 \\
= & \frac{1}{3}(3n^2 - 3n + 1) + \frac{1}{2}(2n -1) + c_3 \\
n^2 = & n^2 - n + \frac{1}{2} + n - \frac{1}{2} + c_3\\
\frac{1}{2} - \frac{1}{3} = & c_3\\
\frac{1}{6} = & \\ \\
& \therefore f(x) = Ȣ\{x^2\} = x^2 + x + \frac{1}{6}\\ \blacksquare
\end{align*}$$

Surprisingly, this result, once you integrate it, is also consistent with known induction formula $\frac{n(n+1)(2n+1)}{6}$.

We can continue on with Sequence Flatlining technique and discover that $$\begin{align*} Ȣx^3 = & x^3+\frac{3x^2}{2}+\frac{x}{2} \\
Ȣx^4 = & x^4 + 2x^3 + x^2 - \frac{1}{30}\\
Ȣx^5 = & x^5 + \frac{5x^4}{2} + \frac{5x^3}{3} - \frac{x}{6}\\
Ȣx^6 = & x^6 + 3x^5 + \frac{5x^4}{2} - \frac{x^2}{2}+\frac{1}{42}\\
Ȣx^7 = & x^7 + \frac{7x^6}{2} + \frac{7x^5}{2}-\frac{7x^3}{6}+\frac{x}{6}\\
Ȣx^8 = &x^8 + 4x^7 + \frac{14x^6}{3}- \frac{7x^4}{3}+\frac{2x^2}{3} - \frac{1}{30} \\
Ȣx^9 = &x^9 + \frac{9x^8}{2}+ 6x^7- \frac{21x^5}{5}+ 2x^3 - \frac{3x}{10}\\
Ȣx^{10} = &x^{10} + 5x^9 +\frac{15x^8}{2}-7x^6+ 5x^4-\frac{3x^2}{2}+\frac{5}{66}
\end{align*}$$ and so on and so on.

Again, let's demonstrate that Essence Transformation works as intended:
$$\begin{align*} \sum_{i=11}^{18}i^5 = & 11^5 + 12^5 + 13^5 + 14^5 + 15^5 + 16^5 + 17^5 + 18^5 \\
& \vdots \\
= & 6,436,376 \\ \\
\int_{11-1}^{18} x^5 + \frac{5x^4}{2} + \frac{5x^3}{3} - \frac{x}{6} dx = & \left[  \frac{x^6}{6} + \frac{x^5}{2} + \frac{5x^4}{12}-\frac{x^2}{12} \right]_{10}^{18}\\
= & (\frac{18^6}{6} + \frac{18^5}{2} + \frac{5*18^4}{12}-\frac{18^2}{12})- (\frac{10^6}{6} + \frac{10^5}{2} + \frac{5*10^4}{12}-\frac{10^2}{12})\\
= & 6,657,201 - 220,825\\
= & 6,436,376
\end{align*}$$
Pretty cool, right?
Notice that as the range gets bigger, calculating the series through Essence would be much simpler and easier than summing up individual terms of the sequence.

Note, also, that through Essence definition of series, something like $\sum_{i=1}^{1.5}i$ becomes computable. This was not possible from original definition of series as incrementing by a unit from lower bound $1$ will never arrive us at the upper bound $1.5$.
But using our new definition of series we find that
$$\begin{align*}\sum_{i=1}^{1.5}i = & \int_{0}^{1.5}x+\frac{1}{2} dx \\
 = & \frac{15}{8} \end{align*}$$

One interesting point to note is that for Power Sequences $Ȣf(x)$ results in curve with same degree as $f(x)$. This means that solving for Essence of Power Sequence can essentially be reduced to finding patterns in each terms' coefficients.
The general expression for the coefficients will be studied in further depth in future articles.

Comments

Popular posts from this blog

Large Polynomial Series using Matrices (Calculating Bernoulli's Number with Pascal Matrix)

Polynomial Series can be easily solved using Power Series Formulas for each term in the polynomial. However, this can be frustrating since not every Power Formula are intuitive to memorize. We would like to find a more elegant and easier-to-recall formula for computing a Polynomial Series. This can be done using matrices. Notations and Equations We will borrow the notations and simple equations from Sequence Curve articles . There, we have extended a series to be continuous through the following identity: $$\sum_{i=m}^n f(i) = \int_{m-1}^nȢ\{f(x)\}dx $$ The $Ȣ\{f(x)\}$ acts as the rate of change of a series relative to its bounds. $$\frac{d}{dt} \sum_{i=m}^nf(i) = Ȣ\{f(n)\}\frac{dn}{dt} - Ȣ\{f(m-1)\}\frac{dm}{dt}$$ Letting $t=m=n$, we find $$\frac{d}{dt} \sum_{i=t}^tf(i) = \frac{d}{dt}f(t)= Ȣ\{f(t)\} - Ȣ\{f(t-1)\} $$ Remebering that $Ȣ$ Transformation is Linear, we can derive a simple identity of $$\frac{d}{dx}f(x+1) = Ȣ\{f(x+1)\} - Ȣ\{f(x)\} = Ȣ\{f(x+1)-f(x)\}$$ This will be use

Partition Counter using Trees, Recursion, Tables, and Algorithm

Partitions are number of ways an integer can be represented as sum of positive integers. We can ask, for example, what is the partition of 5? If we write out every possible combination, $$\begin{align*}  5 &= 1+1+1+1+1 \\  &= 1+1+1+2\\ &= 1+1+3\\ &= 1+4\\ &= 1+2+2\\ &= 2+3\\  &= 5 \end{align*} $$ we can see that partition of 5 is 7. One will immediately notice, however, that this is not the most efficient approach to answering the question: not only does partition grow quickly, attempting to organize and not miss or repeat an algebraic expression becomes messy and impractical. Chapter 1: Partition Tree A cleaner, but still not the best, approach would be to use tree diagrams to represent the partitions of a number. In a Partition Tree of $n$, every path to a terminating node will add up to the number $n$. And also, every child of a node (every nodes below a given parent node) will always be greater than or equal to the parent node in order to not

Vector Plane Rotation

Vector is a useful tool in maths and in physics that describes both magnitude and direction in a space. Vectors are used to describe concepts where direction matter, such as forces and velocity. For example, two cars can have same speed(scalar) but if one if heading north and the other south, they will have different velocity(vector). Vectors are usually written in terms of their components. Components are magnitude of the vector in one dimension. In 2D space, vector can be represented as coordinate of 2 components, <x,y>, and in 3D space, vector can be represented as coordinate of 3 components, <x, y, z>. Breaking down vectors into their components is extremely useful, especially in physics, because it allows one to break down complicated changes in vector into simpler and easier-to-work-with changes in its components. Take a look at projectile motion, for example. Velocity of an object fired up into the sky at an angle changes both its direction and magnitude. Attem

Summation of Harmonic Interval and Simplified Riemann Sum

When you first learn of Summation in Pre-Calc class, it is probably in build-up to Riemann Sum. You are taught the basics of Summation and simple manipulations, such a $\sum_nc*a_n = c*\sum_na_a$ for constant $c$, but you don't get to explore deeper into Summations and really play around with it. One thing that really intrigued me was why Summations only incremented by a unit. When we have a Sum like $\sum_{n=3}^{10}a_n$, it is equal to $a_3+a_4+a_5+\dots+a_9+a_{10}$, always incrementing by a unit( $1$ ). Why couldn't we Sum with interval that is not a unit? How could we have something like $a_3+a_{3.5}+a_4+a_{4.5}+a_5+a_{5.5}+\dots+a_{8.5}+a_9+a_{9.5}+a_{10}$, incrementing by a half instead of a full unit? Sets Of course you technically can use Sets and Set Builder Notation to create a specific Summation. $$\sum_{i \in S}i,S=\{\frac{x}{2}|x\in\mathbb{N_0}\} \, = \, 0+\frac{1}{2}+1+\frac{3}{2}+\dots$$ The above is Sum of all Rational Numbers incrementing by $\frac{1