Skip to main content

Romberg: e^(-x²) from 1 to 2

Short answer

Build trapezoidal values for h = 1, 0.5, 0.25, 0.125 and apply Richardson extrapolation. The successive results agree to three decimals, giving the integral approximately 0.135.

Evaluate ∫ from 1 to 2 of e^(-x²) dx using Romberg method correct up to 3 decimal places.

Romberg integration: build trapezoidal values T(h) for h = 1, 0.5, 0.25, 0.125, then extrapolate with R(i,j) = [4ʲ·R(i,j−1) − R(i−1,j−1)] / (4ʲ − 1), until two successive values agree to 3 decimals.

f(x) = e^(−x²), a = 1, b = 2.

Romberg table (6 decimals kept for accuracy)

hT (order h²)R1R2R3
10.193098
0.50.1492480.134632
0.250.1387200.1352100.135249
0.1250.1361210.1352540.1352570.135257

Check: R2 at h = 0.25 is 0.135249 and R3 at h = 0.125 is 0.135257; both round to 0.135, so the required 3-decimal accuracy is reached.

Answer: ∫₁² e^(−x²) dx ≈ 0.135 (exact value 0.135257).

✓ Answer checked by Claude Opus 5.5 (high effort)

All Numerical Method chapter 5 questions →