TAILIEUCHUNG - Burden - Numerical Analysis 5e (PWS, 1993) Episode 2 Part 4

Tham khảo tài liệu 'burden - numerical analysis 5e (pws, 1993) episode 2 part 4', kỹ thuật - công nghệ, cơ khí - chế tạo máy phục vụ nhu cầu học tập, nghiên cứu và làm việc hiệu quả | 1 -f Ì 318 CHAPTER 5 B Initial-Value Problems for Ordinary Differential Equations Having computed tj tj and Wj we need to determine W 1 the solution to F y y - wj 2 To approximate this solution select usually as Wj and generate wJỈ Ị by applying Newton s method to ALGORITHM wji - Wj - f tj Wj f tj li wji 1 -----------------------2 -- --------------------- 1 - until is sufficiently small. Normally only three or four iterations per step are requừed. The Secant method could be used as an alternative to Newton s method in Eq. 5 67 but then two distinct initial approximations to w7 1 are required. To employ the Secant method the usual practice is to let Ị Wj and obtain from some explicit multistep method. When a system of stiff equations is involved a generalization is requừed for either Newton s or the Secant method. These topics will be considered in Chapter 10. Trapezoidal with Newton Iteration To approximate the solution of the initial-value problem ý ơ y y à a at N 1 equally spaced numbers in the interval a INPUT endpoints ứ -T integer N initial condition a tolerance TC L maximum number of iterations M at any one step. OUTPUT approximation w to y at the N 4- 1 values of t or a message of failure. Step 1 Set ỉ b à IN f a w a OUTPUT t w . Step 2 For i 1 2 do Steps 3-7. h Step 3 Set kỵ w ị f t w w0 7q j 1 FLAG 0. Step 4 While FLAG 0 do Steps 5-6. Stiff Differential Equations 319 j vvo ọ f t w o y - Stef 5 Set w W o------------------------. u. 1 - f fy t - h W0 . Step 6 If I w0 w01 TOL then set FLAG 1 cs else set j j 1 r w0 w . if j M then u OUTPUT MAXIMUM NUMBER OF È ITERATIONS EXCEEDED STOP. Step 7 Sett ar- iff 7 OUTPUT t w . Step 8 STOP. EXAMPLE 3 The stiff initial-value problem y 5e5t y - r 2 -r 1 0 i l y 0 -l has solution y f t e 5t. To show the effects of stiffness the Trapezoidal Algorithm and the Runge-Kutta fourth-order method are applied with N 4 h and N 5 h . The Trapezoidal method performs well in both cases using M . 10 and TOL

TỪ KHÓA LIÊN QUAN
TÀI LIỆU MỚI ĐĂNG
54    151    1    27-12-2024
26    137    2    27-12-2024
TAILIEUCHUNG - Chia sẻ tài liệu không giới hạn
Địa chỉ : 444 Hoang Hoa Tham, Hanoi, Viet Nam
Website : tailieuchung.com
Email : tailieuchung20@gmail.com
Tailieuchung.com là thư viện tài liệu trực tuyến, nơi chia sẽ trao đổi hàng triệu tài liệu như luận văn đồ án, sách, giáo trình, đề thi.
Chúng tôi không chịu trách nhiệm liên quan đến các vấn đề bản quyền nội dung tài liệu được thành viên tự nguyện đăng tải lên, nếu phát hiện thấy tài liệu xấu hoặc tài liệu có bản quyền xin hãy email cho chúng tôi.
Đã phát hiện trình chặn quảng cáo AdBlock
Trang web này phụ thuộc vào doanh thu từ số lần hiển thị quảng cáo để tồn tại. Vui lòng tắt trình chặn quảng cáo của bạn hoặc tạm dừng tính năng chặn quảng cáo cho trang web này.