Recursive sequence models
A recursive rule gives each term from the one before it, plus a start value. Two simple ones:
- Arithmetic: u(n+1) = u(n) + r. Each step adds r. Then u(n) = u(0) + n·r.
- Geometric: u(n+1) = q·u(n). Each step multiplies by q. Then u(n) = u(0)·qn.
Use them for equal yearly savings (arithmetic) or interest and population growth (geometric).
Arithmetico-geometric sequences
Many real models do both: u(n+1) = a·u(n) + b. If a ≠ 1, look for a fixed point L, a value that does not change: L = a·L + b, so L = b ÷ (1 − a).
Now write v(n) = u(n) − L. Then v(n+1) = a·v(n): geometric! So v(n) = an·v(0) and
u(n) = L + an·(u(0) − L)
If |a| < 1 the term an goes to 0, so u(n) tends to L. If a > 1 it grows without limit (unless u(0) = L).
The equation y' = ay + b
Now let time be continuous. The rate of change of y is y' = a·y + b. The solutions are
y(x) = C·eax − b/a (a ≠ 0), where C is found from the start value.
The constant solution y = −b/a is the equilibrium, like the fixed point. For a < 0 every solution tends to it; for a > 0 they move away.
Discrete versus continuous models
A discrete model changes at steps (once a year). A continuous model changes all the time. If u(n+1) = a·u(n) + b is matched with y' = k·y + c where a = ek, the curve passes exactly through all the points. Use discrete when you measure at fixed times (yearly census); use continuous when change is smooth (cooling, decay).
Logistic model
Plain growth can not go on forever. In a limited place the growth slows near a limit K. A simple logistic rule is
u(n+1) = u(n) + r·u(n)·(1 − u(n)/K)
When u is small, the bracket is near 1, so it grows almost geometrically. When u is near K, the bracket is near 0 and growth stops. The graph is an S-shaped curve. The continuous form is y' = r·y·(1 − y/K).
Euler's method
Most equations have no neat formula. Euler's method walks along the slope in small steps h:
y(n+1) = y(n) + h·f(y(n)) for y' = f(y).
Smaller h gives a better answer but needs more steps. The logistic bars in the 3D are exactly Euler steps with h = 1.
Try it: predict then check
In the 3D free play set a = 0.8, b = 4, start = 2. Predict the limit using b ÷ (1 − a). Check the green line. Now try a = 0.5, b = 4: predict, then check. Finally set a = 1. What kind of sequence is it?
Key formulas and definitions
- Arithmetic: u(n) = u(0) + n·r
- Geometric: u(n) = u(0)·q<sup>n</sup>
- Fixed point: L = b ÷ (1 − a), for a ≠ 1
- u(n) = L + a<sup>n</sup>·(u(0) − L)
- y' = ay + b gives y = C·e<sup>ax</sup> − b/a
- Euler: y(n+1) = y(n) + h·f(y(n))
- Logistic: u(n+1) = u(n) + r·u(n)·(1 − u(n)/K)
Worked examples
1. u(0) = 2 and u(n+1) = u(n) + 3. Find u(5).
Arithmetic with r = 3. u(5) = 2 + 5 × 3 = 17.
2. u(0) = 5 and u(n+1) = 2·u(n). Find u(4).
Geometric with q = 2. u(4) = 5 × 2⁴ = 5 × 16 = 80.
3. u(0) = 2 and u(n+1) = 0.8·u(n) + 4. Find the fixed point and u(3).
L = 4 ÷ (1 − 0.8) = 20. u(n) = 20 − 18 × 0.8ⁿ. u(3) = 20 − 18 × 0.512 = 20 − 9.216 = 10.784. Check by steps: u1 = 5.6, u2 = 8.48, u3 = 10.784.
4. You save Rs 1000, earn 5% a year and add Rs 100 each year. Write the rule and find the amount after 2 years.
u(n+1) = 1.05·u(n) + 100. u1 = 1050 + 100 = 1150. u2 = 1.05 × 1150 + 100 = 1207.5 + 100 = 1307.5. Rs 1307.50.
5. Solve y' = 2y + 4 with y(0) = 1, and find y(1).
y = C·e^(2x) − 4/2 = C·e^(2x) − 2. y(0) = 1 gives C − 2 = 1, so C = 3. y(1) = 3e² − 2 ≈ 3 × 7.389 − 2 ≈ 20.17.
6. Use Euler's method with h = 0.5 for y' = y, y(0) = 1, to estimate y(1). Compare with the exact value e ≈ 2.718.
y1 = 1 + 0.5 × 1 = 1.5. y2 = 1.5 + 0.5 × 1.5 = 2.25. The estimate is 2.25; the exact value is 2.718, so Euler underestimates by about 0.47. A smaller h gives a closer value.
Common mistakes
- Using L = b ÷ (1 − a) when a = 1. Then there is no fixed point (unless b = 0); the sequence is arithmetic.
- Thinking every u(n+1) = a·u(n) + b tends to a limit. It does only when |a| < 1.
- Forgetting the minus when solving y' = ay + b: the constant solution is −b/a, not b/a.
- Using a big h in Euler's method and trusting the result. Check with a smaller h.