Phương pháp ner tìm nghiệm gần đúng phương trình vi phân (LV01952)

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI NGUYỄN THỊ KIM THANH PHƯƠNG PHÁP NER TÌM NGHIỆM GẦN ĐÚNG PHƯƠNG TRÌNH VI PHÂN LUẬN VĂN THẠC SĨ TOÁN HỌC Hà Nội-2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM HÀ NỘI NGUYỄN THỊ KIM THANH PHƯƠNG PHÁP NER TÌM NGHIỆM GẦN ĐÚNG PHƯƠNG TRÌNH VI PHÂN LUẬN VĂN THẠC SĨ TOÁN HỌC Chuyên ngành: Toán giải tích Mã số: 60 46 01 02 Người hướng dẫn khoa học: TS Lê Đình Định Hà Nội-2016 Lời cảm ơn Trước trình bày nội dung luận văn, em xin bày tỏ lòng biết ơn sâu sắc tới TS Lê Đình Định, người tận tình hướng dẫn để em hoàn thành luận văn Em xin bày tỏ lòng biết ơn chân thành tới toàn thể thầy cô giáo phòng sau đại học, trường đại học sư phạm Hà Nội dạy bảo em tận tình suốt trình học tập Nhân dịp em xin gửi lời cảm ơn chân thành tới gia đình, bạn bè bên em, cổ vũ, động viên, giúp đỡ em suốt trình học tập hoàn thành luận văn Hà Nội, tháng 06 năm 2016 Học viên Nguyễn Thị Kim Thanh i Lời cam đoan Tôi xin cam đoan, bảo hướng dẫn TS Lê Đình Định, luận văn chuyên ngành Toán giải tích với đề tài: “Phương pháp NER tìm nghiệm gần phương trình vi phân” hoàn thành nhận thức tìm hiểu thân tác giả Trong trình nghiên cứu thực luận văn, tác giả kế thừa kết nhà khoa học với trân trọng biết ơn Hà Nội, tháng 06 năm 2016 Học viên Nguyễn Thị Kim Thanh ii Mục lục Mở đầu 1 Kiến thức chuẩn bị 1.1 1.2 Tổng quan phương trình vi phân 1.1.1 Phương trình vi phân 1.1.2 Bài toán Cauchy phương trình vi phân cấp 1.1.3 Bài toán Cauchy phương trình vi phân cấp n Một số phương pháp giải gần phương trình vi phân 1.2.1 Phương pháp Newton 1.2.2 Phương pháp Euler 17 1.2.3 Phương pháp Euler cải tiến 19 1.2.4 Phương pháp Richarson 21 Phương pháp phối hợp NER giải phương trình vi phân 23 2.1 Một số khái niệm định lý cần dùng 23 2.2 Phương pháp phối hợp NER 29 2.2.1 Biểu diễn nghiệm hệ phương trình vi phân 29 2.2.2 Lập nghiệm theo phương pháp Richarson 36 iii Kết luận 43 Tài liệu tham khảo 44 iv MỞ ĐẦU Lý chọn đề tài Phương trình vi phân phương trình toán học biểu diễn mối quan hệ hàm chưa biết với đạo hàm Phương trình vi phân đóng vai trò quan trọng nhiều lĩnh vực như: kĩ thuật, vật lý, kinh tế, Có nhiều phương pháp giải phương trình vi phân việc tìm nghiệm phương trình vi phân khó khăn Người ta tìm nghiệm vài phương trình vi phân đặc biệt đa số tìm nghiệm xấp xỉ Có hướng tìm nghiệm phương trình vi phân hướng tìm nghiệm giải tích hướng tìm nghiệm gần Các phương pháp giải tích phương pháp tìm nghiệm dạng biểu thức giải tích phương pháp hệ số bất định Các phương pháp giải số phương pháp Euler, phương pháp Runge-Kuta, tìm nghiệm dạng xấp xỉ Dưới hướng dẫn TS Lê Đình Định với mong muốn nghiên cứu tìm hiểu sâu sắc phương trình vi phân, chọn đề tài “Phương pháp NER tìm nghiệm gần phương trình vi phân” làm luận văn thạc sĩ Mục đích nghiên cứu • Nghiên cứu phương trình vi phân hệ thống lại số phương pháp giải phương trình vi phân • Phối hợp phương pháp giải phương trình vi phân thành phương pháp giải gần phương trình vi phân Nhiệm vụ nghiên cứu • Tìm hiểu phương trình vi phân thường số phương pháp giải phương trình vi phân • Phối hợp phương pháp biết đưa phương pháp giải gần phương trình vi phân Đối tượng phạm vi nghiên cứu • Đối tượng nghiên cứu: Phương trình vi phân • Phạm vi nghiên cứu: Nghiên cứu số phương pháp giải phương trình vi phân Phương pháp nghiên cứu • Phương pháp lý luận: Trước tiên đọc nghiên cứu tài liệu, giáo trình có liên quan đến phương trình vi phân số phương pháp giải phương trình vi phân Sau phân hóa, hệ thống kiến thức • Một số phương pháp giải phương trình vi phân Đóng góp luận văn Nội dung luận văn đưa phương pháp tìm nghiệm xấp xỉ phương trình vi phân Song bên cạnh phần kiến thức sở, luận văn đề cập tới số số vấn đề liên quan khác như: Các kiến thức phương trình vi phân, số phương pháp quen thuộc giải phương trình vi phân Những kiến thức đưa thêm vào không phục vụ cho việc làm luận văn mà giúp bạn đọc có thêm kiến thức để học tốt học phần Giải tích Chương Kiến thức chuẩn bị Trong chương này, trình bày số kiến thức giải tích số như: Tổng quan phương trình vi phân, toán Cauchy phương trình vi phân, số phương pháp quen thuộc giải phương trình vi phân ví dụ minh họa, Những kiến thức sử dụng để trình bày phương pháp NER tìm nghiệm gần phương trình vi phân Các kiến thức ta tìm thấy [1], [3] [4] 1.1 1.1.1 Tổng quan phương trình vi phân Phương trình vi phân Phương trình vi phân phương trình liên hệ biến độc lập hàm cần tìm đạo hàm Phương trình vi phân cấp n hệ thức có dạng F (x, y, y , y , , y (n) ) = 0, (1.1) thỏa mãn Khi u∗ nghiệm phương trình (2.25) u(x) nghiệm toán (2.26) ta có lim u(x) = u∗ x→∞ Ở ta áp dụng phương pháp Euler để giải toán (2.26) sau dùng phương pháp trung bình có trọng số Richarson để nâng cao độ xác nghiệm tìm Bằng phép đổi biến t= t x hay x = , ≤ t < 1, 1+x 1−t toán (2.26) trở thành  du −1   =− [F (u)] F (u), ≤ t < 1, dt (1 − t)   u(0) = u (2.27) (2.28) Khi ta có lim u(t) = u∗ , t→1−0 u(t) nghiệm toán (2.28) Bài toán (2.28) có tích phân t F (u) = F (u0 ).e− 1−t Do −t 1 1−t → 0, F (u) = F (u ).e (1 − t)2 (1 − t)2 t → − Vì ta áp dụng trực tiếp phương pháp Euler để giải toán (2.28) 30 Không tính tổng quát ta giải toán sau   du = f (t, u), ≤ t ≤ dt  u(0) = u0 , u0 ∈ X (2.29) với điều kiện f (t, u) ∈ C r {[0, 1] x X} , (2.30) u(t) ∈ C r+1 [0, 1] (2.31) Áp dụng phương pháp Euler bậc ta giải phương trình sai phân   uh (t) = f (t, uh ), t ∈ ωh t (2.32)  uh (0) = u , u ∈ X 0 M Các kết nâng cao độ xác nghiệm gần theo lưới chia ωh đoạn [0, 1] với h = phương pháp trung bình có trọng số Richarson phát triển cho phương pháp Euler bậc hai trường hợp X không gian Banach Giả sử biết nghiệm u(t) toán (2.29) với X không gian Banach ta xây dựng hệ k hàm vk (t) ∈ C r−k+1 [0, 1], k = 0; 1; ; r − phương pháp quy nạp Hệ hàm thỏa mãn hệ phương trình vi phân (2.15), với X không gian Banach Ta áp dụng phương pháp Euler bậc hai để giải gần toán (2.29) điều kiện (2.30) thỏa mãn 31 Giải hệ phương trình sai phân   wh (t) = f (t, wh ) + h utt (t), t ∈ ωh t  wh (0) = u , u ∈ X (2.33) với uh (t) tìm từ hệ phương trình (2.32) utt (t) định nghĩa Bổ đề 2.2 Giả sử biết nghiệm u(t) toán (2.29) với điều kiện (2.30) thỏa mãn, ta lập hệ hàm yk (t) ∈ C r−k+1 [0, 1], k = 0; 1; ; r − , thỏa mãn hệ phương trình vi phân    y0 (t) = u(t)     k+1 [ k+1 k ∂ sf ] v (2s) ds yk+1−s k+1−2s = s yi1 yis + , s  dt  s=1 s! s=1 s! ∂u i1 + +is =k s=1 (2s)!     y (0) = 0, k = 1; 2; ; r − k (2.34) Ta chứng minh hệ hàm tồn phương pháp quy nạp Với k = ta có: y0 (t) = u(t) Với k = phương trình (2.34) có dạng d2 y0 dy1 ∂f v0 (t) + = y + dt2 dt ∂u Khi biết y0 (t) = v0 (t) = u(t), phương trình phương trình vi phân tuyến tính cấp 1, có điều kiện ban đầu y1 (0) = Do hàm f (t, u) thỏa mãn điều kiện (2.30), hàm u(t) thỏa mãn điều kiện (2.31), nên phương trình tồn nghiệm y1 (t) ∈ C r [0, 1] Giả sử biết hàm yn (t) ∈ C r−n+1 [0, 1], n = 0, 1, , k − 32 Thay chúng vào phương trình vi phân (2.34) ta phương trình vi phân tuyến tính cấp yk , với điều kiện ban đầu yk (0) = Với độ trơn giả thiết (2.30) giả thiết quy nạp, chắn phương trình (2.34) tồn nghiệm yk (t) Trên sở hai hệ hàm vk (t) yk (t), k = 0; 1; ; r − ta chứng minh định lý sau Định lý 2.4 Nếu điều kiện (2.30) thỏa mãn, nghiệm toán sai phân biểu diễn dạng r−1 h hk yk (t) + hr µh (t), w (t) = u(t) + (2.35) k=1 với hàm yk (t) thỏa mãn hệ phương trình vi phân (2.34), k = 0; 1; ; r− Đồng thời nghiệm wh (t) uh (t) hệ phương trình sai phân tương ứng (2.32) hệ (2.32) tương ứng bị chặn theo h ωh wh (t) C,h ≤ c3 , (2.36) µh (t) bị chặn theo h ωh µh (t) C,h ≤ c4 , (2.37) số c3 c4 không phụ thuộc vào t h Chứng minh Thay (2.17) vào (2.35) (2.33) ta có r−1 r−1 k r h hk yk + hr µh ), h yk (t)t + h µ (t)t = f (t, k=0 k=0 + h r−1 k=0 h hk (vk )tt + ηtth 33 (2.38) Áp dụng Bồ đề 2.1 Bổ đề 2.2 vào (2.17) với m = r − k + 1, đồng thời r−1 khai triển Taylor hàm f (t, u) hk yk ta có: k=0 r−k+1 r−1 k h s=0 r−1 k=0 =f hs (s+1) yk (t) + hr (s + 1)! r−1 k h yk t, r h +h + k=0 r−1 +h r k=0 r−1 xr−k (t) + hr µht , k=0 r−1 k+1 s=0 k=0 d2s+2 yk+1−s h2s , (2s + 2)! dt2s+2 hr+1 h ρr−k (t) + η tt (2.39) r−k+1 với n = Trong σ h r−1 = 0(h), a = r−1 xr−k (t) b = k=0 ρr−k (t) không phụ k=0 thuộc vào h bị chặn theo t Thay đổi thứ tự lấy tổng ta r−1 r−k+1 k h s=0 k=0 hs (s+1) yk (t) = (s + 1)! ]−1 [ r−k+1 r−1 k+1 h s=0 k=0 r−1 k+1 h k s=1 k=0 (2s+2) (t) v = h2s k−s+1 (2s + 2)! ds yk−s+1 , s! dts [ k+1 ] r−1 k ,h s=1 k=1 (2s) vk+1 − 2s(t) (2s)! r−1 Đồng thời khai triển Taylor tới bậc r hàm f (2.40) t, hk y k (2.41) u = y0 (t) = k=0 u(t) ta có r−1 f r hk y k t, = f (t, u) + k=0 k=1 ∂ k f (t, u) k! ∂uk k r−1 hs yk + hr δ h , s=1 (2.42) δ h = 0(h), vì: r−1 r h hδ r−1 s =0 h ys r r hs−1 ys = h s=1 s=1 34 r Trong vế phải (2.42), ta khai triển lũy thừa bậc k thay đổi thứ tự lấy tổng ta thu kết quả: k r−1 k(r−1) s h ys hs = s=1 yi1 yik i1 + +ik =s s=k (r−1)2 r−1 f k h yk t, k k h = f (t, u) + k=0 s=1 k=1 ∂ s f (t, u) s! ∂us yi1 yis i1 + +is =k + hr δ h (2.43) Ta tách số hạng thứ hai vế trái (2.43) thành r2 k h ∂ s f (t, u) s! ∂us k s=1 k=1 r−1 k k = h s=1 k=1 r ∂ sf s! ∂us k r k−r +h h s=1 k=r yi1 yis i1 + +is =k yi1 yis i1 + +is =k ∂ sf s! ∂us yi1 yis i1 + +is =k Do f ∈ C r {[0, 1] x X} yk ∈ C r−k+1 [0, 1], nên r2 k k−r d= h k=r s=1 ∂ sf s! ∂us yi1 yis , i1 + +is =k bị chặn theo h Thay kết (2.40) (2.43) vào (2.39) ta thu r−1 k+1 k h s=1 k=0 ds yk+1−s + hr a + hr µht = f (t, u) + hr (dh + δ h s s! dt h + σ + b + ηtth ) + r−1 h [ k+1 ] + s=1 k k h k=1 s=1 ∂ s f (t, u) s! ∂us yi1 yis i1 + +is =k (2s) vk+1−2s (t) (2s)! (2.44) 35 Do hàm yk (t), k = 0; 1; ; r − 1, thỏa mãn hệ phương trình vi phân (2.34) nên h µht = b − a + dh + δ h + σ h + ηtt (2.45) Tới ta áp dụng Bổ đề 2.2 ta suy µh bị chặn theo h 2.2.2 Lập nghiệm theo phương pháp Richarson Trên sở định lý 2.4, ta lập nghiệm trung bình có trọng số theo phương pháp Richarson Ta lấy số cố định nguyên dương < N1 < N2 < < Nr Lập lưới cách ωhk với bước lưới tương ứng , k = 1; 2; ; r với M nguyên dương tùy ý Giải toán hk = M.Nk vi phân (2.32) (2.33) lưới ωhk ta nghiệm vi phân uhk whk tương ứng (k = 1; 2; ; r) Lập nghiệm trung bình có trọng số: r H vk whk , u = (2.46) k=1 , hệ số vk , k = 1; ; r, thỏa M mãn hệ phương trình đại số tuyến tính Trên lưới ωh với bước lưới h =        r vk = k=1 r k=1 vk hjk (2.47) = 0, j = 1; 2; ; r − hệ số vk tồn điều kiện N1 < N2 < < Nr Định lý 2.5 Nếu nghiệm uhk whk tìm từ hệ phương trình vi phân (2.32) (2.33) tương ứng điều kiện (2.30) thỏa mãn, 36 nghiệm uH tính theo công thức (2.46) với hệ số vk thỏa mãn hệ phương trình đại số tuyến tính (2.47) nghiệm gần toán (2.29) với sai số: uH − u(t) C,h ≤ C.hr , (2.48) C số không phụ thuộc vào h Chứng minh Thay (2.35) vào điều kiện u(0) = u0 ta có: r−1 = u(t) + Vì hk = k=1 r r vk hjk yj (t) j=1 vk hrk µhk + j=1 k=1 k=1 hjk yj (t) vk vk u(t) + u = r r−1 r r H vk hrk µhk + k=1 (2.49) k=1 , k = 1; ; r µhk bị chặn theo h, ta có: M.Nk vk ≤ + C5 C5 r , với C5 = Nk+1 − ;1 ≤ k ≤ r − Nk Vậy theo (2.47) ta có: r H vk hrk µhk u − u(t) = ≤ C4 k=1 uH − u(t) C,h + C5 C5 r r r h k=1 , Nk ≤ C.hr Trên sở định lý 2.1 định lý 2.2 ta xây dựng phương pháp giải gần phương trình (2.25) gồm bước sau đây: 37 Bước I: Đưa phương trình (2.25) toán (2.28) nhờ định lý 2.1 phép đổi biến (2.27) Bước II: Giải gần toán (2.28) hệ phương trình vi phân (2.32) (2.33) lưới chia ωhk Bước III: Lập nghiệm trung bình có trọng số theo công thức (2.46) với hệ số tìm từ hệ phương trình (2.47) Đây nghiệm gần toán (2.28) với đánh giá sai số (2.48) Nghiệm gần phương trình (2.25) nghiệm theo công M −1 thức (2.46) với t = M Như đưa thuật toán để giải gần phương trình toán tử phi tuyến (2.25) Thuật toán ta gọi phương pháp NER Thuật toán kết hợp phương pháp Newton – đưa thêm tham số liên tục với phương pháp Euler – giải gần phương trình vi phân phương pháp lập nghiệm trung bình có trọng số Richarson Phương pháp NER nhược điểm cồng kềnh đòi hỏi tính trơn cao toán tử phi tuyến F (u) Tuy phương pháp có ưu điểm tính ổn định đảm bảo, tốc độ hội tụ cao, thuật toán dễ đưa vào máy tính điện tử Ví dụ 2.1 Áp dụng phương pháp NER giải gần phương trình tích phân phi tuyến sau s[y(s)]2 ds y(x) = x − 16 + 12 38 (2.50) Giải Đặt s[y(s)]2 ds F (y) = y(x) − x + 16 − 12 Ta có F (y).u(x) = u(x) − 24 s.u(s)ds Để tìm dạng toán tử [F (y)]−1 ta giải phương trình tích phân sau u(x) − 24 s.u(s)ds = f (x) Đây phương trình tích phân tuyến tính với nhân suy biến, nên ta dễ dàng tìm kết thức Nghiệm phương trình cho công thức 24 u(x) = f (x) − 11 s.f (s)ds Do ta có 24 [F (y)] f (x) = f (x) − 11 −1 Khi toán (2.28) có dạng  dy   =   (1 − t)2   dt 24 x + − y(x, t) +  11 11      y(x, 0) = y0 (x), ≤ t < s.f (s)ds s.y(s)ds − 12 11 s.y (s)ds (2.51) Để tìm nghiệm gần với nghiệm xác biết y = x + phương trình (2.50), ta chọn nghiệm gần thỏa mãn điều kiện tương đối gần với nghiệm xác sau giải phương pháp Euler lập nghiệm theo công thức (2.46) 39 Hệ phương trình vi phân toán (2.51) có dạng  h   y(x, t + h) = y(x, t) +   (1 − t)   12 24 x+ s.y(s, t)ds− − y(x, t) + s.y (s, t)ds   11 11 11     y(x, 0) = y (x), t ∈ ω , h (2.52) ωh định nghĩa phần trước Ta giải hệ phương trình vi phân hai trường hợp: Nghiệm gần chọn y = x2 +1, M = 4, N1 = 1, N2 = 2, N3 = Nghiệm gần chọn y = x2 , M = 10, N1 = 1, N2 = 2, N3 = Trường hợp a Với h1 = 0, 25000 ta có y1 (x; 0, 7500) = 0, 00000x2 + 1, 00000x + 1, 01651 b Với h2 = 0, 125 ta có y2 (x; 0, 7500) = 0, 02151x2 + 0, 97849x + 1, 00528 c Với h3 = 0, 08333 ta có y3 (x; 0, 7500) = 0, 03291x2 + 0, 96709x + 1, 00549 Để lập nghiệm theo công thức (2.46), ta giải hệ phương trình đại số (2.47) với r = ν1 = 0, 5; ν2 = −4, 0; ν3 = 4, 40 Nghiệm gần phương trình tích phân (2.50) y(x) = y H (x; 0, 75000) = 0, 0621x2 + 0, 9379x + 1, 0118 Sai số theo lý thuyết nêu 0(h3 ) với h = 0, 25, sai số thực tế so với nghiệm xác 0, 012 Hai sai số hoàn toàn phù hợp với Trường hợp a Với h = 0, 10000 ta có y1 (x; 0, 90000) = 0, 01433x2 + 0, 98567x + 1, 01509 b Với h = 0, 05000 ta có y2 (x; 0, 90000) = 0, 00249x2 + 0, 99751x + 1, 00212 c Với h = 0, 03333 ta có y3 (x; 0, 90000) = 0, 00035x2 + 0, 99965x + 0, 99992 Nghiệm gần phương trình (2.50) theo công thức (2.46) y(x) = y H (x; 0, 90000) = −0, 0012x2 + 1, 0012x + 0, 9987 Sai số theo lý thuyết 0(h3 ) với h = 0, 1000, sai số thực tế 1, 4.10−3 Hai sai số hoàn toàn phù hợp với Qua ví dụ nhận thấy rằng: - Thuật toán trình bày thực hành thực tế - Thuật toán tưởng chừng cồng kềnh thực tế tính toán lại đơn giản dễ thực Trong tình tính toán dùng chung 41 công thức cho tất nghiệm trung gian hai trường hợp, công thức đơn giản dễ tính toán Chúng tính công cụ đơn giản lần tính đạt kết Như vậy, đưa vào máy tính điện tử dễ dàng đòi hỏi độ xác cao mà máy tính điện tử vận hành - Nghiệm gần chọn tốt tốc độ hội tụ nhanh Nhưng chọn chưa tốt trường hợp ta đạt kết theo yêu cầu đề Theo chúng tôi, trình tính toán đơn giản, phương pháp Newton cho phép chuyển sang dạng tuyến tính hóa được, phương pháp Euler đơn giản phương trình vi phân tuyến tính Những nhận xét phần minh họa ưu điểm thuật toán NER 42 Kết luận Luận văn trình bày phương pháp NER giải phương trình vi phân thường Các kết trình bày luận văn bao gồm: Các phương pháp quen thuộc giải phương trình vi phân thường là: Phương pháp Newton, Phương pháp Euler, Phương pháp Richarson, ưu, nhược điểm phương pháp ví dụ minh họa Phương pháp phối hợp NER giải phương trình vi phân thường ví dụ minh họa Mặc dù tác giả cố gắng, song kiến thức hạn chế nên luận văn không tránh khỏi thiếu sót, tác giả mong nhận ý kiến đóng góp Thầy, Cô giáo bạn 43 Tài liệu tham khảo [1] Phạm Kỳ Anh (2005), Giải tích số, Nhà xuất ĐHQG Hà Nội [2] Phan Đức Chính (1978), Giải tích hàm, Nhà xuất đại học trung học chuyên nghiệp Hà Nội [3] Nguyễn Minh Chương, Y.a.D.Mamedov, Khuất Văn Ninh (1992), Giải xấp xỉ phương trình toán tử, Nhà xuất khoa học kĩ thuật Hà Nội [4] Nguyễn Minh Chương, Nguyễn Văn Khải, Khuất Văn Ninh, Nguyễn Văn Tuấn, Nguyễn Tường (2000), Giải tích số, Nhà xuất Giáo dục [5] Phan Văn Hạp, Hoàng Đức Nguyên, Lê Đình Thịnh (1996), Phương pháp tính, Nhà xuất khoa học kĩ thuật [6] Hoàng Tụy (2003), Hàm thực giải tích hàm, Nhà xuất đại học quốc gia Hà Nội 44 [...]... cần tìm, y , y , , y (n) là các đạo hàm của hàm số y = y(x) Ta gọi cấp của phương trình vi phân là cấp cao nhất của đạo hàm có mặt trong phương trình Nghiệm của phương trình vi phân là một hàm số y = y(x) khi thay vào phương trình ta được một đồng nhất thức Khi n = 1 ta có phương trình vi phân cấp 1 1.1.2 Bài toán Cauchy đối với phương trình vi phân cấp 1 a Bài toán Cauchy Tìm nghiệm y = y(x) của phương. .. một nghiệm y = y(x) của phương trình (1.1) thỏa mãn các điều kiện ban đầu 6 1.2 Một số phương pháp giải gần đúng phương trình vi phân 1.2.1 Phương pháp Newton Trong mục này ta sẽ đi sâu nghiên cứu phương pháp Newton để tính nghiệm gần đúng của phương trình dạng tổng quát sau: f (x) = 0, (1.2) trong đó f (x) là một hàm số Hàm f (x) liên tục trên khoảng mà ta cần tìm nghiệm Phương trình (1.2) có nghiệm. .. dụng chương trình Newton để giải phương trình trong ví dụ trên, ta lại nhận được nghiệm x = 0.7 sau ba bước lặp, nhanh hơn khá nhiều so với phương pháp dây cung cùng với điều kiện ban đầu Ví dụ 1.2 Sử dụng phương pháp Newton để tìm nghiệm xấp xỉ của phương trình f (x) ≡ ex − cos(x) − 2 = 0 Giải Ta kết hợp chương trình AROOTB với chương trình Newton để giải phương trình trên, dùng một chương trình chính... [0,1] E + sup ta ∆t−a w(t + ∆t) − w(t) 0≤t