Następnie wyliczamy posługując się na przykład metodą eliminacji Gaussa, otrzymujemy co daje
co jest zgodne z przykładami Wolframu z dokładnością do mnożnika licznika i mianownika.
Wypełnianie macierzy Niech N=m+n+2 będzie rozmiarem macierzy A z normalnym indeksem liczonym od 1 do N.
Czyścimy macierz inicjując ją zerami; for (int i = 0; i <= m; i++)
for (int j = 0; j <= i; j++)
if (j<=n) A[i+1, m+j+2] = c[i – j];
A[i+1, i+1] = -1;
for (int i = 0; i<= n – 1; i++) for (int j = 0; j <= n; j++) A[m + i + 2, m + j + 2] = c[m + n – i – j]; ; końcowe b0 = 1 A[m + n + 2, m + 2] = 1;
Przypisy
Bibliografia
, Jerome Soucy Université Laval
Linki zewnętrzne Module for Padé Approximation by John H. Mathews Padé Approximants by Oleksandr Pavlyk, The Wolfram Demonstrations Project.
Kategoria:Metody numeryczne Kategoria:Funkcje wymierne