рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ
рдпрд╣ рдЯреВрд▓ n-рдмрд┐рдВрджреБ рд╕рд╛рдорд╛рдиреНрдпреАрдХреГрдд рдЧреЙрд╕-рд▓рд╛рдЧреБрдПрд░ рдХреНрд╡рд╛рдбреНрд░реЗрдЪрд░ рдХреЗ рд▓рд┐рдП рдиреЛрдбреНрд╕ (рдПрдмреНрд╕рд┐рд╕рд╛) \(x_i\) рдФрд░ рд╡реЗрдЯреНрд╕ \(w_i\) рдХреА рдЧрдгрдирд╛ рдХрд░рддрд╛ рд╣реИред рдпрд╣ рд╢реБрджреНрдз рдЧрдгрд┐рдд рдХрд╛ рдПрдХ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╕рдорд╛рдХрд▓рди (numerical integration) рдЯреВрд▓ рд╣реИ рдЬреЛ рд╣рд░ рдЬрдЧрд╣ рдПрдХ рдЬреИрд╕рд╛ рдХрд╛рдо рдХрд░рддрд╛ рд╣реИред рдпрд╣ рдирд┐рдпрдо рдЕрд░реНрдз-рдЕрдирдВрдд рдЕрдВрддрд░рд╛рд▓ [0, тИЮ) рдкрд░ рдЙрди рд╕рдорд╛рдХрд▓рдиреЛрдВ рдХрд╛ рд╕рдиреНрдирд┐рдХрдЯрди (approximation) рдХрд░рддрд╛ рд╣реИ рдЬрд┐рдирдореЗрдВ рднрд╛рд░ рдлрд▓рди (weight function) \(x^{\alpha}e^{-x}\) рд╣реЛрддрд╛ рд╣реИ:
$$\int_{0}^{\infty} x^{\alpha} e^{-x} f(x)\, dx \approx \sum_{i=1}^{n} w_i\, f(x_i)$$рдиреЛрдбреНрд╕ рд╕рд╛рдорд╛рдиреНрдпреАрдХреГрдд рд▓рд╛рдЧреБрдПрд░ рдмрд╣реБрдкрдж \(L_n^{(\alpha)}(x)\) рдХреЗ рдзрдирд╛рддреНрдордХ рдореВрд▓ (zeros) рд╣реЛрддреЗ рд╣реИрдВ, рдФрд░ рдпрд╣ рдирд┐рдпрдо рддрдм рдкреВрд░реНрдгрддрдГ рд╕рдЯреАрдХ рд░рд╣рддрд╛ рд╣реИ рдЬрдм f рдЕрдзрд┐рдХрддрдо \(2n-1\) рдШрд╛рдд рдХрд╛ рдмрд╣реБрдкрдж рд╣реЛред
рдЗрд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ
рдХреНрд░рдо n (рдмрд┐рдВрджреБрдУрдВ рдХреА рд╕рдВрдЦреНрдпрд╛, 2 рд╕реЗ 100) рдЪреБрдиреЗрдВ, рдШрд╛рддрд╛рдВрдХ рдкреИрд░рд╛рдореАрдЯрд░ ╬▒ рджрд░реНрдЬ рдХрд░реЗрдВ (тИТ1 рд╕реЗ рдмрдбрд╝реА рдХреЛрдИ рднреА рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╕рдВрдЦреНрдпрд╛; рдХреНрд▓рд╛рд╕рд┐рдХ рдЧреЙрд╕-рд▓рд╛рдЧреБрдПрд░ рдирд┐рдпрдо рдореЗрдВ \(\alpha = 0\) рд╣реЛрддрд╛ рд╣реИ), рдФрд░ рддрдп рдХрд░реЗрдВ рдХрд┐ рдЖрдк рдХрд┐рддрдиреЗ рд╕рд╛рд░реНрдердХ рджрд░реНрд╢рд╛рдП рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рдЕрдВрдХ (display digits) рдЪрд╛рд╣рддреЗ рд╣реИрдВред рдкрд░рд┐рдгрд╛рдо рдореЗрдВ рд╣рд░ рдиреЛрдб рдФрд░ рдЙрд╕рд╕реЗ рдЬреБрдбрд╝рд╛ рд╡реЗрдЯ \(x_i\) рдХреЗ рдмрдврд╝рддреЗ рдХреНрд░рдо рдореЗрдВ рджрд┐рдЦрддрд╛ рд╣реИ, рд╕рд╛рде рд╣реА рдПрдХ рдЕрдВрддрд░реНрдирд┐рд╣рд┐рдд рд╕реНрд╡-рдЬрд╛рдБрдЪ (self-check) рднреАред
рд╕реВрддреНрд░ рдФрд░ рд╡рд┐рдзрд┐
рд╣рд░ рд╡реЗрдЯ рдЗрд╕ рдмрдВрдж рд╕реВрддреНрд░ рдХрд╛ рдкрд╛рд▓рди рдХрд░рддрд╛ рд╣реИ: $$w_i = \frac{\Gamma(n+\alpha+1)\cdot x_i}{n!\cdot\left[(n+1)L_{n+1}^{(\alpha)}(x_i)\right]^{2}}$$ рдЖрдВрддрд░рд┐рдХ рд░реВрдк рд╕реЗ рд╣рдо рд╕рдорддреБрд▓реНрдп рдФрд░ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд░реВрдк рд╕реЗ рд╕реНрдерд┐рд░ рдЧреЛрд▓рдм-рд╡реЗрд▓реНрд╢ (Golub-Welsch) рд╡рд┐рдзрд┐ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреЗ рд╣реИрдВ: рд╡рд┐рдХрд░реНрдг \(a_k = 2k+\alpha+1\) рдФрд░ рдЙрдк-рд╡рд┐рдХрд░реНрдгреЛрдВ \(b_k = \sqrt{k(k+\alpha)}\) рд╡рд╛рд▓рд╛ рд╕рдордорд┐рдд рддреНрд░рд┐рд╡рд┐рдХрд░реНрдгреАрдп (tridiagonal) рдпрд╛рдХреЛрдмреА рдореИрдЯреНрд░рд┐рдХреНрд╕ рдмрдирд╛рдПрдБред рдЗрд╕рдХреЗ рдЖрдЗрдЧреЗрдирдорд╛рди (eigenvalues) рд╣реА рдиреЛрдбреНрд╕ рд╣реЛрддреЗ рд╣реИрдВ, рдФрд░ рд╣рд░ рд╡реЗрдЯ \(\mu_0\cdot(\text{рдкрд╣рд▓реЗ рдЖрдЗрдЧреЗрдирд╡реЗрдХреНрдЯрд░ рдШрдЯрдХ})^{2}\) рдХреЗ рдмрд░рд╛рдмрд░ рд╣реЛрддрд╛ рд╣реИ, рдЬрд╣рд╛рдБ \(\mu_0 = \Gamma(\alpha+1)\) рд╢реВрдиреНрдпрд╡рд╛рдБ рдЖрдШреВрд░реНрдг (zeroth moment) рд╣реИред рдпрд╣ рдмрдбрд╝реЗ рдХреНрд░рдордЧреБрдгрд┐рддреЛрдВ (factorials) рд╕реЗ рд╣реЛрдиреЗ рд╡рд╛рд▓реЗ рдУрд╡рд░рдлрд╝реНрд▓реЛ рд╕реЗ рдмрдЪрд╛рддрд╛ рд╣реИред
рд╣рд▓ рдХрд┐рдпрд╛ рд╣реБрдЖ рдЙрджрд╛рд╣рд░рдг
n = 2, ╬▒ = 0 рдХреЗ рд▓рд┐рдП: \(L_2^{(0)}(x) = (x^2-4x+2)/2\), рдЕрддрдГ рдореВрд▓ рд╣реИрдВ \(x = 2 \pm \sqrt{2}\), рдЬрд┐рд╕рд╕реЗ \(x_1 = 0.5857864\) рдФрд░ \(x_2 = 3.4142136\) рдорд┐рд▓рддреЗ рд╣реИрдВред рд╡реЗрдЯреНрд╕ рд╣реИрдВ $$w_1 = \frac{2+\sqrt{2}}{4} = 0.8535534 \qquad w_2 = \frac{2-\sqrt{2}}{4} = 0.1464466$$ рдЗрдирдХрд╛ рдпреЛрдЧ \(1 = \Gamma(1)\) рд╣реЛрддрд╛ рд╣реИ, рдЬреЛ рдкрд░рд┐рдгрд╛рдо рдХреА рдкреБрд╖реНрдЯрд┐ рдХрд░рддрд╛ рд╣реИред
рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рдкреНрд░рд╢реНрди
╬▒ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ? рдпрд╣ рднрд╛рд░ \(x^{\alpha}\) рддрдп рдХрд░рддрд╛ рд╣реИ; \(\alpha = 0\) рдорд╛рдирдХ рдЧреЙрд╕-рд▓рд╛рдЧреБрдПрд░ рджреЗрддрд╛ рд╣реИ, рдЬрдмрдХрд┐ \(\alpha > 0\) рднрд╛рд░ рдХреЛ рдореВрд▓ рдмрд┐рдВрджреБ (origin) рд╕реЗ рджреВрд░ рдЕрдзрд┐рдХ рдорд╣рддреНрд╡ рджреЗрддрд╛ рд╣реИред рдпрд╣ тИТ1 рд╕реЗ рдмрдбрд╝рд╛ рд╣реЛрдирд╛ рдЬрд╝рд░реВрд░реА рд╣реИред
рдпрд╣ рдХрд┐рддрдирд╛ рд╕рдЯреАрдХ рд╣реИ? n-рдмрд┐рдВрджреБ рдирд┐рдпрдо \(2n-1\) рдШрд╛рдд рддрдХ рдХреЗ рдмрд╣реБрдкрджреЛрдВ рдХрд╛ рдкреВрд░реНрдгрддрдГ рд╕рдЯреАрдХ рд╕рдорд╛рдХрд▓рди рдХрд░рддрд╛ рд╣реИ; рдЪрд┐рдХрдиреЗ (smooth) рдлрд▓рди рддреЗрдЬрд╝реА рд╕реЗ рдЕрднрд┐рд╕рд░рд┐рдд (converge) рд╣реЛрддреЗ рд╣реИрдВред
рдЖрдЙрдЯрдкреБрдЯ рдХреА рдЬрд╛рдБрдЪ рдХреИрд╕реЗ рдХрд░реВрдБ? рд╕рднреА рд╡реЗрдЯреНрд╕ рдХрд╛ рдпреЛрдЧ рд╣рдореЗрд╢рд╛ \(\Gamma(\alpha+1)\) рдХреЗ рдмрд░рд╛рдмрд░ рд╣реЛрддрд╛ рд╣реИ, рдЬреЛ рд╢реВрдиреНрдпрд╡реЗрдВ рдЖрдШреВрд░реНрдг (zeroth moment) рд╡рд╛рд▓реА рдкрдВрдХреНрддрд┐ рдореЗрдВ рджрд┐рдЦрд╛рдпрд╛ рдЬрд╛рддрд╛ рд╣реИред