The Legacy of Peter Wynn

: After the death of Peter Wynn in December 2017, manuscript documents he left came to our knowledge. They concern continued fractions, rational (Padé) approximation, Thiele interpolation, orthogonal polynomials, moment problems, series, and abstract algebra. The purpose of this paper is to analyze them and to make them available to the mathematical community. Some of them are in quite good shape, almost ﬁnished, and ready to be published by anyone willing to check and complete them. Others are rough notes, and need to be reworked. Anyway, we think that these works are valuable additions to the literature on these topics and that they cannot be left unknown since they contain ideas that were never exploited. They can lead to new research and results. Two unpublished papers are also mentioned here for the ﬁrst time.


Introduction
Peter Wynn (1931Wynn ( -2017 was a mathematician, a numerical analyst, and a computer scientist (see in Figure 1 a photo of him taken in 1975). In his scientific life he produced 109 publications (see References ), and he translated two books from Russian [110,111]. He is mostly known for his discovery of the ε-algorithm [3], a recursive method for the implementation of the Shanks transformation for scalar sequences [112], for its extensions to the vector, matrix, and confluent cases [12,13,23,24], and for his numerous reports and papers on Padé approximants and continued fractions. His works influenced a generation of pure and numerical analysts, with an important impact on the creation of new methods for the acceleration of scalar, vector, matrix, and tensor sequences, on the approximation of functions, and on iterative procedures for the solution of fixed point problems. Volume 80,No. 1 of the journal Numerical Algorithms was dedicated to him with his full biography. More recently, a complete analysis of all his works was provided in Reference [113], together with those of other scientists who worked and are still working on these domains. Thus, one can wonder why it was necessary to publish an additional paper on Wynn's work. During the last years of his life, Peter Wynn was living in Zacatecas, Mexico. Each year, he had to come back to the United States for some administrative reasons. On one occasion, he was visiting friends in San Antonio and left them boxes containing mathematical documents he did not want, for some unknown reason, to keep with him in Mexico. In January 2020, C.B. was contacted by F.A.N., a colleague of these friends, who informed him of the existence of these documents. This is how Wynn's legacy came to light. Then, the authors of this paper decided to analyze these unpublished works.
As everyone can understand, it was a quite difficult task. Only a part of the documents have been extracted from the boxes and studied. The handwritten lists made by Wynn for indicating the contents of the boxes show that he put together several kinds of documents. What he named "rough notes" are very difficult to read and understand. In these lists, he often indicated what he called "notes" and, in this case, they are usually well written and understandable. The best documents are designated as "paper". Moreover, Wynn sometimes wrote new notes on the back of another document! In addition, he also made Xerox copies of documents, and inserted them into the boxes. However, his lists help us to try to identify the kind of document we were considering. But, sometimes, this was quite impossible. Moreover, most of the documents we found have no date and, often, pages are not numbered. When unnumbered sheets of paper are stacked on top of each other, without any separation, it has sometimes been difficult to know where a document begins and where it ends. It is also possible that Wynn himself mixed up some texts. For these reasons, certain groupings of pages may be questionable. Thus, we apologize in advance for all possible mistakes contained in this paper.
Let us motivate the potential reader by giving an idea of the main themes covered by Wynn. The most important documents left by Wynn, which are almost complete and in a good shape are the following. One of them concerns Bürmann series over a field; they generalize Taylor series and are used in the reversion of series. Another one (187 pages) is on stability and F-functions that play a role in the solution of the differential equation y (t) = Ay(t), variation diminishing functions, interpolating rational functions, exponential fitting forms. The Hamburger-Pick-Nevanlinna problem is treated in a document of 179 pages. There is a document on continued fraction transformations of the Euler-Mclaurin series that has 202 pages; it contains applications to various series. The convergence of associated continued fractions, and truncation error bounds for Thiele's continued fractions are the topics of another document. Then, we analyze various documents on interpolation.
The first one is about functional interpolation, in which a recursive algorithm, which seems to be new, is given for constructing interpolating rational functions. Interpolation by the use of rational functions is studied in another document. Wynn gave two recursive algorithms for their computation. They also seem to be new. A document extends a report of Wynn on the abstract theory of the ε-algorithm [69]. The ε-algorithm is applied to sequences of elements of a ring. An interesting pedagogical document is on iterated complex number spirals. We present some numerical experiments illustrating Wynn's ideas. There is a document that looks like a book project on extrapolation, Padé approximation, continued fractions, and orthogonal polynomials. It can certainly serve as a basis for lectures on continued fractions, Padé approximation, and the ε-algorithm. Other documents are on algebra. One of them (266 pages) is on S-rings. It is quite theoretical, without any application nor reference to the literature. A second one treats factorisations of a triangular matrix. There are also unfinished manuscripts on various topics, which can be of interest. Two unpublished papers of Wynn are also mentioned here for the first time. The documents left by Wynn show the intellectual process leading to the elaboration of new results until their publication. They are also a testimony on the human side of research by describing the friendship and the collaboration between researchers, and their mutual influence.
The purpose of this paper is to make this legacy available to the international mathematical community. It contains a description of the unpublished manuscripts of Wynn. They offer many new results and developments. Despite the fact that not all documents have been sorted yet, we decided to propose them immediately since we think that the research must go on. We hope to encourage some readers to resume the work of Wynn and bring it to an end. We are sure that several of his ideas are worth pursuing.
All the documents extracted up to now from the boxes left by Wynn have been digitized, and they can be downloaded from the following URL (Legacy Archive: Peter Wynn). The main material is mentioned in this paper. The unusable or incomplete documents we found are not listed here, but they are also inserted in the same website together with a small description. Other information on Peter Wynn can also be found in the site Mathematics Research of the Department of Mathematics at the University of Texas at San Antonio (UTSA) https://mathresearch.utsa.edu/Legacy/Peter-Wynn/ (accessed on 30 April 2021).
The history of the discovery of these documents is told in Section 2 with the testimony of F.A.N. and of the friends of Wynn who inherited them. For readers who are not familiar with the topics touched upon by Wynn, a short mathematical introduction is provided in Section 3. More details can easily be found in the literature, particularly in Reference [113].
Some of the documents left by Wynn in San Antonio are analyzed and commented in Section 4. Some others will be analyzed in a second paper if they are of interest.
All quotations from Wynn are in italics. Inside a quotation, our own comments (when necessary) are in roman characters into square brackets. When Wynn mentions a reference given at the end of his text, we replaced it by that of our own bibliography numbered in arabic figures between square brackets [·]. The documents of Wynn are numbered in bold italics arabic figures. For referring to them in the text, these numbers are placed between curly brackets {·} to distinguish them from the bibliographical references. Concerning the references of Wynn, we were able to update them, by inserting new DOI, MR (Mathematical Reviews) and Zbl (Zentralblatt reviews) numbers, and by adding newly discovered papers. Let us also mention that in 1960, the journal Mathematical Tables and Other Aids to Computation (known in short as MTAC) changed its name and became Mathematics of Computation, and that from its Volume 5, the journal Revue française de traitement de l'information, Chiffres became simply Chiffres. Several papers and communications of Wynn, when he was at the Stichting Mathematisch Centrum (now Centrum Wiskunde & Informatica) in Amsterdam, from 1960 to 1964, can also be found at the CWI's Institutional Repository at the address https://ir.cwi.nl/ (accessed on 30 April 2021).

The Discovery
On 14 January 2020, C.B. received the following message from F.A.N.: I have just now read your recent remembrance of Peter Wynn appearing in Numerical Algorithms. I was charmed-thank you for that. But I am also very much interested because, serendipitously, I have just today "inherited" from a retired colleague several rather heavy boxes of Wynn's papers which had been left with my colleague for safekeeping (or convenience) some years earlier. I have not yet opened the boxes, but I will begin examining the contents soon. The reason I am writing you is to ask if you would be willing to answer a couple of questions, which I pose below.
As I begin looking at the material, I could find (1) drafts of papers that might be of interest to mathematicians, (2) personal items that might be of interest to family or friends, (3) miscellany that would be of little interest to anyone, or (4) items for which I am at a loss to know whether to keep or discard [. . . ] Here are my questions to you: In case I find material that looks as if it might have some mathematical import, would you be willing to take look at it, or could you suggest someone who would be?
The colleague he was speaking about is Manuel Philip Berriozábal. He was born in 1931. He was awarded the Bachelor of Science degree in mathematics from Rockhurst College in 1952, a Master of Science degree in mathematics from the University of Notre Dame in 1956, and a Ph.D. in mathematics from the University of California at Los Angeles (UCLA) in 1961. After serving for one year as a lecturer at UCLA, he joined the faculty at Tulane University as an Assistant Professor. Four years later, he moved to the University of New Orleans as an Associate Professor. He was promoted to Professor six years later. In 1975, Manuel Berriozábal married Maria Antonietta Rodriguez (see Figure 2). He joined the faculty at the University of Texas at San Antonio (UTSA) in 1976, and in 1979, he started the now nationally recognized Prefreshman Engineering Program (PREP) at UTSA. San Antonio PREP received a Presidential Award for Excellence in Science, Mathematics and Engineering Mentoring and a La Promesa Program Award from the National Latino Children's Institute. Several years ago, TexPREP (Texas Prefreshman Engineering Program) received a special commendation from the Texas Senate. These accomplishments caught the attention of the Washington, DC-based Quality Education for Minorities Network and resulted in Berriozábal being named one of the six Giants in Science at a conference held in February 1998. In May 1998, he was a recipient of the San Antonio "I Have a Dream" Foundation Endeavors Award. PREP has also been replicated on eight college campuses in eight states outside of Texas. It was during his professorship at the Louisiana State University in New Orleans that he became a close friend of Peter Wynn. Of course, C.B. and M.R.-Z. were very much interested in the contents of these boxes since they were in the final preparation of their book which contains, among others, an analysis of all the publications of Peter Wynn [113]. F.A.N. began to look at them, to scan a lot of their material, and to send it to us. The correspondence with him was completely taken up by M.R.-Z., since the same day C.B. received F.A.N.'s first message, he had to enter a hospital for a health problem.
On 15 January, F.A.N. wrote to C.B. and M.R.-Z.: I have opened two of the four boxes, looking thoroughly through the contents of one of them. FIRST BOX The first box was a smaller one containing several hundreds of pages of handwritten notes and a few other random items. In this and the next message, I have attached several of the smaller items that might be indicative of material to come. Nothing is dated so I have no idea whether this material simply anticipates papers that have been completed and published or represents new work. In addition, not much is paginated, so I generally presume that the order in which pages appear is the order in which they were completed and assembled. Some of this material was left with Dr. Berriozábal here in San Antonio after Peter had visited Mexico and there are some suggestions (e.g., written accounts of receipts) that suggest that some of the material was produced while he was living in Mexico. Whether this was after his last publication in 1981, I do not know.
To Do List: The first document above appears to be a list of various topics he wished to address, mainly through the construction of relevant notes.
Projects: This document lists a number of projects that he had planned and may well have completed. Whether these "projects" have found their way into the literature as published papers, I do not know. I've copied a snapshot from the document and you will note that some of the items have been denoted with a "D__" . Perhaps this is Peter's shorthand for "Done" or fait accompli. I didn't see in his publications any that specifically included, for example, "stratified commutative ring" or "Bürmann series", so perhaps this represents some work that hasn't yet been published, but will appear somewhere in his collection of notes. I did find the beginnings of what looks to be a monograph on Interpolation Theory-the first in his list of projects and one that doesn't seem to have been completed at that time.
Duplicate List: This may reflect a list of some of the items that Peter included in his boxes of documents.
The next messages will include 2 larger documents, one a bibliography and the other a more extensive list of documents.
Other items in this first box include the aforementioned piece of a monograph on Interpolation Theory, as well as a couple of hundred pages of paginated notes of another document. In that case, I see pages 70-300+ so I can't be sure of the initial title. There is another paginated treatment of a topic (I don't recall exactly what at the moment), a book review, and some hundreds of other pages of notes that do not appear to fit together well. SECOND BOX I'll mention here that the second box, one that originally held a dozen reams of office paper, now holds about 20 kg of handwritten notes. I have not yet examined these notes in detail, and when I do, I want to be careful not to disturb the assembly of any documents.
Bibliography: This is likely not anything that you would use. I imagine Peter kept this bibliography simply for reference purposes.
List of documents: This seems to be a list of documents and/or things he felt he needed to do, such as copy this or Xerox that. I don't know if it has any intrinsic value, but it might reflect documents included among his papers.
In a message to F.A.N., 17 January 2020, Maria Antonietta Berriozábal wrote: We liked his visits. Each time he came it was without much notice and we would just go out to eat. For me I was so very busy for so many years and so was Manny [her husband Manuel] that a surprise visit took creativity to get the three of us together, but we did it for Peter. He liked nice places. And I loved hearing he and Manny talk old times. He would also talk about living in Mexico which he loved.
He was such a mystery to me because he was pretty much alone and he liked it that way. I do not recall meeting Peter when we were in New Orleans which was August 1975 to June 1976 [. . .] What I do recall is that after we married I was sorting and clearing papers for Manny from our apartment which he had lived in for many years and I came across a letter from Peter written around the time Manny and I met. In that letter I gathered Peter was a bachelor friend of Manny who was not the marrying kind and he was commenting that Manny had met "the one". He was happy for Manny but it seemed to me Peter was happy being alone. Now that I think of it I don't think I ever got to know who Peter was. He had a wonderful smile. He liked to laugh and joke and was very very blond.
C.B. and M.R.-Z. wrote to Maria Berriozábal and sent her the photo of Wynn inserted above (also see References [113,114]). She answered on 21 January 2020: [. . .] The photos are valuable. I met Peter when he was a much older man but still had that incredibly beautiful smile and lots and lots of hair! Each time he came to the US for many years, he would call Manny and we would go to dinner. If I recall correctly, it was mostly on Peter's nickel and as I stated to Sandy [F.A.N.] we went to nice restaurants, although Manny may have won some times and we would go to our famous Luby's Cafeteria-always the frugal Manny.
I do not recall Peter in New Orleans when I lived there with Manny the first nine months of our marriage. I left San Antonio to join him as he was waiting and hoping to join UTSA. The New Orleans scholars community from the four universities and colleges there met socially and regularly at people's homes but I do not recall Peter in any of that circle. It could have been that he was no longer in New Orleans. This was August 1975 to June 1976.
This morning when I read your note the words that rang in my head is how youth is wasted on the young. I wish I had paid more attention and had retained the conversations with Peter, but at that time my life was incredibly full and so was Manny's. There were times when I had to make major changes in my schedule to join Manny and Peter and I always did and looked forward to those dinners. My questions to Peter were always about Mexico and his life there. That is what I wish I could recall more. I just know that he was very happy there. I do not recall what his conversations with Manny were all about and now Manny is forgetting a lot of things.
On one of these trips which had to occur every six months, I believe, because that is how he kept his US citizenship active he said he had some boxes with him and asked if we could store them for him for safekeeping. I recall his taking them out of his car. When I looked at them again seriously in recent years I saw how big the boxes were and soon learned how heavy they were. Yet, I seem to recall Peter carrying them by himself. Maybe a Dolly helped. I just know he trusted us with them.
One of my last year's resolutions was to clean the garage and I did. It took some weeks but with some help I did it. Peter's boxes were the only things that were still in there without any resolution. We had been telling each other for a couple of years that Peter had probably died since we did not hear from him anymore. He quit coming. Manny and I would comment from time to time that we still had those boxes and that we needed to do something about them.
A year ago when the boxes were the last unsettled things in the garage I told our gardener who was helping me with the project the story of the boxes. He said what if it is money. Since it had been so long that Peter had left the boxes with us I gave myself permission to open them. Hilario, my helper, even wondered if they were full of money! I only opened two boxes and all I saw were reams of papers with math problems. Pages and pages and pages. When Sandy [F.A.N.] helped us with Manny's papers after he retired I decided to ask him to go over the papers and determine what they were. Manny looked at them but had no idea what the material was. When Sandy [F.A.N.] took the boxes out I saw how at least one of them had folders of some kind so the papers were not only sheets of math problems. I hope someone who has Peter's interest and his career uppermost in their mind will continue to review it and possibly record it for posterity.
To close, thank you again for honoring Peter's work.

Mathematical Background
Some documents written by Wynn are pure algebra or complex analysis. They require, at least for some of them, quite a good knowledge of these fields. However, it is impossible to give herein a full account of the definitions and notions necessary for their complete understanding. We will restrict ourselves to the most specialized topics addressed by Wynn in the domain of numerical analysis, namely the Shanks transformation and the ε-algorithms, Padé approximation, continued fractions, and rational interpolation. In this section, we only present the definitions and the main algebraic properties that are sufficient for understanding most of the documents analyzed. The fundamental question of convergence is not addressed. We refer to References [113,[115][116][117][118][119] for more details.

The Shanks Transformation and the ε-Algorithms
When a sequence (S n ) of numbers is slowly converging to its limit S, and when one has no access to the process building it, it can be transformed into a new sequence (or a set of new sequences), which, under some assumptions, converge(s) faster to the same limit. Many such sequence transformations exist and have been studied; see References [115,[117][118][119]. Among them, one of the most well known, studied, and used is due to Shanks [112]. It consists in transforming (S n ) into a set of sequences denoted {e k (S n )}, indexed by k and n, and defined by These numbers can be recursively computed by the scalar ε-algorithm of Wynn [3], whose rule is ε with ∀n, ε 2k+1 are only intermediate results with no interest for our purpose. In fact, quantities with a different parity can be eliminated from (2), thus leading to the cross-rule of Wynn [48] which only links those of the same parity It can be proved that ∀n, ε (n) 2k = S if and only if the sequence (S n ) satisfies a linear difference equation of order k a 0 (S n − S) + a 1 (S n+1 − S) + · · · + a k (S n+k − S) = 0, n = 0, 1, . . . , where the coefficients a i are such that a 0 a k = 0 and a 0 + · · · + a k = 0. In other words, if and only if S n has the form where A i , B i and C i are polynomials in n such that, if d i is the degree of A i plus 1 for i = 1, . . . , p, and the maximum of the degrees of B i and C i plus 1 for i = p + 1, . . . , q, one has with the conventions that the second sum vanishes if there are no complex zeros, and m = −1 if there is no term in δ in (Kronecker's symbol). The set of such sequences is named the kernel of the transformation. Since many sequences produced by iterative methods have this form or are close to it, it explains the success of this algorithm.
In Reference [24], Wynn extended his algorithm to the case where the S n are vectors or square matrices. In the matrix case, the significance of the power −1 in (2) is clear. For vectors, the inverse y −1 of a vector y is defined as y −1 = y/(y, y), thus leading to the vector ε-algorithm.
Similarly, when a function f (t) is slowly converging to its limit S when t tends to infinity, it can be transformed into a set of functions converging faster to S under certain assumptions. For that purpose, Wynn extended his algorithm to that case by proposing the first confluent form of the ε-algorithm whose rule is, for all t, ε k+1 (t) = ε k−1 (t) + (ε k (t)) −1 , k = 0, 1, . . . , with ε −1 (t) = 0 and ε 0 (t) = f (t). Again the ε 2k+1 (t) are intermediate computation. It can be proved that, for all t, ε 2k (t) = S if and only if f satisfies the differential equation of order k a 0 ( f (t) − S) + a 1 f (t) + · · · + a k f (k) (t) = 0, with a 0 a k = 0, that is, in other words, where A i , B i and C i are polynomials in t such that, if d i is the degree of A i plus 1 for i = 1, . . . , p, and the maximum of the degrees of B i and C i plus 1 for i = p + 1, . . . , q, one has Moreover, the ε 2k (t) can be expressed by a ratio of determinants quite similar to (1), but in which the derivatives of f are replacing the powers of the difference operator ∆ (see, for example, Reference [113] (p. 24) or Reference [115] (p. 257)).

Padé Approximation
Let f be a formal power series in which the coefficients c i and the variable t can be complex. A Padé approximant of f is a rational function with a numerator of degree p at most and a denominator of degree q at most such that its power series expansion agrees with f as far as possible, that is up to the degree p + q inclusively. Such an approximant is denoted [p/q] f and, by construction, it holds Then, linearizing the conditions of the definition as with the convention that c i = 0 for i < 0 which allows to treat simultaneously the cases p ≤ q and p ≥ q. Setting b 0 = 1 allows to solve the system formed by the preceding last q equations for the coefficients b 1 , . . . , b q . Knowing the b i 's, the first p + 1 equations directly provide the a i 's.
It holds where f m denotes the partial sum of f up to the term of degree m inclusively. It is easy to see from (1) that, applying the ε-algorithm to S n = ∑ n i=0 c i t i leads to ε (n) 2k = [n + k/k] f (t). Let g be the reciprocal series of f defined by f (t)g(t) = 1 (it exists if and only if c 0 = 0). If the ε-algorithm is applied to the sequence of the partial sums of g, then ε

Continued Fractions
A continued fraction is an expression of the form . For evident typographical reasons, it is written as The numbers a k and b k are called the kth partial numerator and partial denominator, respectively, a k /b k is the kth partial quotient, and C n = b 0 + a 1 b 1 + · · · + a n−1 b n−1 + a n b n is called the nth convergent of the continued fraction C (even if the sequence (C n ) does not converge). A continued fraction is said to converge if the sequence (C n ) converges when n tends to infinity. After reducing to the same denominator, C n can be written as C n = A n /B n . It can be computed by the recurrence relationships Let (C n ) be the sequence of convergents of the continued fraction C and let (C p n ) be a subsequence. The continued fraction whose convergents satisfy C n = C p n is given by a n = C p n−1 − C p n C p n−1 − C p n−2 , b n = C p n − C p n−2 C p n−1 − C p n−2 , with b 0 = C p 0 , b 1 = 1 and a 1 = C p 1 − C p 0 . This operation is called a contraction of the continued fraction. Usually, p n = 2n. The analytic theory of continued fractions is concerned with continued fractions whose partial numerators and/or denominators are functions of the complex variable z. Let us consider the continued fraction From the recurrence relations, we see that A 2k−1 , A 2k and B 2k are polynomials of degree k in z and that B 2k−1 is a polynomial of degree k − 1. The expansions of C k and C k−1 in ascending powers of z agree up to the term of degree k − 1 inclusively. It is possible to choose b 0 , a 1 , a 2 , . . . so that the expansion of C k agrees with that of a given series f (z) = c 0 + c 1 z + c 2 z 2 + · · · up to the term of degree k. This continued fraction is called the continued fraction corresponding to the series f . By a contraction of this continued fraction, with p k = 2k as explained above, we obtain a continued fraction whose convergent C k agrees with that of f up to the term of degree 2k. This is the continued fraction associated to the series f . Thus, by the uniqueness property of Padé approximants, we have

Rational Interpolation
Consider the continued fraction k−2 for k = 1, 2, . . ., and α 2k can be expressed as a ratio of two determinants.
The kth convergent C (n) k of this continued fraction satisfies the interpolation conditions C k+1 are the reciprocal differences of f . They formed the -algorithm used by Wynn for rational extrapolation at infinity [2].
Let us now consider the confluent reciprocal differences of a function f defined by with −1 (t) = 0 and 0 (t) = f (t). This formula will be used as an extrapolation method for functions by Wynn, and named the confluent form of the -algorithm [12]. Thiele's expansion of a function f consists in the continued fraction with α k (t) = k (t) − k−2 (t) for k = 1, 2, . . . Replacing t by 0 and h by x, we get The successive convergents C k (x) = A k (x)/B k (x) of this continued fraction are such that f (x) − C k (x) = O(x k+1 ). Since Padé approximants are uniquely defined, it holds

The Legacy
Let us now describe the various documents contained in the boxes left by Wynn at his friends' house in San Antonio. All, except two, are handwritten.

Main Documents
We tried to sort the documents by themes. However, our classification is only an attempt since there are many connections between the topics, the documents are not dated, and, maybe, some of them contain pages coming out from various sources since they are not numbered and, maybe, inserted by Wynn in disorder into the boxes.

1
Bürmann series over a field These titled notes (52 pages) are devoted to Bürmann series. They are mentioned in documents {26-28}, and were probably written at the same time as {14}, since both titles are listed together.
The Bürmann series of a function f is a generalized form of a Taylor series in which, instead of a series in powers of z − z 0 , we have a series in powers of the analytic function φ(z) − φ(z 0 ). It is used in the reversion of series which consists, starting from z = f (w), in expressing w as w = g(z The general results derived are illustrated by application to a problem concerning the transformation of asymptotic relationships. Notation. Let n > 0 be fixed finite integer. (1) With w a fixed point in C, the finite part of the complex plane, Ω(w) an open set of points in C\{w} with limit point at w,Ω(w) ≡ {w} Ω(w), and p, q mappings ofΩ(w) into C, (2) Let M be a nonvoid set of points in C. For each w ∈ M, let Ω(w),Ω(w) be as in (1) and means that for each w ∈ M, p(w, z) (w,Ω) q(z). [We have]Ω = w∈MΩ (w) and The main problem considered has a simple form as follows: let a j , f j ∈ C (j : n) [This notation means that j runs from 1 to n. The notation (j : 0, n) means that j goes from 0 to n] with a 1 = 0 and, with w ∈ C fixed, let a, f :Ω(w) → C be such that In the general form of the problem a j , f j : M → C (j : n) are mappings with a 1 (w) = 0 for all w ∈ M; for each w ∈ M, a(w, ·) :Ω(w) → C is such that are to be determined.
Then, Wynn described three variants of the problem corresponding to various properties of the function a. The particular case of polynomials is treated. The text ends by: The function c(z) = ln(1 + z) is inverse to a(z) = e z − 1. c(e iθ ) = ln{2(1 + cos θ)} + iθ/2 maps the segment −π ≤ θ ≤ π onto a curve C, symmetric about the real axis, enclosing its nonpositive part, containing the real point ln(z), the imaginary points ±iπ/3, and having as asymptotes the lines z = ±iπ/2, D being the open domain bounded by C, a maps D bijectively onto the unit open disc.

On Stability Functions
The handwritten paper On stability functions (mentioned in {26} as "Paper, latest version") contains 187 pages with a bibliography at the end. Its first section is an introduction and a presentation of the notations: This paper is concerned with functions of the form for ν = 0, 1, . . . exist. A function f of this form will be called an F-function. If the context permits, the notation f ∈ F or, where convenience dictates, f (z) ∈ F, will be used. The function s in the representation (1) plays a significant role in the theory of the function f . A function of the form (2) with ψ as described will be called an S-function; again the notations s ∈ S or s(z) ∈ S will be used. The mapping properties of F-and S-functions, in particular, will be investigated.
Then, Wynn stated that an F-function f is real for finite negative values of z with f (z) > 0 for all sufficiently small negative values of z, is asymptotically represented as z tends to zero over an open set in the sector π/2, 3π/2 with the limit point 0, and that the function f can be asymptotically represented by a series of the form f (z) = ∑ ∞ ν=0 c ν z ν (formula scratched) which generates an associated continued fraction whose even convergents map the closed infinite left half-plane L into the closed unit disc D if and only if f is an F-function. The proof is given with a reference to his paper (Reference [98]) for details. An example of an F-function is e z , and Wynn added: The study of F-functions was motivated by the following consideration: With A a bounded linear operator, the solution of the differential equation for n = 0, 1, . . . If 0 < h < ∞ and the eigenvalues of A lie in L, those of exp(Ah) lie in D. y(nh) remains bounded and, indeed, decreases to zero as n increases indefinitely. An approximation y * (t) to the solution of Equation (9) may be obtained by use of a Taylor series method based upon use of an approximate identity Setting m = n = r and taking the a i and b i to be the denominator and numerator coefficients of powers of z in C 2r (hz), where C 2r (z), with r ≥ 1 fixed, is a convergent of the continued fraction associated with the exponential series, and setting t = nh, the special form of the approximate identity (11) applied to Equation (10) yields the relationship for n = 0, 1, . . . As a consequence of the mapping properties of C 2r (t) described above, the eigenvalues of C 2r (Ah) lie in D, and the remarks concerning the behaviour of y(nh) apply with equal force to y * (nh) : the exact and approximate solutions of Equation (9) behave in the same way. The practical details of the way in which relation (12) is implemented are not of immediate concern; any method for the approximate solution of Equation (9) based upon use of recursion (12) is stable.
Let us remember that the convergent C 2r of the associated continued fraction to a series is its [r/r] Padé approximant. Since the computation of Padé approximants to the exponential function are highly numerically unstable (see References [123,124]) one can doubt the practical usefulness of the procedure mentioned.
Other examples of F-functions are given. Wynn claimed that they open up the possibility of constructing stable schemes for the approximate solution of certain nonlinear differential equations. A characterisation of F-functions which is independent of the continued fraction theory is given. Based on the proof of this result, Wynn asserted that it is possible to demonstrate the existence of functions with mapping properties less specific than those of F-functions. Indeed, a number of F-functions can be derived from a given F-function. This remark reminds the way some totally monotone sequences can be derived from a totally monotone one as explained in Reference [125]. Then, Wynn proved that F-functions are closed with respect to multiplication. Meromorphic F-functions are then considered. The next part of the work deals with variation diminishing functions. Let (x i ) be a sequence of real numbers. The transformation (x i ) −→ (y j ) given by y j = ∑ n i=0 a n−i x i , j = 0, 1, . . . (it seems that n should be replaced by j), is said to be variation diminishing if the number of changes of sign of the y j is less than or equal to the number of changes of sign of the x i . Wynn wrote: Transformations of this type underly the theory of many smoothing operations, and also occur in the numerical solution of certain partial differential equations by iterative methods, and he referred to Reference [126]. Several results are proved. The next section of this manuscript is devoted to interpolation. Wynn wrote: It is possible to construct a rational F-function whose derivative values agree with those of a generating F-function up to prescribed orders not only at the origin but at a prescribed sequence of points in L; furthermore this rational function may be derived by the use of purely algebraic methods of rational function interpolation.
A theorem is proved and recurrences for the coefficients occurring in the interpolating rational F-function are given. A long discussion, where orthogonal polynomials play a role, follows. The algebraic problem of determining a rational function with denominator and numerator of degrees equal to m and m − 1 respectively, which satisfies the interpolation conditions may be solved by a recursive process which is described and justified in his Appendix 3 Extremal solutions of the Pick-Nevanlinna problem. It could be of interest to code and test this algorithm, and to compare it with the other existing ones. Exponential fitting forms the topic of the following section. Since differential systems of the form (9) constitute a very restricted class, Wynn now considered the system dy(x)/dx = f (y(x)). Linearizing it yields dy(x)/dx = f (0) + Jy(x) where J is the Jacobian matrix of f . Assuming its nonsingularity, we have Assuming that the eigenvalues λ j of J are distinct, and . . . what follows is not clearly stated and it seems that one page of the manuscript is missing. Anyway, Wynn constructed a rational function such that its derivatives of prescribed orders agree with those of e Jh at the points λ j h, and he wrote: It is eminently desirable that a method for the construction of general rational functions mapping L into D and having prescribed orders of contact with the exponential function at the origin and other specified points in L should be made available. This is precisely the service offered by Section [not identified] of this paper. In Theorem 4 [not identified], the F-function producing the interpolation data is taken to be the exponential function; the results of that theorem then show that the required function is obtained simply by means of rational function interpolation.
A section on meromorphic F-functions follows in which the properties of the corresponding continued fractions they generate are studied. The even and odd convergents are examined. The first Appendix has the title The asymptotic expansion of positive real functions. Appendix two is on The construction of functions belonging to certain classes. It deals with a general theory of the derivation of functions of the form from others of the same form. In a previous Appendix (without a number) Wynn already treated the same problem for S-functions. The third Appendix, already mentioned above, considers the Pick-Nevanlinna problem in a wider setting. Let us remember that this problem consists in finding a holomorphic function ϕ that interpolates the data λ 1 , . . . , λ n ∈ D subject to the constraint |ϕ(λ)| < 1 for all λ ∈ D; see Reference [127]. Wynn wrote: It is clear from the above conspectus of results from the Pick-Nevanlinna theory, that the solution of the problem of determining a rational function which satisfies the mapping and interpolation conditions described above differs from that of constructing a rational function which satisfies interpolation conditions alone in at least two respects: in the solution of the first problem a combination of function-theoretic and algebraic methods is involved, while in that of the second, algebraic methods are exclusively deployed; furthermore, each stage of the solution of the first problem results, not in the construction of a single interpolating function as is the case for the second problem, but in that of a family of functions with the required properties.
Developments and theoretical results follow. A bibliography of 26 items ends the paper. One can wonder why Wynn never published it.

3
The Hamburger-Pick-Nevanlinna problem This manuscript contains 179 pages but with many portions scratched (it probably contains a mixing of a draft of a paper, notes and rough notes all cited in {26,28}). Let p, q, r be finite positive integers, let a ν (ν : 0, 2p − 2) be finite real numbers, let x k (k : q) be distinct real argument values, T(1, k) (k : q) be the corresponding finite positive integers, and b k,ν (k : q | ν : 0, 2T(1, k) − 1) sets of finite real valued coefficients, let λ k (k : r) be distinct argument values in L, T(2, k) (k : r) be the corresponding finite positive integers, and c k,τ (k : r | τ : 0, T(2, k) − 1) sets of finite complex valued coefficients. The problem is to determine a function G which satisfies the asymptotic relationship as λ tends to infinity in a sector of the form π < γ ≤ arg(λ) ≤ δ < 2π, and satisfies also the asymptotic relationship for k : q as λ tends to x k over an open set contained in L with limit point at x k , for k : r, and moreover satisfies the interpolation conditions λ τ G(λ k )/τ! = c k,τ for τ : 0, T(2, k) − 1, and finally maps L into D (or U?). According to Wynn The proof of the principal theorem of this section is largely based upon results, due to M. Riesz, in the theory of linear functionals.
Since this part is a draft in quite bad shape, we will not pursue its analysis.

The document contains another interesting section
Matrix criteria In this section Akhiezer's treatment of the simple Pick-Nevanlinna problem is extended to the diminished Hamburger-Pick-Nevanlinna problem; conditions that are necessary and sufficient for the solubility of the latter problem are established and, assuming this problem to be nondegenerate, inclusion discs for its solutions are located.
As a preliminary, a Hermitian matrix is constructed from the data which leads to an extension of a theorem on the existence of inclusion discs for the values of all solutions to the Hamburger problem and the simple Pick-Nevanlinna problem. Wynn concluded: The results of the above theorem are exclusively concerned with the diminished version of the Hamburger-Pick-Nevanlinna problem. It is possible to extend the method of proof to the treatment of the Hamburger-Pick-Nevanlinna problem itself, to examine the structure of the rational function solutions to this problem and its diminished form, and to describe the relative positions of the inclusion discs deriving from a sequence of subordinate problems. These matters are, however, more conveniently dealt with methods described in the following section, in which explicit expressions for the general solutions to the interpolation problems concerned are described.
Unfortunately, this following section of the manuscript is not under a form which allows to give a clear account of it.

Continued fraction transformations of the Euler-Mclaurin series
The document is 202 pages long. It dates from December 1976. The first part of the document is in a quite good shape (probably a draft of the paper mentioned in {26}, with an abstract and a bibliography) with some corrections done with a pencil. Not all pages have been written with the same pen. Some of them are missing or not in the right order since all of them are not numbered (they could be some additional notes on integral transform and analytic continuation). Moreover, references to some formulas are missing, and some authors quoted in the text are not listed into the bibliography given at the end. An in-depth study of this document is needed to fully understand it if possible to exploit it. It begins by a long abstract from which we extract the main points.
Results concerning the convergence of forward diagonal sequences of quotients in the Padé table are given. In particular, it is shown that, if ( * ) f ν = ∞ 0 t ν dσ(t) (ν = 0, 1, . . .), σ being a bounded nondecreasing real valued function such that all moments ( * ) exist, and ( * * ) f ν = O{(χν)!ξ ν } (0 < χ ≤ 2, 0 < ξ < ∞) then all forward diagonal sequences of Padé quotients derived from the series ∑ ∞ ν=0 f ν z ν converge uniformly over any bounded region in the z-plane not containing any point of the nonnegative real axis to ( * * * ) f (z) = ∞ 0 (1 − zt) −1 dσ(t), f (z) being the le Roy or (B, χ) sum of the given series for all finite z in the sector χπ This result extends results given by various authors, and it can also be extended to the case where the lower limit in ( * ) and ( * * * ) is replaced by −∞. These results are applied to the delayed Euler-Mclaurin series (j ≥ 0, and D being d/dµ) regarded as a series expansion in ascending powers of k 2 . Convergence results for the Padé approximants of this series are derived, and also for the same series in which single zeros are inserted between its successive terms. Applications to Stirling's asymptotic expansion of the logarithm of the gamma function, and to the asymptotic expansion of the generalised Riemann zeta-function are presented.

Convergence and truncation error bounds for associated continued fractions
This is a short document (13 pages), probably only the notes mentioned in {26}.
In this work, he proposes to give a convergence theorem for the functions 2r {ψ(µ)}, where, presumably, 2r is the 2rth convergent of Thiele's continued fraction. He begins to prove other results (not reproduced herein since they contain too many erasures) upon which the proof of the theorem is based. It is: be a fixed interval of the finite real axis, and ξ(s) be a bounded nondecreasing real valued function for α ≤ s ≤ β, and not a step function with a finite number of salti. Let µ ∈ (−∞, ∞) × [α, β] be fixed and Then According to what precedes, it seems that the sign that looks like × in µ ∈ (−∞, ∞) × [α, β] has to be replaced by / ∈ [α, β]. We do not know if this result had been later rediscovered by other researchers. Thus, it was interesting to reproduce it here.

Functional Interpolation
This very well handwritten manuscript has 52 pages and seems to be the paper mentioned several times in all the list of documents, projects and activities, with a bibliography with a last reference dating from 1984. The first section is Interpolation and extrapolation with a subsection 1.1-Procrustean technique. This word was already used by Wynn in the title of the published paper of Reference [2] where he introduced a particular form of the -algorithm for extrapolation at infinity by a rational function in n (see Reference [113]). It describes situations where different properties are fitted to an arbitrary standard. In the Greek myth, Procrustes was a son of Poseidon. He compelled travelers to lie on a bed, he cut off their legs that were longer than bed, and stretched the feet of those who were too small. The manuscript begins by: Most general theories arise from investigations of particular problems, and in this respect the theory to be described is not exceptional. By way of motivation, the problem of deriving an extrapolation method from an interpolatory formula and its converse are considered. It is first supposed that an interpolatory function of complex variables is available for which (a) F (m) r (d | y; y ν ) = d ν (ν ∈ I m,m+r ; throughout the paper I i,j is the sequence i, i + 1, . . . , i + j; I i is the sequence i, i + 1, . . .; I is I 0 ) and (b) for certain distributions of the d ν and y ν G (m) exists. The function F (m) r serves as the basis of an extrapolation method: given the sequence   2i (d | y) gives the extrapolated limit. These numbers can be recursively computed by Thiele's reciprocal difference algorithm (see (3) in Section 3.4), and they can be displayed in a two dimensional array for which Wynn does not use the usual notation but a new one, and he writes: [The process] involves numbers ε i,j which may be set at the intersections of the full rows and columns and of the half rows and columns of a chipped triangular array in which the row index i ranges overĪ −1/2 (Ī k is the sequence k, k + 1/2, k + 1, . . . ;Ī isĪ 0 ) and the column index j over the range I i , the number ε −1/2,−1/2 being missing. The numbers ε i,j are constructed from the initial vales ε 0,j = S j (j ∈ I), ε 0,j = 0 (j ∈ I) (the dash is used to indicate a displacement operation acting upon numbers with two suffixes, whose effect is illustrated by the relationships ε 0,j = ε −1/2,j+1/2 and ε i,j = ε i−1/2,j+1/2 ) by use of the relationship And he concludes with his personal sense of humor: We have called an extrapolation method of the above type a Procrustean technique [2], although in fitting function values to a sequence, i.e. the bed to the victim, we are a little kinder than Procrustes is reputed to have been.

He continues:
If in formula (1) λ is fixed and one of the y ν is very large, e ν [not defined] is approximately equal to lim F  There should be an error in what Wynn wrote. Since he states that λ is fixed, λ → ∞ should be replaced by y ν → ∞ (remark of one of the reviewers). Wynn concludes that this property is, for example, satisfied by the ratio of determinants expressing G (m) r (S | y) which, thus, tends to S ν as y ν tends to infinity.
The second subsection is named 1.2-Interpolatory functions. He claims that The steps taken in the above derivation of an extrapolation method from an interpolatory formula may be reversed. Then, he shows how to reverse the recursive rule given in the preceding subsection.
[. . . ] set d ν = d(y ν ), e ν (z) = e(y ν , z) and f ν (z) = f (y ν , z) (ν ∈ I). Suppose that a function of complex variables  i,j (z) = d ν as z → y ν (ν ∈ I m,m+i+j ). Imposing the further condition it follows that The interpolatory function derived from the extrapolation estimate i,i (z). He proves that lim r (m) i,j (z) = d ν as z → y ν (ν ∈ I m,m+i+j ). Then Setting now ε 0,j = 0, ε 0,j = e j (z) (j ∈ I) and w i, In this result, the conditions of the simple case in which (8) e(y, z) = d(y), f (y, z) = z − y may be imposed upon e and f . Now e ν (z) = d ν is a constant, independent of z (ν ∈ I) and f ν (z) is the difference z − y ν (ν ∈ I). In this case, r i,i is the quotient of two ith degree polynomials in z, the rational function from which the extrapolation limit (m) 2i (S | y) was derived. The above process now reduces to an algorithm for rational function interpolation due to Brezinski [128] of which a generalisation has been proposed by Cordellier [129]. The discussion of interpolatory formulae and extrapolation methods is terminated by the remark that under appropriate conditions cyclic derivation of extrapolation methods from interpolatory formulae and conversely may be repeated indefinitely.
Section 2 of this document is entitled Approximants of general order. In the simple case i,j is a ratio of two polynomials of degree i. In this Section, he proposes to study the more general system of approximants r Recursive relations for the L i (z)'s are also given. Using divided differences, they are also expressed in Newton form. Determinantal formulae are related to the recursive rules given. It is showed that particular cases for the e ν and the f ν give back the usual Lagrange interpolation formula, the Neville-Aitken scheme, and Newton series. The work of Stoer on interpolation by rational functions [131] and the variant of the ε-algorithm due to Claessens [132] are also recovered as special cases. The next subsection is named Termination. When e(y, z) is a polynomial or rational function of f (y, z), termination of the algorithms previously given is proved. The following subsection treats Confluence that is when some argument values coincide. Wynn examines what happens to the previous formulae and recursions. In that case, for n ∈ I 2 , and formulae for the computation of the coefficients c ν and d ν are given. They implement a truncated composition of polynomials. The d ν 's are the coefficients of the Newton series representations of the corresponding Lagrange forms, which are confluent forms of the functional divided differences. They are related to the b ν , but all details are too complex to be given herein.
Then, comes a subsection on Zero finding algorithms. Under suitable conditions, the above algorithm for the truncated composition of polynomials can be used for the inversion of formal power series. Moreover In so doing it serves as the basis of a number of algorithms for determining the zeros of a function and motivates the use of the approximants r (m) i,j (z) for the same purpose. Setting e(y, z) = z − y and f (y, z) = φ(z) − φ(y), φ being the function under treatment, relationships of the form (31,32) hold with x(y, z) = z − y, c ν (y, z) = φ (ν) (y)/ν! (ν ∈ I 1 ) are Taylor series coefficients, d(y) = 0, b 1 (y, z) = 1 and b ν (y, z) = 0 (ν ∈ I 2 ).
For example, if y 1 = y 2 , L (1) 1 (z) leads to the method regula falsi z = (y 2 φ 1 − y 1 φ 2 )/(φ 1 − φ 2 ). If y 0 = y 1 , then Wynn obtains two combinations of Newton's method and regula falsi In the same way, an nth order single point iteration process can be obtained. The subsection ends by The more general theory yields multipoint processes (for a further application of the εalgorithm to the problem of finding the zeros of a function, see [133,134]).
The last subsection of this document is entitled Extensions of the Lagrange-Bürmann expansion. Wynn claims that The above treatment of the confluent case offers an interpretation of the theory of this paper. He first gives the Lagrange-Bürmann expansion of d(z) in powers of φ(z) − φ(y µ ) where y µ belongs to a close contour in the complex plane. Then, he obtains an algorithm for determining the coefficients in the Lagrange-Bürmann expansion and in an asymptotic version of it from the Taylor series coefficients of d(z) − d(y) and φ(z) − φ(y) at the point y = y µ . The case of confluence is also treated. In addition, in this case, we do not know why Wynn never published this work since it was ready to be submitted. 7 Interpolation by the use of rational functions This handwritten complete paper of 90 pages with a bibliography of 27 items, the last one dated 1979, is present in the projects and in the Lists of documents, and it was never published. It seems to be related to the previous manuscript {6}.
In the Section 1 of this document, titled The Thiele-Nörlund interpolation theory, Wynn reminds how to construct rational interpolating functions C (m) 2s in which the numerator and the denominator have degree s, and C (m) 2s+1 with a numerator of degree s + 1 and a denominator of degree s and such that C The successive numerators and denominators are recursively computed by Wynn asserts that conditions to ensure that all rational functions C (m) r can be constructed by the above scheme and that they have the required interpolation properties can be formulated in terms of determinants, and that determinantal formulae can also be given for the numbers and the functions involved, and he claims that Such formulae are made more concise by the use of a special notation. It takes 6 pages to establish these notations. After quoting a remark in German by Nörlund [130] (Ch. 15, p. 420) that Wynn finds perhaps a little exuberant, he writes that: It is the principal purpose of this paper to point out that, using another very simple relationship (namely, in particular, that, if x i is replaced by (λ − x i ) −1 , (m) 2r becomes C (m) 2r (λ)) many results suggested by the behaviour of reciprocal differences may be obtained for convergents.
These procedures are described in Section 2 of this document, titled The σ-and µalgorithms. After explaining how to obtain them, Wynn writes  2s (λ) has already produced the σ-algorithm of relationship (16), the simplest and most economical method, subject to the stated conditions, for evaluating these functions; the observation also leads directly, as will be shown below, to new interpolatory theory. Once made, the observation is trivial, and its implications are not difficult to work out; perhaps its most interesting feature is that it has not been made before. The simple relationship between reciprocal differences and rational interpolating functions was not so much discovered as forced upon the author's attention while working out the consequences of principles underlying the process of interpolation and the transformation of divergent series. These two subjects have recently became increasingly important in computational mathematics; new theory in what once might have been considered dead subjects is constantly being developed (mention may be made of recent generalisations of polynomial interpolation described in [135][136][137][138]); it is highly probable that the principles concerned will find further applications, and for this reason they are now outlined.
It is difficult to summarize what follows without quoting large parts of the document. Moreover, some notations and their inferences are not clear enough. Basically, Wynn comes back to the link between interpolation and extrapolation already discussed in the document {6}, that he named Functional Interpolation. In particular, he considered the following extended ε-algorithm (which contains the -algorithm and some other extensions [128]) He considers the particular case γ (m) r = 1, which corresponds to the ε-algorithm. He reminds that, when applied to the partial sums of a formal power series, this algorithm furnishes the Padé approximants belonging to the lower half, diagonal included, of the Padé table, and that he derived (no reference) various determinantal formulae from a more general form of approximating fractions given by Jacobi [139]. He also mentions that, in its special form, the problem was also studied by Frobenius [140] and Padé [141], and he adds the following remark in which he explains how he obtained his ε-algorithm While idly investigating the formulae resulting from the choice γ (m) r = 1 in the relationship (25) [. . . ] the author noted that expressions obtained for the numbers ε (m) 2r were equivalent to extrapolatory determinantal expressions, simplified versions of those due to Jacobi and used by Frobenius, previously published by Schmidt [142] and republished by Shanks [112]. In this way the ε-algorithm was discovered.
Then, Wynn notices that (25) has been used by Claessens with γ (m) r = (λ − x m+r+1 ) −1 , and applied to the partial sums of the Newton interpolation series, for obtaining rational interpolating functions [132], and he explains the theoretical basis of this algorithm Claessens was led to the discovery of the extended ε-algorithm by the consideration of interpolatory continued fractions not of Thiele form, but of a form introduced by Kronecker [143] in connection with a process initiated by Rosenhain [144] and Borchardt [145] for constructing the resolvent of two polynomials from systems of their numerical values.
If all the x i tend to a common value x, the interpolation fractions tend to the Padé approximants and the extended ε-algorithm (25) tends to the usual one. When γ (m) r = m + r + 1, (25) gives back the -algorithm studied by Wynn in Reference [2].
The following Section is on Lozenge algorithms, which are algorithms relating quantities located at the four corners of a lozenge in a table where the lower index indicates a columns and the upper index a descending diagonal. The ε-algorithm and its generalizations [146], the -algorithm, and the algorithms numbered by Wynn (23) (see the manuscript {6}), (16) and (25) (see above in this manuscript) enter into this class. They share some algebraic properties when, instead of applying them to a sequence (S n ), they are applied to (aS n + b), when γ (m) r in (25) is multiplied by a constant, and when a fractional linear transformation is applied to the elements with an even lower index (property named homographic invariance). In these algorithms, the quantities with an odd lower index are only intermediate results with no interest for their purpose. They can be eliminated and a new rule relating five quantities disposed at the center and at the extremities of a cross are obtained. The first algorithm to have been treated in that way is the ε-algorithm for which Wynn obtained this cross rule [48]. Concerning this rule, he writes in the document under consideration Despite the author's pianissimo protests, this relationship and others of its kind have been referred to by various writers as Wynn identities. This extended cross rule and the homographic invariance property were independently presented by Cordellier in his Doctoral Thesis [147]. They were also published by the same author in Reference [129] without, of course, knowing this unpublished document by Wynn. On the contrary, Wynn was knowing Cordellier's results through his correspondence with C.B., and he indeed refers to Reference [129] for the extension of his cross rule. The homographic invariance for the σ's and the µ's presented above is treated. In the same document, Wynn also discusses the possible breakdown of these relations when a denominator becomes zero. Singular rules for continuing the computations in the case of only one isolated singularity are given. They extend the singular rules he gave in References [20,36] for the usual ε-algorithm. All these results are gathered into Theorem 2 in which Wynn also gives the cross rule for formulae (16) and (25). If the relationship γ 2s+1 is satisfied between the γ, then the cross rule also holds as described with the five ε's replaced by (αε + β)(γε + δ) −1 . Then, if an extended cross rule holds between five numbers as explained above, Wynn looks at which conditions a recurrence relationship of the form (25) holds between them. The following Section deals with Invariance properties. Denoting by (m) r (x, f ) the reciprocal differences obtained by Thiele's formula (that is the -algorithm given by (3)) from the argument and functions values x i , f i , i = 0, 1, . . ., and setting y i = Ax i + B and g i = (α f i + β)/(γ f i + δ) with αδ − βγ = 0, Wynn proves that If γ = 0, a property for (m) 2s (y, g) is also given. These invariance properties have their counterparts for the convergents C (m) s . Section 4 is on Confluent algorithms. Up to now, all the x i have been assumed to be distinct. Now Wynn sets x i = x + ih and allows h to tend to zero. The functions C (m) r (λ) evolve, in particular, to a confluent form C r (λ, x) and become the successive convergents of the continued fraction corresponding to the Taylor series The coefficients of this continued fraction may be expressed in terms of Thiele's reciprocal derivatives. He refers to Reference [130] (Ch. 15) for the corresponding theory, and explains that the confluent forms of the algorithms of formulae (16) are with σ −1 (λ, x) = 0 and σ 0 (λ, x) = f (x). Determinantal formulae for these σ's are given. The functions µ are treated similarly. Complete proofs are given.

Sequel to the abstract theory of the ε-algorithm
In February 1971, when he was working at the Centre de Recherches Mathématiques of the Université de Montréal, Wynn published a report on the abstract theory of the ε-algorithm [69]. In the boxes he left with his friends, we found a document with the title Notes on sequel to the abstract theory of the ε-algorithm. All topics treated are in the noncommutative case. It probably corresponds to the rough notes indicated in {26}. It consists of 153 pages including numerous corrections and additions as well as strikethrough passages, and pages from Reference [69] intercalated. In the lists, a corresponding paper is also mentioned, but, unfortunately, we have not found it yet. The document begins with a discussion about the application of Euclid's algorithm to rational functions in a ring with a zero and a unit element. Another one on noncommutative continued fractions follows. Recurrence relations for the numerators and the denominators of their successive convergents are given. They are the same as the usual ones. Wynn claims that: Then, he discusses Euclid's algorithm for formal power series with coefficients over a ring. Wynn proves that they have two sided reciprocal series. It is difficult to extract some results from this part without much work. Then, he considers systems of continued fractions derived from a single power series. Notions on C-regularity and semi-normality are introduced. Applying a recursive algorithm whose rule is that of r ) −1 , but without giving the definition of the inverse, to the partial sums of a semi-normal series he proves that E (m) 2r = D(z) −1 N(z) for a certain value of r. Other theoretical results are given, but they are difficult to exploit. Then, pages 93-96 of Reference [69] are reproduced which leads to think that what precedes is a complement to that report. Wynn then shows that both halves of the Padé table can be constructed via this algorithm under certain assumptions. A cross-rule is given for each half. Pages on vector continued fractions follow. Arithmetic operations upon vectors and formal power series with vector coefficients are first explained. But the pages, the concepts (vector valued rational functions, McLeod isomorphism [148], . . . ), and the results are so entangled that it seems difficult to extract something coherent. There are pages on Euclid's algorithm for formal power series with vector valued coefficients, a vector qd-algorithm, and the vector ε-algorithm. However, an interesting result is given but it needs to be verified. Let us give it after a simplification of the notations. Summations are over ν: Theorem. If the finite dimensional vectors of the sequence { f ν } satisfy a system of relationships of the form where the {d ν } are real numbers, then the series ∑ f ν z ν is generated by a vector values rational function f (z). If vector value rational functions ε for r = 0, . . . , 2r − 1, m = 0, 1, . . ., then ε Looking at the proof, which is not complete, one can understand that the words the series ∑ f ν z ν is generated by the rational function f (z) means that where ∑ r 0 d ν f r−ν = η r (r = 0, . . . , r − 1). The vector valued Padé table is presented. Then, Wynn gives some elementary results in the theory of ε-algorithm as applied to sequences of elements of a ring. In particular, a cross rule is obtained. Rational functions with coefficients over a ring are discussed. Then, Wynn studies the application of the ε-algorithm to a sequence of numbers satisfying an inhomogeneous linear difference equation. He proves the following result (notations slightly changed) Theorem. Let the numbers S ν satisfy the recursion ∑ r 0 d ν S m+ν = H, m = 0, 1, . . . and no recursion of a similar form with r replaced by a smaller integer; set ∑ r 0 d ν = D. It is not possible that both D and H are both factors of zero.
It is then proved that, for such a sequence, ε (m) 2r = D −1 H for all m. Other results of a lower interest are given. Then, Wynn looks at the properties of the hierarchies of ε-arrays, as defined in Reference [79], that are obtained for such sequences. Since pages 23-28 of his report from the Université de Montréal corresponding to the paper of Reference [79] are inserted here between the pages of the manuscript, and preceded by If followed by a difference equation, what precedes can be considered as an addition to them. Using the vector-matrix isomorphism of McLeod [148], Wynn writes that, for recursions involving vectors, it is a relatively simple matter to derive results analogous to those of the preceding section. Further notes end the document. 9 How to find the centre of a spiral This handwritten document contains 49 pages that correspond to the slides for a talk. There are parts in the text separated by numbers in bold from 1 to 30 in large square boxes (see {28}). Its title could be How to find the centre of a spiral, and it was delivered by Wynn after 1987 since a book dating 1987 is quoted. After there are 30 pages where the same content is detailed and illustrated very carefully, probably some notes to copy and distribute. The document gives an interesting new interpretation of extrapolation by the ε-algorithm.
We start with a fixed complex number λ. If the modulus of λ is less than 1, λ lies inside the unit circle. If the modulus of λ is greater than 1, λ lies outside the unit circle. Next, we form the successive powers λ 0 , λ 1 , λ 2 , . . . of λ. They lie on a spiral in the complex plane. The argument of λ 2 is twice arg(λ), the argument of λ 3 is three times arg(λ) and so on. If the modulus of λ is less than 1, the powers of |λ| tend to zero and the spiral contracts to the origin. If the modulus of λ is greater than 1, the spiral expands away from the origin. Then, we multiply each of the points λ 0 , λ 1 , λ 2 , . . . by a constant complex number b. This does no more than rotate the spiral through arg(b) and magnify or diminish it by a factor of |b|. But a contracting spiral remains contracting and an expanding spiral remains expanding. Lastly we add a constant C to each of the transformed points and produce members of a first order spiral sequence. In the case in which |λ| < 1, the sequence converges to C and the spiral upon which its members lie contracts toward C. In the case in which |λ| > 1, the spiral expands away from C. If lim i→∞ s i is finite, lim i→∞ s i = C. Already we may pose a problem in connection with first order spiral sequences. It is known that three successive members S k , S k+1 and S k+2 of a sequence have the form but the values of C, b, λ and k are unknown. Find the value of C alone. If we can solve this problem then what we are able to do is fit a first order spiral sequence to the subsequence S k , S k+1 and S k+2 and find its centre.
Then, Wynn shows that the problem can be solved by computing ε (k) 2 by a rule which is that of the ε-algorithm (without mentioning its name). He gave the formulae for ε (k) 1 and ε (k+1) 1 . Then, he computes ε (k) 2 , and shows that it holds ε (k) 2 = C. Thus, the spiral has been fitted without the determination of the values of b and λ. Then, he considers three other members S k+1 , S k+2 and S k+3 having the form Another spiral can be fitted through these three numbers, and a similar computation is identical to the one computed before. The same process can be continued with S k+2 , S k+3 and S k+4 , and so on. For each value of the upper index, the center of the spiral can be computed by the rule . . This is the ε-algorithm [3]. Wynn mentions that, if S 0 , S 1 , S 2 , . . . are successive members of the same first order spiral sequence, that is S i = C + bλ i for i = 0, 1, . . ., then ε 2 ) which is, in fact, Aitken's ∆ 2 process. Then, Wynn considers a second spiral entangled into the first one.
The members of a second order spiral sequence have the representation . This sequence has two forms: either (a) only one geometric term accompanied by a linear function of i is present and or (b) two geometric terms accompanied by constant coefficients are present and In the latter case a second order spiral sequence may be regarded as a first order spiral sequence whose members are represented by b 2,0 λ i 2 whose center moves on a further first order spiral upon which the points C + b 1,0 λ i 1 lie.
In these forms, the kernel of the transformation (S m ) −→ (ε (m) 4 ) is recognized. Then, Wynn considers nth order spirals with the s i satisfying the same formula but now with the condition ∑ h j=1 {τ(j) + 1} = n, and he discusses their various special forms. The center C of such spirals can be obtained via the ε-algorithm which, starting from S k , . . . , S k+2n , delivers ε (k) 2n = C. Since he writes that: There is no time to give a proof of this result this confirms that these notes were written for lectures. The same process is repeated with S k+1 , . . . , S k+2n+1 , and it gives ε (k+1) 2n = C . As before, the two schemes have some ε in common, and, if S 0 , S 1 , . . . are successive members of the same nth order spiral, then, for all m, ε (m) 2n = C. This is a known property of the ε-algorithm since the sequence (S n ) belongs to the kernel of the transformation (S m ) −→ (ε (m) 2n ).
If the parameters C, b and λ in the equation of the first order spiral are real, the s i = C + bλ i , i = 0, 1, . . ., still lie on a spiral in the complex plane, but the numbers s i themselves lie on the real axis. In that case, the ε are also real. Then, Wynn explains how the process for finding the centres of spirals which we have described serve as the basis of a method for obtaining estimates of the limit of a sequence.
Starting from S 0 , S 1 and S 2 , the center ε (0) 2 of the corresponding first order spiral is an estimate of the limit. Then, from S 1 , S 2 and S 3 , a second estimate ε (1) 2 is obtained, and so on. Now, from S 0 , . . . , S 4 , the center ε (0) 4 of the second order spiral passing through these numbers furnishes another estimate of the limit. The process can be repeated by increasing the number of terms of the sequence (S m ) used in the construction of the successive spirals. And Wynn writes This method of estimating the limit of a sequence is known as the ε-algorithm [...]A question now arises: does this method work? Wynn gives two numerical examples: partial sums of the series 1 − 1/2 + 1/3 − 1/4 + · · · which converge to ln 2, and those of the series ∑ i C i −1/2 1/(i + 1) (C i −1/2 represents the binomial coefficient) which tends to 2(2 1/2 − 1). He tentatively explains the success of the process: Perhaps a happier interpretation is to suggest that a great deal of information lies lurking in the first few members of a sequence: we have only to think of a way of getting it out. This interpretation serves us a little better when we come to consider the transformation of divergent sequences.
In saying that, Wynn is too much optimistic. Indeed, the behavior of a sequence can completely change after a certain number of terms, and change again later. As proved in Reference [149], an algorithm to accelerate the convergence of all sequences cannot exist. It is the same even for restricted classes of sequences [150]. Then, Wynn gives examples of divergent series and explains that the reasons why the ε-algorithm works in these examples are found in the classical theory of continued fractions. After that, he describes the use of the algorithm in the solution of the fixed point problem x = f (x). He explains graphically the convergence of the Picard iterates, and what the notion of order of convergence of a sequence is. Since Newton's and higher order methods require the use of the derivatives of f , he explains how this drawback can be avoided. He starts from an initial estimate C (0) of the fixed point, sets S 0 = C (0) , performs the iterations S i+1 = f (S i ) for i = 0, . . . , 2r − 1, and finds the center C (1) = ε (0) 2r of the spiral constructed from these iterates, and restarts the Picard iterates from S 0 = C (1) . When r = 1, this is exactly Steffensen's method which has order 2. Wynn claims that under suitable conditions C (k+1) − x = O((C (k) − x) r ). Let us mention that this claim has never been proved. Then, Wynn explains that the ε-algorithm can be applied to sequences of functions, of vectors, of a field or of any mathematical system over which addition, subtraction and the formation of an inverse are defined, and that the theory of continued fractions has been developed for such cases. The continuous case of the preceding method is considered. It is aimed at the estimation of lim t→∞ S(t). In what precedes, i is replaced by t and λ i j by e −α j t , and he names an rth order spiral function, the expression Wynn proposes to estimate lim t→∞ S(t) by fitting a spiral function s(t) to S(t) at one value x of t, and to take its center as this estimate. For that purpose, he assumes that lim t→x S (i) (t) are given for i = 0, . . . , 2n, but that the value of x itself is unknown. For finding the center C of the nth order spiral function, Wynn gives the following lozenge algorithm (defined in Reference [88,89]): Since Texas is an oil-producing state the last problem we consider is that of the owner of a well producing oil with cast price c per barrel who wishes to determine (a) the price P par barrel which maximizes his profit and (b) his profit at the optimal price P. We now describe a solution and hasten to say that the economic part of the argument is taken directly from the mid nineteenth century French economist Cournot who considered the problem of the owner of a spring producing mineral water at zero cost per bottle.
After giving the solution of this problem, Wynn proposes some references for those who are interesting in studying the matter further. Among them, the book of Antoine Augustin Cournot (1801-1877) [151]. As said before, at the end, there are 30 particularly nicely written pages with figures about spirals (see Figure 3), implementation schemes for the ε and the ω algorithms, an illustration of Picard iterates, numerical results, illustrations for the order of convergence, details about the problem of oil pricing policy, and references. It seems that these pages were ready to be copied for distribution to participants. Since the application was about an oil problem in Texas, this document could have been prepared for a seminar at the University of Texas at San Antonio. Let us mention that the interpretation of the scalar and vector Aitken's process and ε-algorithm by means of spirals was independently rediscovered in 2014 by Berlinet [152], together with other interesting geometrical analyses of them. We did some numerical experiments with the spirals introduced by Wynn, and, in particular, we wanted to test his drawing reproduced in Figure 3. Remember that, according to Wynn's nomenclature, a first order spiral has the form (our notation) S n = C + a 1 λ n 1 , and that a second order spiral can have two forms (a) S n = C + λ n 1 (a 1 + b 1 n) or (b) S n = C + a 1 λ n 1 + a 2 λ n 2 . If the ε-algorithm is applied to a first order spiral, then, ∀n, ε (n) 2 = C, which is Aitken's ∆ 2 process. Applied to a second order spiral of form (a) or (b), this algorithm gives ∀n, ε (n) 4 = C, while the iterated ∆ 2 process (that is the ∆ 2 process applied to the sequence (ε (n) 2 )) does not possess this property. Wynn's drawing of Figure 3 corresponds to the case (b). He wanted to illustrate the fact that S n = C + a 1 λ n 1 + a 2 λ n 2 can be written as S n = C n + a 2 λ n 2 , with C n = C + a 1 λ n 1 , which he described as the spiral a 2 λ n 2 with center moving on spiral C + a 1 λ n 1 . In our numerical experiments, we always took C = 1.5 + i. In Figure 4, we show two first order spirals with a 1 = 1 − i: the one on the left is obtained with λ 1 = 0.8 − 0.5i, and it converges to C, the one on the right is with λ 1 = −0.8 + 0.6i and does not converge or diverge since |λ 1 | = 1. Now, we consider second order spirals of forms (a) and (b). We took, for the form (a), a 1 = 2 − 4i, b 1 = 2.3 − 0.3i, λ 1 = 0.92 + 0.2i, and, for the form (b), a 1 = 1 − i, λ 1 = 0.8 + 0.5i, a 2 = 1.2 + 2i, λ 2 = 0.7 + 0.65i. All spirals are converging to C since |λ 1 | < 1 and |λ 2 | < 1. Figure 5 shows the corresponding second order spirals: (a) on the left and (b) on the right. In Figure 6, we see the results obtained by the ∆ 2 process applied to second order spirals: (a) on the left and (b) on the right. Comparing these curves with those of Figure 5 shows that, in fact, Aitken's process acts as if it was suppressing one of the spirals of the curve (which contains two of them) since the sequence (ε (n) 2 ) still looks like another spiral, and is not reduced to a single point as in the case where only one spiral was present. Indeed, if (S n ) has the form (b) with |λ 2 | > |λ 1 |, we have S n − C = O(λ n 2 ). It is easy to see that which shows the gain brought by Aitken's process by almost suppressing the role played by λ 2 , the most important one. If (S n ) has the form (a), that is S n = C + (a 1 + b 1 n)λ n , then S n − C = O(nλ n ). Aitken's process leads to ε (n) which shows the acceleration brought by suppressing the leading error term. In Figure 7, we see the results obtained by the iterated ∆ 2 process applied twice to second order spirals: (a) on the left and (b) on the right. Thus, the iterated ∆ 2 process does not gives the exact value of C, which should be a single point and not a curve although (ε (n) 2 ) is a spiral (see Figure 6), while (ε (n) 4 ), in exact arithmetic, gives exactly C. However, the iterated ∆ 2 process has a faster convergence, and it also has the form of a spiral. Spirals with different values of the parameters can have various shapes. It is sufficient, for instance, to change one sign in λ 1 or λ 2 . We intend to deepen this subject in a forthcoming paper.

Project for a Book
We separated this document from the others since it treats a large range of topics and seems to be an advanced skeleton for a book.

Book on ε-algorithm
This 120 pages document is clearly a preparatory work for a book. Looking at its contents, it is probably the Resumé of book on the ε-algorithm, mentioned by Wynn in {26}, and also at the end of his document {28}, with an O (ongoing) in front of the name. In this manuscript, he stated in detail what he intended to treat. It contains an incredible succession of theorems whose proofs are not given, and comments and explanations between them are missing. It is divided into chapters.
Chapter 1: Theorem concerning determinants. It contains identities for compound determinants, Schweins' lemmata, simple and extended Hankel determinants. Chapter 2: Prediction based upon a linear model. Its sections are on exponential extrapolation (that is, the Aitken ∆ 2 process), and extrapolation using a linear recursion (that is the Shanks transformation). There is a section on Hadamard's theory of the Taylor series, and two sections on the convergence of the rows and columns of the Padé table. Another section deals with the construction and convergence of the even order epsilon-array. Chapter 17: The first confluent form of the epsilon algorithm: the rational function limit and the continued fraction limit. This chapter has two sections where these topics are treated. Chapter 18: The second confluent form of the epsilon algorithm: the definition of an integral as the limit of a continued fraction. After defining this second confluent form, Wynn explains the connection between both forms. Then, he treats the rational function and the continued fractions integrals. Euler integration ends the chapter. Chapter 19: The rational function integral. The degenerate theory in treated in the first section, and the general theory in the second one. A special convergence result is then given. Chapter 20: The continued fraction integral. The chapter begins by the degenerate theory before coming to a remarkable equivalence. Chapter 21: The third and fourth confluent forms of the epsilon algorithm. The fourth one is given but not the third one. Chapter 22: A partial differential equation associated with the epsilon algorithm. The first section is about a so-called φ-array. Then, partial differential equations are derived for it. The cases of the epsilon algorithm and that of the Padé surface are then treated. The following section is about self conjugate systems of partial differential equations for symmetric algorithms. Adjoint partial differential equations are discussed, and special solutions of the partial differential equation of the Padé surface are given. Chapter 23: Error analyses of the epsilon algorithm. The chapter begins by a perturbation analysis. Then, the convergence and stability of the epsilon algorithm are studied. The singular rule for the algorithm is given. This preparatory work for a book is in an advanced form since the results are well positioned in it. It only remains to link them by explanations and to give the proofs of the theoretical results (or the corresponding references to the literature).

Theory of stratified commutative field
This document presents a total of 266 pages. These notes, mentioned in many of the lists of Wynn, seem in a final form because they contain no erasures. After defining his notations for the ranges of the indexes, Wynn writes (all indexes are in Z): Definition. A stratified commutative ring, or S-ring, is a system S of numbers arranged into strata S{j} such that (i) the members of S{j} form an additive Abelian group; (ii) commutative and associative multiplication is defined between elements of the various S{j}: with a ∈ S{j}, b ∈ S{k}, ab = ba ∈ S{j + k} and (ab)c = a(bc) for such products; (iii) multiplication is distributive with respect to addition within appropriate strata: with a, b ∈ S{j}, c ∈ S{k}, (a + b)c = ac + bc, addition on the left being within S{j}, that on the right within S{j + k}. A number belonging to the above system is an S-number and the S{j} to which it belongs is its stratum. S is the class of S-rings. It is not assumed in the above definition that addition and subtraction between disparate strata should be possible, although in special cases this may be so. It is so, for example, when all S{j} are the same additive Abelian group; S is then simply a commutative ring. More generally, it is also possible between congruent strata if S{j} ≡ S{k} when j = k mod m, m being a prescribed integer.
Two special S-rings are described, and then: The above examples serve to draw attention to a salient feature of theory to follow. It will be shown that the end results of certain computations involving S-numbers are members of fixed strata-for example that certain bilinear forms with S-number coefficients are expressible as linear sums of squares with coefficients that are S-numbers.
The rest of the document cannot be summarized without going into detail. It is a purely theoretical work. No practical application is given, and there are no references to the literature.

Stratified field, determinantal identities and LU decomposition in such field
This document of 105 pages probably follows the preceding document {11}. It contains theoretical notes about determinantal identities and LU decomposition over stratified commutative ring. He defines the so-called P-numbers that, in fact, have the same definition of the S-number of the document {11}. He defines the P-array and the Pmatrices that are different in their definition. Several result were proved, concerning, for instance Jacobi's theorem on the adjugate, product of determinants, Sylvester's and Schweins' determinantal identities, and so on. A LU decomposition is also proposed. This document, purely theoretical, is in almost good shape, and it seems to be complete since the numbered pages correspond to what Wynn indicates in {26}.

Equations in field extensions
This theoretical document, without title, contains 37 pages numbered from 1 to 37.
Probably it is also intended to follow the work {11}, since inside there are references to that theory. It is perhaps what Wynn called Last notes (see {26}).

Factorisations of a triangular matrix
This document contains 43 pages of well written notes with this title. This work is mentioned in {26-28}, together with the notes {1} on Bürmann series. As usual in many of his works, Wynn begins by introducing a plethora of notations which makes the results more difficult to apprehend. However, we need to give some of them.
A bunch of results are then presented but without any explanation on their purpose. However, it seems that the final goal is the following one: Notation. A being a lower triangular n × n matrix with elements a i,j ∈ K, a(i, k; j) is, for (i, k : n | j : min(n − i + 1, n − k + 1), the jth order determinant formed from the array whose τth row (τ : j) contains the elements a i+τ−1,k+ν−1 (ν : j). Also, a(i, k; 0) = 1 (i : Theorem. Let a i,j (i : n | j : i) be the elements of a nonsingular lower triangular matrix, and c i,j (i : n | j : i) those of its inverse. Then Proof. The formula is obtained from the system of equations The result is correct when j = 1. A short inductive proof based upon the use of the formula, suffices to show that the result is true as stated.
The last page of the document is a list of possible applications, not all of them being intelligible. But Wynn also mentions: Examine band matrix decomposition for a i,j deriving from differential equation; Partial sum transformation

Unpublished Typewritten Documents
We found two typewritten documents which seem to be ready for publication. They are reproduced in the website. constant (slowly varying), and −ψ(µ) = f (µ)J(µ) + f (µ)J (µ), and he assumes that the sequence (e ν ) such that ψ(µ)/ f (µ) = ∑ e ν ( f (µ)) ν is known. Then, he introduces many notations and conditions upon them, but after 20 pages, no clear conclusion follows from this work. In the same document, there is the original review of 11 pages written by Wynn himself of his paper in Reference [106]. It appeared, shortened, in Zentralblatt MATH, as Zbl 0531.40002. Wynn almost surely put these pages together since they both concern the computation of integrals of the form ∞ µ ψ(µ ) dµ .

Extraction of totally monotone sequences from convergents of a continued fraction
This document is a partial copy (pages 31-35) of original notes. There are some pages on an attempt to extract a totally oscillating sequence from the successive convergents of a continued fraction If ∀i, a i = ( − 1) with ∈ (1, ∞), then C i tends to − 1, and the sequence (C 2i+1 − C 2i+3 ) is totally monotone. No other interesting result can be extracted from this partial document. Then, there are 2 pages of rough notes on "Auxiliary sequence transformation before application of the ε-algorithm", which could be related to what precedes, and 3 pages on "Transformation of monotonic sequences by means of the ε-algorithm", which seem to be independent of the topic.

Interpolation Theory
This is an unfinished document of 60 pages, not numbered but in a good shape, entitled Interpolation theory. It is probably a part of a monograph (an important project) referenced in several documents (see {26,28}). Wynn wrote: The following notes concern interpolation in a field by the use of polynomials and rational functions, the interpolatory argument values being assumed to be discrete. The notes serve both as a basis for subsequent more general theory concerning interpolation in the presence of confluent arguments and as a framework for the theory of the transformation of Schweins' series. I. Notations, definitions and classes of matrices and mappings II. Formulae from the calculus of finite differences III. Matrix formulations of finite difference formulae IV. Interpolation by the use of polynomials and rational functions Only the part I. has been found. It is a list of notations and it is too vague to be of interest. A Schweins' series is related to minor identities for quasi-determinants of noncommutative matrices; see, for example, Reference [154]. Schweins obtained it in 1825 for the quotient of two n-rowed determinants which differ only in one column [155].

Linear equations in a commutative ring
This untitled document of 38 pages (numbered from the 9th one) begins by a first list of 7 topics he wanted to treat. Then, 8 pages that look like two tables of contents with several topics covered follow. After that, these items are detailed and, hopefully, Wynn also inserted in front of them the number of the pages where each subitem can be found. In the first page, Wynn inserted also a sketch of the following notes {23}. The first section begins by 1. Linear equations in a commutative ring. It is proposed to study the systems of solutions of equations of the form bx = a mod I where, W being a prescribed commutative ring, a, b, x ∈ W, and I being a system of numbers in W, the above equation is to be interpreted in the sense that u ∈ I exists for which bx − a = u.
The other main sections are: 2. Systems of numbers 3. Divisibility 4.
The ideal reduction A complete analysis of this first section, also if it looks quite complete, seems to be difficult.

Mapping
This only theoretical handwritten draft of 81 pages (there is a page 18 followed by page 18') and is about what Wynn named Aggregates. A sketch of these notes can be found in {22}. It is quite difficult to understand what it is about since the manuscript is not entirely well written, and seems not to have been finished. Let R, S ⊆ W. A nonvoid aggregate σ in (R, S) is composed of a nonvoid set I(σ) ⊆ R and a system of nonvoid sets M(σ, a) ⊆ S defined for each a ∈ I(σ). Then, other definitions and properties follow. Parts of the document are on mappings and classes of aggregates. There is no indication of the purpose of this work, and, searching on internet, nothing similar or related was found.

Factor relationships of the form c = Ed
This is a long theoretical text of 115 pages, entitled Factor relationships of the form c = Ed, where c, E, and d are mapping systems. At the beginning, there is a sketch of the topics that are developed in the sequel. After introducing a lot of notations, Wynn treated the existence of domains of constancy and intersection mapping systems. Then, there is a section on invariance of spaces with respect to sequence rearrangement, one on the properties of domains of constancy, another one on properties of prequotient spaces with respect to domains of constancy and intersection systems, a section giving complete factorisation results, and one on spaces of ordered pairs of mapping systems.

Curriculum vitae of Peter Wynn
This document is a typewritten curriculum vitae by Wynn himself (born 1 September 1931; Hoddesdon, Herts (U.K.)). Then, he reports his education, his professional experience from 1952 to 1980 (thus, it does not contain his positions after that). In this curriculum, we can see that he occupied several positions in different countries in the World. Moreover, we were able to confirm his date and place of birth with certainty. This curriculum contains a bibliography with publications from 1956 to 1981. In this bibliography, Wynn mentioned two references that we did not know: a paper On rational approximations to the exponential function, 1981, to appear, which was found in the boxes and is reproduced at the end of this paper, and a monograph, Numal in

Lists of documents
In each of the two opened boxes, we found a list of documents (one of them with the title Duplicate list. They are probably the lists of documents (xerox copies, handwritten notes and papers, reviews, and so on) that Wynn holds, and that he inserted into the boxes, since we also found a (sometime) partial correspondence mixed with the documents found and analyzed in this paper. There are plenty of what Wynn named "rough notes". Among these rough notes, we find the ones showed in Figure 8, that lead us to think he was looking for a determinantal expression for the even vectors obtained by the vector ε-algorithm.
Since the text could be difficult to read, we reproduce below what Wynn wrote: Thus, one can wonder if Wynn was trying to express the vectors computed by his vector ε-algorithm as ratios of some kind of (generalized) determinants. Such a result would be of great interest for the understanding of this algorithm, but, unfortunately, it has not been obtained yet. In any case, these rough notes are difficult to read and understand!

List of activities
This is a list of activities that have to be accomplished: construct notes and Resumé on different topics. The activities are structured in subitems of an itemized list (from (A) to (J)). Are they separated subjects or the skeleton of a book or for other papers? It is difficult to say. The itemized main topics are: Preparatory algebraic theory; Polynomials and rational functions; Interpolation theory; Moment problem and stability functions; Anti-derivative theory; An array of functions; Functions defined by an integral transform; the Euler-Maclaurin and Boole series; Sequence transformations (a subtopic is about auxiliary transformation before application of the ε-algorithm, a problem addressed in Reference [156]). The item (J) Low priority seems to contains activities having another scope (documentation for Algol programmes, preparation of talks, and another Resumé). This document was probably written after 1981 since Reference [109] is mentioned.

List of projects
There are 3 pages containing a list of projects. Some of them were probably realized (those preceded by a D) and some of them, with an O, perhaps meaning Ongoing. But we are not sure that some other remained in the state of projects, although some of them are preceded by a D. Among this list, we find: theory of stratified commutative ring, functional interpolation, variants of the remainder terms in the Euler Maclaurin and Boole series, auxiliary sequence transformation before application of the ε-algorithm (such a pre-processing or preconditioning as already evoked in References [56,72]), extension of determinantal identities and algorithmic recursions to noncommutative and nonassociative forms by use of linear algebraic equations, numerical experiments in nonassociative algebras, numerical experiments in optimisation, vector ε-algorithm, etc., the analytic continuation of functions defined by an integral transform. And, finally, there is a mention of a book on the ε-algorithm.

Bibliography of various authors
It is a list of references by various authors, probably a list of Xerox copied papers. SSR Ser.A7, 6-8 (1984), also in Russian. This is not strange since Wynn knew Russian (and also German, Duch, Spanish, and some French) and translated two books in this language [110,111]. There are also four handwritten reviews written for zbMATH. The first two are related to the papers (Zbl 0579.40001, Zbl 0571.40003) on continued fractions previously indicated, and the last one (Zbl 0554.65003) to a paper by Naoyuki Tokuda (A new application of Lagrange-Bürmann expansions. I. General principle, Z. Angew. Math. Phys. 34, 697-727 (1983)). Wynn, during his stay in Mexico, was a very active mathematical reviewer for zbMATH (145 reviews) by analyzing papers and books published from 1971 to 2011, some of them in Russian. The last one, very detailed, contains a review of a paper on rational function interpolation and the ε-algorithm, but it is not present in zbMATH. The first three handwritten reviews have been shortened when published in zbMATH.