Difference between revisions of "Length of a curve using numerical methods"
Line 24: | Line 24: | ||
and the answer is thus = (1/6)(g(0) + 4g(1/2) + 2g(1) + 4g(3/2) + 2g(2) + 4g(5/2) + 2g(3) + 4g(7/2) + 2g(4) + 4g(9/2) + g(5)) | and the answer is thus = (1/6)(g(0) + 4g(1/2) + 2g(1) + 4g(3/2) + 2g(2) + 4g(5/2) + 2g(3) + 4g(7/2) + 2g(4) + 4g(9/2) + g(5)) | ||
+ | |||
+ | = (1/6)(1 + |
Revision as of 19:04, 1 April 2021
Simpson's Rule is used to calculate integrals numerically.
From the Length of curve page, the length can be calculated using the integral: ∫sqr(1+f '(x)2)dx
---
Instead of evaluating this integral with a table, a numerical method called Simpson's rule can be used:
the length of the curve after n iterations = Ln = (Δx/3)(f(x0) + 4f(x1) + 2f(x2) + 4f(x3) + 2f(x4) + ...+ 2f(xn-2) + 4f(xn-1) + f(x0))
where Δx = (b-a)/n with a and b being the boundaries and n the number of iterations of the calculation.
---
The length of curve equation used in the example was ∫sqr(1+(30e−0.3t)2)dt
using boundary conditions of 0 and 5 minutes, and using 10 iterations, the problem becomes:
g(t) = sqr(1+(30e−0.3t)2
Δt = (5-0)/10 = 1/2
the subintervals are {0, 1/2, 1, 3/2, 2, 5/2, 3, 7/2, 4, 9/2, 5}
and the answer is thus = (1/6)(g(0) + 4g(1/2) + 2g(1) + 4g(3/2) + 2g(2) + 4g(5/2) + 2g(3) + 4g(7/2) + 2g(4) + 4g(9/2) + g(5))
= (1/6)(1 +