nou

polynom coeff integ

2 mins read ·

Integration of polynomial: Coefficients operation

polynomial as list

A polynomial in the form of

y=i=0nbiti(1)\tag{1} y = \sum_{i=0}^n b_i t^i

can be presented in a column matrix

y[b0  b1  b2  b3    bn1  bn](2)\tag{2} y \equiv [b_0 \ \ b_1 \ \ b_2 \ \ b_3 \ \ \dots \ \ b_{n-1} \ \ b_n]

or simply a list in Python.

differential

Let us define other polynomial

x=j=0majtj(3)\tag{3} x = \sum_{j=0}^m a_j t^j

and use it to represent differential of (1)

dydt=i=1nibiti1.(4)\tag{4} \frac{dy}{dt} = \sum_{i=1}^n ib_i t^{i-1}.

Equate (3) and (4) will give

aj=(j+1)bj+1(5)\tag{5} a_j = (j+1) b_{j+1}

since

x=[b1  2b2  3b3    (n1)bn1    nbn ]=[a0   a1   a2           an2         an1].(6)\tag{6} \begin{array}{rcl} x & = & [b_1 \ \ 2 b_2 \ \ 3b_3 \ \ \dots \ \ (n-1) b_{n-1} \ \ \ \ n b_n \ ] \newline & = & [a_0 \ \ \ a_1 \ \ \ a_2 \ \ \ \dots \ \ \ \ \ \ \ \ a_{n-2} \ \ \ \ \ \ \ \ \ a_{n-1}]. \end{array}

Notice that in this case m=n1m = n - 1.

integral

Again another polynomial

z=j=0mcjtj(7)\tag{7} z = \sum_{j=0}^m c_j t^j

to represent integral of (1)

y dt=c0+i=0n1i+1biti+1.(8)\tag{8} \int y \ dt = c_0 + \sum_{i=0}^n \frac{1}{i+1} b_i t^{i+1}.

Equate (7) and (8) will give

cj=1j1bj1(9)\tag{9} c_j = \frac{1}{j-1} b_{j-1}

since

z=[c0  b0  12b1  13b2    (n1)bn1    nbn ]=[c0   c1   c2    c3           cn         cn+1].(10)\tag{10} \begin{array}{rcl} z & = & [c_0 \ \ b_0 \ \ \frac12 b_1 \ \ \frac13 b_2 \ \ \dots \ \ (n-1) b_{n-1} \ \ \ \ n b_n \ ] \newline & = & [c_0 \ \ \ c_1 \ \ \ c_2 \ \ \ \ c_3 \ \ \ \dots \ \ \ \ \ \ \ \ c_n \ \ \ \ \ \ \ \ \ c_{n+1}]. \end{array}

Tags: