Quote:
Originally Posted by Ether
Yeah; I had already prepared Question 6 back when everyone was using a closed-form analytical solution for Q(t). You beat me to the punch by setting up your numerical integration script.
|
Makes sense.
Quote:
|
Originally Posted by Ether
Question 6 solution:
I think you are showing more decimal places than are warranted for the accuracy of your solution method. |
I never bothered to do any accuracy analysis, so that's probably the case.
After lowering the tolerances and calculating bounds for the error (total error should be less than 1e-7), I got the following revised figures:
New parameters and error analysis:
Relative and absolute tolerances set to 10^-12 (previously 10^-6): The error at each step should not exceed max(10^-12, |x * 10^-12|), where x is the current value of the ODE solution
Number of ODE solver steps: 1022
Maximum element of the state: 60.8668
Upper bound for error: max(10^-12, |60.8668 * 10^-12|) * 1022 = 6.2267e-08