Skip to main content

Angular velocity from rod data

Short answer

Assuming a corrected θ(0.6)=1.123, use Newton's forward formula with p=1. Differentiating yields angular velocity ω(0.2) ≈ 1.224 rad/s and angular acceleration α(0.2) ≈ 6.225 rad/s².

The following table gives the angle in radians (θ) through which a rotating rod has turned for various values of time in seconds (t). Find the angular velocity and angular acceleration at t = 0.2.

t00.20.40.60.8
θ00.1220.4930.1232.022

Data note. The printed θ at t = 0.6 is 0.123, which breaks the smooth trend (0, 0.122, 0.493, …, 2.022; the second differences are almost constant only if θ(0.6) = 1.123). It is taken as a misprint of 1.123.

The time steps are equal (h = 0.2) and t = 0.2 is the second tabulated point, so Newton's forward formula with t₀ = 0 and p = 1 is used. Angular velocity ω = dθ/dt and angular acceleration α = d²θ/dt².

Difference table

xyΔyΔ²yΔ³yΔ⁴y
00.0000.1220.2490.0100.000
0.20.1220.3710.2590.010
0.40.4930.6300.269
0.61.1230.899
0.82.022

Take t₀ = 0, h = 0.2, so p = (0.2 − 0)/0.2 = 1.00. Differentiating Newton's forward formula with respect to t gives

dθ/dt  = (1/h)[Δθ₀ + (2p−1)/2·Δ²θ₀ + (3p²−6p+2)/6·Δ³θ₀ + (4p³−18p²+22p−6)/24·Δ⁴θ₀]
d²θ/dt² = (1/h²)[Δ²θ₀ + (p−1)Δ³θ₀ + (6p²−18p+11)/12·Δ⁴θ₀]

First derivative

TermCoefficientValue of differenceProduct
Δθ₀1.00000.1220.1220
Δ²θ₀0.50000.2490.1245
Δ³θ₀-0.16670.010-0.0017
Δ⁴θ₀0.08330.0000.0000

Sum = 0.2448, so dθ/dt = 0.2448/0.2 = 1.2242

Second derivative

TermCoefficientValue of differenceProduct
Δ²θ₀1.00000.2490.2490
Δ³θ₀0.00000.0100.0000
Δ⁴θ₀-0.08330.0000.0000

Sum = 0.2490, so d²θ/dt² = 0.2490/0.2² = 6.2250

Answer: ω(0.2) ≈ 1.224 rad/s, α(0.2) ≈ 6.225 rad/s²

(With the printed 0.123 taken literally the same method gives ω ≈ 3.724 rad/s and α ≈ -2.108 rad/s², which is not physically consistent with the rest of the table.)

✓ Answer checked by Claude Opus 5.5 (high effort)

All Numerical Method chapter 5 questions →