Thread Rating:
  • 0 Vote(s) - 0 Average
  • 1
  • 2
  • 3
  • 4
  • 5
Observations on power series involving logarithmic singularities
#10
[updated]

jaydfox Wrote:Gottfried, I must admit that I haven't put enough study into your matrix methods, having focussed all my attention on Andrew's solution.

Can you point me to a discussion that mentions how to compute the coefficients of S2, P, etc., so that I can be sure I'm looking at the right matrices? If there is a relationship between your matrices and the ones I'm using, I'd like to understand it.

Jay,

I give a short list of the used matrices. I'm preparing a code-pad for Pari/gp with which you may then experiment with these matrices. Perhaps tomorrow or friday.

In general
V(x)~ = [1,x,x^2,x^3,...]
A prefix d declares this as diagonalmatrix
So, for instance
V(2)~ = [1,2,4,8,16,...]

F = [0!,1!,2!,...], dF arranges this as diagonal dF^-1 contains the reciprocals (to construct the exponential-series, for instance)

P (binomialmatrix)


Two properties of this matrix are important here.

1) Application of the binomial-rules, when postmutiplied by a formal powerseries


2) Derivatives

Y contains then scaled derivatives for



St2 Confusedtirling kind 2, version 1 (Abramowitsch&Stegun)


St2 Confusedtirling kind 2, version 2, Wikipedia


I use this version here




S2 : factorial scaled version of St2: dF^-1 * St2 * dF


This version performs U-exponentiation for a powerseries
V(x)~ * S2 = V(exp(x)-1)~
(see Abramowitsch & Stegun)

Since Input and output are of the form of a powerseries, one can iterate to implement the "x-> (exp(x)-1)"" - iteration







St1 : Stirlingnumbers 1'st kind = inverse of St2



S1 : factorial scaled St1, inverse of S2


V(x)~ * S1 = V(log(1+x))~

Since this is the inverse of S2, it performs x->log(1+x), and since Input and Output are of the form of a powerseries, this can be iterated
(See also A&S)






B : my base-matrix for T-iteration


This matrix can be understood in two ways:

As

just the same way as in your code-snippet

Or




Application is

and since input and output are of the form of powerseries, this can be iterated.

Note that since B = S2 * P~ we have, that
V(x)~ * S2 = V(e^x-1)~
and (see binomial-rules using P -transposed as above))
V(e^x-1)~ * P~ = V((e^x-1)+1)~ = V(e^x)~
we have
V(x)~ * S2 * P~ = V(e^x-1)~ *P~ = V(e^x)~
which is the same as
V(x)~ * B = V(e^x)~



To apply another base s, different from e, such that s<>e, this matrix must be premultiplied by powers of logarithms of s. I use lambda here for brevity. The decomposed description

Bs


Application is

and since input and output are of the form of powerseries, this can be iterated.




Bs numerically



The terms, that you use can -if at all- mostly be found in the second-column of the matrices or result, since you need only the final scalar result, and not the additional powers, which occur in the next coumns.

So much for short

Gottfried
Gottfried Helms, Kassel
Reply


Messages In This Thread
RE: Observations on power series involving logarithmic singularities - by Gottfried - 10/30/2007, 06:29 AM

Possibly Related Threads...
Thread Author Replies Views Last Post
  Perhaps a new series for log^0.5(x) Gottfried 0 83 12/05/2019, 04:35 PM
Last Post: Gottfried
  A Notation Question (raising the highest value in pow-tower to a different power) Micah 8 2,111 02/18/2019, 10:34 PM
Last Post: Micah
Question Taylor series of i[x] Xorter 12 11,088 02/20/2018, 09:55 PM
Last Post: Xorter
  Functional power Xorter 0 1,269 03/11/2017, 10:22 AM
Last Post: Xorter
  2 fixpoints related by power ? tommy1729 0 1,396 12/07/2016, 01:29 PM
Last Post: tommy1729
  Taylor series of cheta Xorter 13 11,965 08/28/2016, 08:52 PM
Last Post: sheldonison
  Inverse power tower functions tommy1729 0 1,808 01/04/2016, 12:03 PM
Last Post: tommy1729
  Further observations on fractional calc solution to tetration JmsNxn 13 13,443 06/05/2014, 08:54 PM
Last Post: tommy1729
  Remark on Gottfried's "problem with an infinite product" power tower variation tommy1729 4 5,017 05/06/2014, 09:47 PM
Last Post: tommy1729
  [integral] How to integrate a fourier series ? tommy1729 1 2,454 05/04/2014, 03:19 PM
Last Post: tommy1729



Users browsing this thread: 1 Guest(s)