Post by Stachu Chebelf(z)=sin(z)/(z*(exp(z)-1))
Rozwini?cie w szereg Taylora sin(z)/z oraz exp(z)-1 nie stanowi problemu.
Jak to teraz posk?ada? (za przeproszeniem) do kupy?
exp(z) - 1 = z*((exp(z) - 1)/z) = z*h(z)
szereg Taylora h jest latwy do wylizenia. Teraz
f(z) = z^(-2)*(sin(z)/h(z))
z^(-2) to proste przeindeksowanie. Pierwszy wspolczynnik Taylora
h jest niezerowy, wiec dzielenie to w miare prosta rekursja.
Doklaniej, jak
g(z) = \sum a_nz^n
h(z) = \sum b_nz^n
i b_0 jest niezorowe, to
g(z)/h(z) = \sum c_nz^n
gdzie c_n spelnia
a_n = \sum_{k=0}^n b_kc_{n-k} = b_0c_n + \sum_{k=1}^nb_kc_{n-k}
czyli
c_n = (a_n - \sum_{k=1}^nb_kc_{n-k})/b_0
Jak masz juz obliczone c_k dla k < n to mozesz obliczac prawa
strone i dostaniesz c_n.
Ta metoda jest uzywana do komuterowego obliczania c_n.
--
Waldek Hebisch