Page 16-70
These results indicate that (Δx)
next
= 0.34049…
Function RRKSTEP
This function uses an input list similar to that of function RRK, as well as the
tolerance for the solution, a possible step Δx, and a number (LAST) specifying
the last method used in the solution (1, if RKF was used, or 2, if RRK was used).
Function RRKSTEP returns the same input list, followed by the tolerance, an
estimate of the next step in the independent variable, and the current method
(CURRENT) used to arrive at the next step. Thus, the input stack looks as
follows:
4: {‘x’, ‘y’, ‘f(x,y)’}
3: ε
2: Δx
1: LAST
After running this function, the stack will show the lines:
4: {‘x’, ‘y’, ‘f(x,y)’}
3: ε
2: (Δx)
next
1: CURRENT
Thus, this function is used to determine the appropriate size of a time step
((Δx)
next
) to satisfy the required tolerance, and the method used to arrive at that
result (CURRENT).
The following screen shots show the RPN stack before and after application of
function RRKSTEP: