MATHEMATICA TUTORIAL, Part 1.3: Heun Methods
www.cfm.brown.edu › am33 › Mathematicap n + 1 = y n + h f ( x n, y n), y n + 1 = y n + h 2 [ f ( x n, y n) + f ( x n + 1, p n + 1)], n = 0, 1, 2, …. Therefore, one of the versions on the Heun method in Mathematica is as follows. Clear [y] f [x_, y_] := x^2 + y^2. y [0] = 1; h:=0.1; Do [k1 = h f [h n, y [n]]; k2 = h f [h (n + 1), y [n] + k1];
Numerical Methods--Heun's Method
calculuslab.deltacollege.edu/ODE/7-C-2/7-C-2-h.htmlThat's what we'll do with Heun's method! We'll use Euler's method to roughly estimate the coordinates of the next point in the solution, and once we have this information, we'll re-predict (or correct) our original estimate of the location of the next solution point by using the method of averaging the slopes of the left and right tangent lines that we so carefully developed above.
Topic 14.2: Heun's Method (Examples)
ece.uwaterloo.ca › 14IVPs › heunTopic 14.2: Heun's Method (Examples) Introduction Notes Theory HOWTO Examples Engineering Error Questions Matlab Maple Example 1 Given the IVP y (1) ( t) = 1 - t y ( t) with y (0) = 1, approximate y (1) with one step. First, let t0 = 0, y0 = 1, and h = 1. Thus, we calculate K0 = f (0, 1) = 1 K1 = f (1, 2) = -1 ½ ( K0 + K1) = 0