рдпрд╣ рдХреИрд▓рдХреБрд▓реЗрдЯрд░ рдХреНрдпрд╛ рдХрд░рддрд╛ рд╣реИ
рдпрд╣ рдЯреВрд▓ \(y' = F(x, y)\) рд░реВрдк рдХреЗ рдХрд┐рд╕реА рдкреНрд░рдердо-рдХреНрд░рдо рд╕рд╛рдзрд╛рд░рдг рдЕрд╡рдХрд▓ рд╕рдореАрдХрд░рдг рдХреЛ рджреА рдЧрдИ рдкреНрд░рд╛рд░рдВрднрд┐рдХ рд╢рд░реНрдд \(y(x_0) = y_0\) рдХреЗ рд╕рд╛рде, \(x_0\) рд╕реЗ \(x_n\) рддрдХ рдХреЗ рдЕрдВрддрд░рд╛рд▓ рдкрд░ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд░реВрдк рд╕реЗ рд╣рд▓ рдХрд░рддрд╛ рд╣реИред рдЗрд╕рдореЗрдВ рдХреНрд▓рд╛рд╕рд┐рдХрд▓ рдЪреМрдереЗ-рдХреНрд░рдо рд░реБрдВрдЧреЗ-рдХреБрдЯреНрдЯрд╛ рд╡рд┐рдзрд┐ (RK4) рдХрд╛ рдЙрдкрдпреЛрдЧ рд╣реЛрддрд╛ рд╣реИ, рдЬреЛ рд╕рдВрдЦреНрдпрд╛рддреНрдордХ рд╡рд┐рд╢реНрд▓реЗрд╖рдг рдореЗрдВ рд╕рдмрд╕реЗ рдЕрдзрд┐рдХ рдкреНрд░рдЪрд▓рд┐рдд рдФрд░ рднрд░реЛрд╕реЗрдордВрдж рдПрдХрд▓-рдЪрд░рдг рдЗрдВрдЯреАрдЧреНрд░реЗрдЯрд░реЛрдВ рдореЗрдВ рд╕реЗ рдПрдХ рд╣реИред рдкрд░рд┐рдгрд╛рдо рдХреЗ рд░реВрдк рдореЗрдВ рдЖрдкрдХреЛ \((x_i, y_i)\) рдмрд┐рдВрджреБрдУрдВ рдХреА рдПрдХ рд╕рд╛рд░рдгреА рдорд┐рд▓рддреА рд╣реИ рдЬреЛ рд╡рд╛рд╕реНрддрд╡рд┐рдХ рд╣рд▓ рдХрд╛ рд╕рдиреНрдирд┐рдХрдЯ рдорд╛рди рджреЗрддреА рд╣реИ, рд╕рд╛рде рд╣реА рдЕрдВрддрд┐рдо рдорд╛рди \(y(x_n)\) рднреАред рдпрд╣ рд╡рд┐рд╢реБрджреНрдз рдЧрдгрд┐рддреАрдп рдЯреВрд▓ рд╣реИ, рдЗрд╕рдХрд╛ рдХрд┐рд╕реА рджреЗрд╢ рдпрд╛ рдЗрдХрд╛рдИ рд╕реЗ рдХреЛрдИ рд╕рдВрдмрдВрдз рдирд╣реАрдВ рд╣реИред
рдЗрд╕рдХрд╛ рдЙрдкрдпреЛрдЧ рдХреИрд╕реЗ рдХрд░реЗрдВ
рджрд╛рд╣рд┐рдиреА рдУрд░ рдХрд╛ рдкрдж \(F(x,y)\) рдХреЛ x рдФрд░ y рдореЗрдВ рдПрдХ рд╡реНрдпрдВрдЬрдХ рдХреЗ рд░реВрдк рдореЗрдВ рдбрд╛рд▓реЗрдВ (рдЙрджрд╛рд╣рд░рдг рдХреЗ рд▓рд┐рдП 1-y^2, x+y, x*y, рдпрд╛ sin(x)+y)ред рд╕рдорд░реНрдерд┐рдд рдСрдкрд░реЗрдЯрд░ рд╣реИрдВ + - * / ^ рдФрд░ рдлрд▓рди рдЬреИрд╕реЗ sin, cos, tan, exp, log, ln, sqrt, abs, tanh, рд╕рд╛рде рд╣реА рд╕реНрдерд┐рд░рд╛рдВрдХ pi рдФрд░ eред рдкреНрд░рд╛рд░рдВрднрд┐рдХ рдорд╛рди \(x_0\), рдкреНрд░рд╛рд░рдВрднрд┐рдХ \(y_0\), рдЕрдВрддрд┐рдо \(x_n\) рд╕реЗрдЯ рдХрд░реЗрдВ, рдФрд░ рдЪреБрдиреЗрдВ рдХрд┐ рдХрд┐рддрдиреЗ рд╕рдорд╛рди рдЙрдкрд╡рд┐рднрд╛рдЬрди \(n\) рд▓реЗрдиреЗ рд╣реИрдВред рдЬрд┐рддрдиреЗ рдЕрдзрд┐рдХ рдЙрдкрд╡рд┐рднрд╛рдЬрди рд╣реЛрдВрдЧреЗ, рд╕рдЯреАрдХрддрд╛ рдЙрддрдиреА рд╣реА рдЕрдзрд┐рдХ рд╣реЛрдЧреА, рдХреНрдпреЛрдВрдХрд┐ RK4 рдХреА рд╡реИрд╢реНрд╡рд┐рдХ рддреНрд░реБрдЯрд┐ \(h^4\) рдХреЗ рдЕрдиреБрдкрд╛рдд рдореЗрдВ рдШрдЯрддреА рд╣реИред
рд╕реВрддреНрд░ рдХреА рд╡реНрдпрд╛рдЦреНрдпрд╛
рдЕрдВрддрд░рд╛рд▓ рдХреЛ \(n\) рд╕рдорд╛рди рдЪрд░рдгреЛрдВ рдореЗрдВ рдмрд╛рдБрдЯрд╛ рдЬрд╛рддрд╛ рд╣реИ, рдЬрд┐рдирдореЗрдВ рд╕реЗ рдкреНрд░рддреНрдпреЗрдХ рдХреА рдЪреМрдбрд╝рд╛рдИ \(h = (x_n - x_0)/n\) рд╣реЛрддреА рд╣реИред рд╣рд░ рдЪрд░рдг рдкрд░ RK4 рдврд╛рд▓ (slope) рдХрд╛ рдЪрд╛рд░ рдмрд╛рд░ рдорд╛рди рдирд┐рдХрд╛рд▓рддрд╛ рд╣реИ: рдПрдХ рдмрд╛рд░ рд╢реБрд░реБрдЖрдд рдореЗрдВ (\(k_1\)), рджреЛ рдмрд╛рд░ рдордзреНрдп-рдмрд┐рдВрджреБ рдкрд░ (\(k_2, k_3\)), рдФрд░ рдПрдХ рдмрд╛рд░ рдЕрдВрдд рдореЗрдВ (\(k_4\))ред рдЕрдЧрд▓рд╛ рдорд╛рди рдПрдХ рднрд╛рд░рд┐рдд рдФрд╕рдд рд╣реЛрддрд╛ рд╣реИ, $$y_{n+1} = y_n + \frac{h}{6}\left(k_1 + 2k_2 + 2k_3 + k_4\right)$$ рдЗрд╕рд╕реЗ рдЪреМрдереЗ рдХреНрд░рдо рддрдХ рдХреА рддреНрд░реБрдЯрд┐ рдХреЗ рдкрдж рдХрдЯ рдЬрд╛рддреЗ рд╣реИрдВ, рдЬрд┐рд╕рд╕реЗ рд╕реНрдерд╛рдиреАрдп рдХрдЯреМрддреА рддреНрд░реБрдЯрд┐ \(O(h^5)\) рдФрд░ рд╡реИрд╢реНрд╡рд┐рдХ рддреНрд░реБрдЯрд┐ \(O(h^4)\) рдорд┐рд▓рддреА рд╣реИред
рд╣рд▓ рдХрд┐рдпрд╛ рдЧрдпрд╛ рдЙрджрд╛рд╣рд░рдг
\(x_0 = 0\), \(y_0 = 0\), \(x_n = 1\) рдФрд░ \(n = 10\) (\(h = 0.1\)) рдХреЗ рд╕рд╛рде \(y' = 1 - y^2\) рдХреЛ рд╣рд▓ рдХрд░реЗрдВред рд╕рдЯреАрдХ рд╣рд▓ \(y = \tanh(x)\) рд╣реИред рдкрд╣рд▓рд╛ RK4 рдЪрд░рдг \(y_1 = 0.0996679\) рджреЗрддрд╛ рд╣реИ, рдЬреЛ \(\tanh(0.1) = 0.0996680\) рд╕реЗ рдореЗрд▓ рдЦрд╛рддрд╛ рд╣реИред рд╕рднреА рджрд╕ рдЪрд░рдгреЛрдВ рдХреЗ рдмрд╛рдж \(y(1) = 0.7615942\) рдорд┐рд▓рддрд╛ рд╣реИ, рдЬреЛ \(\tanh(1) = 0.7615942\) рд╕реЗ рд╕рд╛рдд рдЕрдВрдХреЛрдВ рддрдХ рдореЗрд▓ рдЦрд╛рддрд╛ рд╣реИред
рдЕрдХреНрд╕рд░ рдкреВрдЫреЗ рдЬрд╛рдиреЗ рд╡рд╛рд▓реЗ рдкреНрд░рд╢реНрди
RK4, рдСрдпрд▓рд░ рд╡рд┐рдзрд┐ рд╕реЗ рдмреЗрд╣рддрд░ рдХреНрдпреЛрдВ рд╣реИ? рдСрдпрд▓рд░ рд╡рд┐рдзрд┐ рд╣рд░ рдЪрд░рдг рдореЗрдВ рдХреЗрд╡рд▓ рдПрдХ рдврд╛рд▓ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рддреА рд╣реИ (рддреНрд░реБрдЯрд┐ \(O(h)\))ред RK4 рдЪрд╛рд░ рдврд╛рд▓реЛрдВ рдХрд╛ рдЙрдкрдпреЛрдЧ рдХрд░рдХреЗ рдЙрдирдХрд╛ рдФрд╕рдд рдирд┐рдХрд╛рд▓рддреА рд╣реИ, рдФрд░ рдЙрд╕реА рдЪрд░рдг-рдЖрдХрд╛рд░ рдкрд░ \(O(h^4)\) рд╕рдЯреАрдХрддрд╛ рд╣рд╛рд╕рд┐рд▓ рдХрд░рддреА рд╣реИ тАФ рдЗрд╕рд▓рд┐рдП рдХрд┐рд╕реА рдирд┐рд░реНрдзрд╛рд░рд┐рдд рд╕рдЯреАрдХрддрд╛ рдХреЗ рд▓рд┐рдП рдХрд╣реАрдВ рдХрдо рдЪрд░рдг рдЪрд╛рд╣рд┐рдПред
рдореБрдЭреЗ рдХрд┐рддрдиреЗ рдЪрд░рдг рдЪреБрдирдиреЗ рдЪрд╛рд╣рд┐рдП? 50 рд╕реЗ рд╢реБрд░реВ рдХрд░реЗрдВред рдпрджрд┐ рд╣рд▓ рд╕рд╣рдЬ (smooth) рд╣реИ, рддреЛ рдЖрдорддреМрд░ рдкрд░ рдЗрддрдирд╛ рд╣реА рдкрд░реНрдпрд╛рдкреНрдд рд╣реЛрддрд╛ рд╣реИ; рддреЗрдЬрд╝реА рд╕реЗ рдмрджрд▓рдиреЗ рд╡рд╛рд▓реА рдпрд╛ рд▓рдЧрднрдЧ рд╕реНрдЯрд┐рдлрд╝ (near-stiff) рд╕рдорд╕реНрдпрд╛рдУрдВ рдХреЗ рд▓рд┐рдП рдЗрд╕реЗ 100, 200 рдпрд╛ 500 рддрдХ рдмрдврд╝рд╛рдПрдБред
рдЕрдЧрд░ рдореБрдЭреЗ NaN рдпрд╛ Infinity рдорд┐рд▓реЗ рддреЛ рдХреНрдпрд╛ рдХрд░реЗрдВ? рд╕рдВрднрд╡ рд╣реИ рдХрд┐ рд╣рд▓ рдЕрдкрд╕рд░рд┐рдд (diverge) рд╣реЛ рдЧрдпрд╛ рд╣реЛ, рдпрд╛ \(F(x,y)\) рдиреЗ рдХрд┐рд╕реА рдЕрдорд╛рдиреНрдп рд╕рдВрдХреНрд░рд┐рдпрд╛ рдХрд╛ рд╕рд╛рдордирд╛ рдХрд┐рдпрд╛ рд╣реЛ (рдЬреИрд╕реЗ рдЛрдгрд╛рддреНрдордХ рд╕рдВрдЦреНрдпрд╛ рдХрд╛ log рдпрд╛ рд╢реВрдиреНрдп рд╕реЗ рднрд╛рдЧ)ред рд╡реНрдпрдВрдЬрдХ рдХреА рдЬрд╛рдБрдЪ рдХрд░реЗрдВ рдФрд░ рдЫреЛрдЯрд╛ рдЕрдВрддрд░рд╛рд▓ рдпрд╛ рдЕрдзрд┐рдХ рдЪрд░рдг рдЖрдЬрд╝рдорд╛рдПрдБред