MHB Calculations - absolute stability

AI Thread Summary
The discussion focuses on determining the order of accuracy and absolute stability for a specific Runge-Kutta method applied to a system of differential equations. The method is identified as having an order of accuracy of 2. The region of absolute stability is defined as S = { z ∈ ℂ: |(z²/2) + z + 1| < 1}, which is crucial for determining the appropriate step size, h, for stability in calculations. Participants inquire about the definition of absolute stability and the method for finding it, referencing the stability function r and its relationship to the method's tableau. Understanding these concepts is essential for ensuring accurate and stable numerical solutions.
evinda
Gold Member
MHB
Messages
3,741
Reaction score
0
Hello! (Wave) The following Runge-Kutta method is given.

$$ \begin{array}{c|ccccc}
\tau_1 =0 & a_{11}=0 & a_{12} = 0\\
\tau_2 =\frac{5}{2} & a_{21} = \frac{5}{2} & a_{22} = 0\\
\hline
& b_1 = \frac{4}{5} & b_2 = \frac{1}{5} & \
\end{array} $$

I have to determine the order of accuracy and I found that it is $2$.


Then if the method is applied at the system $\\y_1'=-80y_1+20y_2 \\
y_2'=20y_1-80y_2$ what $h$ should we take so that the calculations get done with absolute stability?

How can we find such an $ h$?
I found that the region of absolute stability is $S=\{ z \in \mathbb{C}: |\frac{z^2}{2}+z+1|<1\}$. Does this help?
 
Mathematics news on Phys.org
Hi! (Mmm)

What do your notes say absolute stability is exactly? And how to find it? (Wondering)

Wiki says that the stability function $r$ is:
$$r(z) = 1 + z b^T (I-zA)^{-1} e = \frac{\det(I-zA+zeb^T)}{\det(I-zA)}$$
where $e$ stands for the vector of ones, and $b^T$ is the bottom row of the tableau. (Thinking)
 
Insights auto threads is broken atm, so I'm manually creating these for new Insight articles. In Dirac’s Principles of Quantum Mechanics published in 1930 he introduced a “convenient notation” he referred to as a “delta function” which he treated as a continuum analog to the discrete Kronecker delta. The Kronecker delta is simply the indexed components of the identity operator in matrix algebra Source: https://www.physicsforums.com/insights/what-exactly-is-diracs-delta-function/ by...
Suppose ,instead of the usual x,y coordinate system with an I basis vector along the x -axis and a corresponding j basis vector along the y-axis we instead have a different pair of basis vectors ,call them e and f along their respective axes. I have seen that this is an important subject in maths My question is what physical applications does such a model apply to? I am asking here because I have devoted quite a lot of time in the past to understanding convectors and the dual...
Back
Top