Intuition The ONE core idea
When several quantities change together and each one's rate of change depends on all the others, we stack them into a single list called a vector and treat that list as one moving object. Then the same four-slope averaging trick (RK4) that works for one quantity works for the whole list, with nothing new invented.
This page assumes zero prior notation. We build every symbol the parent note 4.8.27 leans on, one brick at a time, so that when you read the recipe there, every mark on the page already means something to you.
t and a "quantity that changes"
t is plain time . A quantity like y is a number that has a different value at each moment — think of the height of a swing. We write y ( t ) : "the value of y at time t ".
Picture a single dot sliding along a vertical line as the clock ticks. At t = 0 it sits at some height; a moment later it has moved. That motion is what all of this is about.
Why we need this: the whole topic is about predicting where the dot will be next given where it is now and how fast it is moving.
y ′ (read "y prime")
y ′ means "==how fast y is changing right now=="— its rate of change . If y is a height, y ′ is the vertical speed (up = positive, down = negative).
y ′ is literally a slope
Plot y against time t on a graph. The curve rises and falls. At any single instant, zoom in until the curve looks like a straight line — the steepness (slope) of that tiny line is y ′ . Steep-up = big positive y ′ ; flat = zero; steep-down = negative.
Why we need this: an ODE is a rule that tells us the slope. If we know the slope everywhere, we can march forward and trace the whole curve — that marching is the numerical method.
f ( t , y )
f is a machine : feed it the current time t and current value y , and it hands back the slope y ′ at that moment. So the equation y ′ = f ( t , y ) reads: "the rate of change equals whatever the machine says, given where and when we are."
Worked example Reading a slope-rule out loud
y ′ = − y says "the faster you are high, the faster you fall back" (slope is the negative of your value). y ′ = y says "the higher you are, the faster you climb" — that runaway growth is exactly e t .
Why we need this: an ODE (Ordinary Differential Equation) is exactly one such rule. "Solving" it means finding the moving value y ( t ) whose slope always matches what f demands.
Definition Starting point
t 0 is the start time (often 0 ). y 0 is the known value there . Together y ( t 0 ) = y 0 pins the dot to one spot before it starts moving.
Intuition Why one slope-rule is not enough
The rule f gives the shape of the motion but not where it begins . Many curves have the same slope pattern shifted up or down — like parallel hillsides. The initial value picks which one you are on. A slope-rule plus a starting point is called an IVP (Initial Value Problem) .
h
h is the size of one time-hop we take. Instead of knowing y at every instant, we settle for knowing it at t 0 , t 0 + h , t 0 + 2 h , … — like footprints spaced h apart.
Small h = tiny careful steps = more accurate but slower. The subscript n in y n means "the value after n steps", so y n ≈ y ( t 0 + nh ) .
Why we need this: computers cannot store infinitely many instants, so we walk the curve in discrete hops of length h . Choosing h trades accuracy against work — the theme of Local vs Global Truncation Error .
y
When we track several quantities together — call them y 1 , y 2 , … , y n — we stack them into one column list :
y = ( y 1 , y 2 , … , y n ) T .
The bold y signals "this is a whole list, not one number". The little superscript T (transpose) just means "written as a standing-up column".
list and not separate numbers
Predator and prey, or position and velocity, move together : each one's rate depends on the other. Keeping them in one list lets us say "advance the whole state by one hop" in a single stroke — and, crucially, forces us to advance them at the same time , never one while the other is stale.
The number n is how many quantities there are (the dimension ). n = 1 is the ordinary scalar case; n = 2 is predator–prey or a pendulum.
Definition Vector slope-rule
f is a machine that takes the time t and the whole list y , and returns a list of slopes — one for each quantity:
f ( t , y ) = ( f 1 , f 2 , … , f n ) , y i ′ = f i ( t , y 1 , … , y n ) .
The key point: each f i may look at all the y j 's. That cross-dependence is called coupling .
Common mistake "The equations don't really touch, so solve them one at a time."
Why it feels right: y 1 ′ and y 2 ′ are written on separate lines. Why it's wrong: because f 1 may contain y 2 , freezing y 2 while you advance y 1 feeds the wrong value. Fix: always advance the entire list together , step by step. This is the single most important habit the parent note drills.
Definition Order of an ODE
The order is the highest prime-count. y ′ alone → first-order . y ′′ (the slope's own slope, i.e. acceleration) → second-order .
Intuition Why we only ever need first-order
A second-order rule like y ′′ = g ( t , y , y ′ ) can be repackaged as a first-order system by naming the derivative as a brand-new variable:
u 1 = y , u 2 = y ′ .
Then u 1 ′ = u 2 (the value's rate is the velocity) and u 2 ′ = g ( t , u 1 , u 2 ) (the velocity's rate is the acceleration). One 2nd-order equation becomes two 1st-order ones — see Reducing higher-order ODEs to first-order systems .
Why we need this: it means "RK4 for systems" secretly solves every ODE, of any order, once you reduce it.
Definition Small notation, spelled out
≈ means "approximately equal to " — our hops give near-answers, not exact ones.
∫ t t + h … d s is a definite integral : the exact accumulated change over one hop, the area under the slope curve between t and t + h . RK4 estimates this area by sampling the slope at a few clever points (the same spirit as Simpson's Rule ).
O ( h 4 ) is big-O : "the error shrinks in step with h 4 ." Halve h and the error drops by 2 4 = 16 . See Local vs Global Truncation Error .
2 h , 6 h are just fractions of the step — the half-hops and the weighting divisor in the recipe.
Recall Why
O ( h 4 ) is the headline number
Question: if you halve the step h , roughly how much smaller does RK4's global error get?
Answer ::: About 2 4 = 16 times smaller — that is the meaning of order-4 accuracy.
k i vectors
Each k i is a trial slope-list — the vector slope f evaluated at a guessed state. RK4 collects four of them (start, two middles, end) and takes a weighted average 6 1 ( k 1 + 2 k 2 + 2 k 3 + k 4 ) .
Intuition Why four peeks, weighted
1 , 2 , 2 , 1
Estimating the hop from the starting slope alone (that's the cheaper Euler's method for systems ) drifts because the slope changes during the hop. Peeking again in the middle and at the end, and trusting the middles more (weight 2 ), cancels far more error — matching a Taylor expansion out to h 4 . Because each k i is a list , you must fill the whole list before moving to the next peek.
vector rule f bold: coupling
reduction of higher order
exact integral over one hop
four probe slopes k1 to k4
Everything on the left is a single number or idea; the two boldings (y , f ) turn the scalar story into the system story that feeds RK4.
Read what a symbol means it maps to a picture; if not, reread that section.
y ( t ) the value of a moving quantity at time t — a dot sliding along a line.
y ′ the rate of change of y = the slope of its graph at an instant.
f ( t , y ) the slope-rule machine: input time and value, output the slope.
y ( t 0 ) = y 0 the starting point that selects which curve you're on (an IVP).
h the size of one time-hop; smaller = more accurate, more work.
y (bold)a stacked list of several quantities treated as one moving object.
f ( t , y ) the vector slope-rule; each component may depend on all others (coupling).
Coupling each equation's rate depends on the other unknowns, so you must advance all together.
Order of an ODE the highest prime-count; second-order reduces to a first-order system.
Reduction trick set u 1 = y , u 2 = y ′ so u 1 ′ = u 2 , u 2 ′ = g ( t , u 1 , u 2 ) .
k i a trial slope-list; RK4 averages four of them with weights 1 , 2 , 2 , 1 over 6 .
∫ t t + h f d s the exact accumulated change over one hop, estimated by sampling slopes.
O ( h 4 ) halving h shrinks global error about 16 × .
Parent topic — Hinglish
RK4 for a single ODE — the scalar method these symbols generalise.
Euler's method for systems — one-peek cousin, good for building intuition.
Reducing higher-order ODEs to first-order systems
Local vs Global Truncation Error — meaning of O ( h 4 ) .
Simpson's Rule — the sampling idea behind the 1 , 2 , 2 , 1 weights.
Stiff systems and stability — where step size h bites back.