error_answers

MATLAB:

• expsin_deriv.m, defines the right hand side of the EXPSIN ODE.
• expsin_rk1_double.m, estimates the EXPSIN ODE approximation error by comparing rk1 estimates with n and 2n steps.
• expsin_rk1_extrapolate.m, estimates the EXPSIN ODE approximation error by comparing rk1 estimates with n and 2n steps, and using extrapolation.
• expsin_rk1_higher.m, estimates the EXPSIN ODE approximation error by comparing rk1 and rk2 estimates.
• expsin_solution.m, defines the exact solution of the EXPSIN ODE.
• normal_deriv.m, defines the right hand side of the NORMAL ODE.
• normal_rk1_double.m, estimates the NORMAL ODE approximation error by comparing rk1 estimates with n and 2n steps.
• normal_rk1_extrapolate.m, estimates the NORMAL ODE approximation error by comparing rk1 estimates with n and 2n steps, and using extrapolation.
• normal_solution.m, defines the exact solution of the NORMAL ODE.
• rk1.m, a Runge-Kutta ODE solver of order 1 (this is just the Euler method).
• rk2.m, a Runge-Kutta ODE solver of order 2.
• rms.m, evaluates the RMS norm of a vector.

Images:

• expsin_rk1_double_20.png, displays the n=20 expsin_rk1_double results.
• expsin_rk1_double_error_20.png, displays the estimated error for the n=20 expsin_rk1_double results.
• expsin_rk1_extrapolate_20.png, displays the n=20 expsin_rk1_extrapolate results.
• expsin_rk1_extrapolate_error_20.png, displays the estimated error for the n=20 expsin_rk1_extrapolate results.
• expsin_rk1_higher_20.png, displays the n=20 expsin_rk1_higher results.
• expsin_rk1_higher_error_20.png, displays the estimated error for the n=20 expsin_rk1_higher results.
• normal_rk1_double_10.png, displays the n=10 normal_rk1_double results.
• normal_rk1_double_error_10.png, displays the estimated error for the n=10 normal_rk1_double results.
• normal_rk1_extrapolate_10.png, displays the n=10 normal_rk1_extrapolate results.
• normal_rk1_extrapolate_error_10.png, displays the estimated error for the n=10 normal_rk1_extrapolate results.