Bài giảng Phương pháp số - Bài 2: Nghiệm của các phương trình phi tuyến

Nội dung text: Bài giảng Phương pháp số - Bài 2: Nghiệm của các phương trình phi tuyến

1. BÀI 2 NGHIỆM CỦA CÁC PHƢƠNG TRÌNH PHI TUYẾN
2. GIẢI PHƯƠNG TRÌNH f(x) = 0 KHÁI NIỆM CHUNG Bài toán Cho hàm f(x) liên tục trên đoạn [a, b] hoặc trên khoảng vô hạn và đƣờng cong y = f(x) chỉ có các nghiệm cô lập, tức là tồn tại các khoảng rời nhau chứa các không điểm của f(x) Các bƣớc giải 1- Tách nghiệm hay tìm khoảng cách li nghiệm (a, b) - chỉ chứa một nghiệm của phƣơng trình f(x) = 0 2- Kiện toàn nghiệm: tính gần đúng nghiệm với độ chính xác cho trƣớc Cơ sở của phƣơng pháp tách nghiệm Nếu hàm f(x) xác định và liên tục trên [a, b], f(a)f(b) < 0 và f’(x) giữ dấu trên (a, b) thì tồn tại duy nhất một nghiệm thực x* ∊ (a, b) của phƣơng trình f(x) = 0 PHƢƠNG PHÁP SỐ-Bài 2 2
3. GiẢI PHƯƠNG TRÌNH f(x) = 0 PHƢƠNG PHÁP TÁCH NGHIỆM Lập bảng xét dấu của đạo hàm cấp một f‘(x) rồi tìm các khoảng (a, b) thỏa mãn các điều kiện trên Ví dụ: Tìm các khoảng chứa các nghiệm cô lập của phƣơng trình f(x) = x3 – x – 1= 0 Giải: f‘ (x) = 3x2 – 1, lập bảng xét dấu sau x -∞ -2 -1 0 1 1.5 2 ∞ y’ + + + 0 _ 0 + + + y -7 5 0.875 -1 Vậy phƣơng trình trên có một nghiệm cô lập x1∊(1 ; 1.5) 3 PHƢƠNG PHÁP SỐ-Bài 2
4. CÁC PHƢƠNG PHÁP LẶP (1) Bắt đầu 1 PHƢƠNG PHÁP CHIA ĐÔI Giả thiết Cho f(x) liên tục trên (a, b) Nhập a, b, ε và f(a) f(b) < 0 m=(a+b)/2 y=f(x) đ f(a)f(m)<0 s b=m a=m a b x* (a+b)/2 s b-a<ε đ Kết thúc PHƢƠNG PHÁP SỐ-Bài 2 4
5. CÁC PHƢƠNG PHÁP LẶP (2) 1 PHƢƠNG PHÁP CHIA ĐÔI (tiếp) Thuật toán: Lặp với n = 0, 1, 2, , , cho đến khi tìm đƣợc nghiệm đúng x* hoặc nghiệm gần đúng xn đạt đƣợc độ chính xác mong muốn - Đặt m = (an + bn) / 2, nếu f(m) = 0 dừng (m là nghiệm đúng) - Nếu f(an) f(m) < 0, đặt an + 1 = an, bn + 1 = m - Trái lại, đặt an + 1 = m, bn + 1 = bn Vậy f(x) luôn luôn có không điểm trong khoảng [an + 1 ; bn + 1] . Sự hội tụ và sai số: Sử dụng phƣơng pháp chia đôi liên tiếp ta nhận đƣợc dãy khoảng lồng nhau {(an ; bn)} hữu hạn nếu x* là điểm giữa của khoảng thứ n, n hay vô hạn co lại: an < x* < bn f(an).f(bn) < 0, bn – an = (b – a) / 2 Khi n→ , do sự liên tục của f(x) nên lim bn= lim an = x* và an bn bn an b a x * xn xn 2 2 2n 1 PHƢƠNG PHÁP SỐ-Bài 2 5
6. CÁC PHƢƠNG PHÁP LẶP (3) 1 PHƢƠNG PHÁP CHIA ĐÔI (tiếp) Giải PT x3-x-1=0 trên đoạn [1; 1.5] với độ chính xác ε = 0.0005 n an bn c f(an) f( c) f(an) f( c) Sai số 0 1.0000 1.5000 1.2500 -1.0000 -0.2969 0.29688 0.2500 1 1.2500 1.5000 1.3750 -0.2969 0.2246 -0.06668 0.0625 2 1.2500 1.3750 1.3125 -0.2969 -0.0515 0.01529 0.0156 3 1.3125 1.3750 1.3438 -0.0515 0.0826 -0.00426 0.0039 4 1.3125 1.3438 1.3281 -0.0515 0.0146 -0.00075 0.0010 5 1.3125 1.3281 1.3203 -0.0515 -0.0187 0.00096 0.0002 9 1.3242 1.3252 1.3247 -0.0021 0.0000 0.00000 0.000001 6 PHƢƠNG PHÁP SỐ-Bài 2
7. CÁC PHƢƠNG PHÁP LẶP (4) 1. PHƢƠNG PHÁP CHIA ĐÔI (tiếp): CÀI ĐẶT double chiaDoi (double a, double b, double epsilon) { int lanlap = 0; // Khoi tao so lan lap double m ; do { lanlap++ ; m = (a + b) / 2.0 ; if (f(m) == 0) break; // m la nghiem dung else if (f(a)*f(m) > 0) a = m; else b = m; } while (! ( (b − a) / 2.0 1000) ); cout <<“So lanlap = " << lanlap << endl; return m; } PHƢƠNG PHÁP SỐ-Bài 2 7
8. CÁC PHƢƠNG PHÁP LẶP (5) 2 PHƢƠNG PHÁP THỬ SAI: Ý tưởng: Tăng tốc độ hội tụ của phƣơng pháp chia đôi bằng việc kiểm tra f(x) tại điểm trung bình có trọng số gần 0 hơn | f(b)|.a | f(a)|.b f(b).a f(a).b w | f(b)| | f(a)| f(b) f(a) Thuật toán: Lặp với n = 0, 1, 2, , , cho đến khi tìm đƣợc nghiệm đúng x* hoặc nghiệm gần đúng xn đạt đƣợc độ chính xác mong muốn Tính w = [f(bn)an – f(an)bn] / [f(bn) – f(an)], nếu f(w) = 0 dừng - Nếu f(an)f(w) < 0, đặt an + 1 = an, bn + 1 = w - Trái lại, đặt an + 1 = w, bn + 1 = bn 8 PHƢƠNG PHÁP SỐ-Bài 2
9. CÁC PHƢƠNG PHÁP LẶP (6) 2 PHƢƠNG PHÁP THỬ SAI (tiếp) f(b )a f(a )b w n n n n f(bn ) f(an ) w là điểm mà tại đó đƣờng thẳng cắt trục Ox và đi qua các điểm [an, f(an)] [bn, f(bn)] (một dây cung của f(x)) NHẬN XÉT: Phƣơng pháp này không ƣớc lƣợng đƣợc khoảng chứa nghiệm PHƢƠNG PHÁP SỐ-Bài 2 9
10. CÁC PHƢƠNG PHÁP LẶP (7) 2 PHƢƠNG PHÁP THỬ SAI (tiếp): Giải PT x3-x-1=0 trên đoạn [1; 1.5] n an bn f(an) f(bn) w f(w) f(an) f(w) 0 1.0000 1.5000 -1.0000 0.8750 1.2667 -0.2344 0.2344 1 1.2667 1.5000 -0.2344 0.8750 1.3160 -0.0370 0.0087 2 1.3160 1.5000 -0.0370 0.8750 1.3234 -0.0055 0.0002 3 1.3234 1.5000 -0.0055 0.8750 1.3245 -0.0008 0.0000 4 1.3245 1.5000 -0.0008 0.8750 1.3247 -0.0001 0.0000 PHƢƠNG PHÁP SỐ-Bài 2 10
11. CÁC PHƢƠNG PHÁP LẶP (8) 3 PHƢƠNG PHÁP DÂY CUNG: Cho hàm f(x) liên tục và 2 điểm x−1, x0. Coi xi+1 là nghiệm xấp xỉ thứ i+1 trên đoạn [xi -1 ; xi] ta nhận đƣợc công thức f(xi )xi-1 f(xi-1)xi xi 1 f(xi ) f(xi-1) xi xi 1 xi 1 xi f(xi ) f(xi ) f(xi 1) NHẬN XÉT: f(xi -1) và f(xi) không nhất thiết phải trái dấu, thậm chí nếu f(xi -1) = f(xi) thì không tính đƣợc xi+1 PHƢƠNG PHÁP SỐ-Bài 2 11
12. CÁC PHƢƠNG PHÁP LẶP (9) VÍ DỤ: Giải PT x3-x-1=0 trên đoạn [1; 1.5] với độ chính xác ε = 0.0005 bằng phƣơng pháp dây cung Điều n xn-1 xn xn - xn-1 f(xn-1) f(xn) f(xn) - f(xn-1) chỉnh Sai số -1 1.0000 1.5000 0.5000 -1.0000 0.8750 1.8750 0.2333 0.9375 0 1.5000 1.2667 -0.2333 0.8750 -0.2344 -1.1094 -0.0493 0.4375 1 1.2667 1.3160 0.0493 -0.2344 -0.0370 0.1973 -0.0093 0.0924 2 1.3160 1.3252 0.0093 -0.0370 0.0021 0.0392 0.0005 0.0173 3 1.3252 1.3247 -0.0005 0.0021 0.0000 -0.0021 0.0000 0.0009 4 1.3247 1.3247 0.0000 PHƢƠNG PHÁP SỐ-Bài 2 12
13. CÁC PHƢƠNG PHÁP LẶP (10) SAI SỐ CỦA CỦA PHƢƠNG PHÁP DÂY CUNG Giả sử 0 < m ≤ |f′(x)| ≤ M với mọi x [a, b], từ công thức f(xn )xn-1 f(xn-1)xn f(xn ) f(xn 1) x n 1 f(xn ) (x n 1 x n ) f(xn ) f(xn-1) x n x n 1 Theo công thức số gia hữu hạn, giả sử ξ là nghiệm đúng f(xn) f(xn 1) (ξ xn)f '(ξn) f(ξ) f(xn) (xn 1 xn) xn xn 1 f '(ηn)(xn 1 xn) hay (ξ xn 1 xn 1 xn)f '(ξn) f '(ηn)(xn 1 xn) | f '(ηn) | | f '(ξn)| M m |ξ xn 1| |x n 1 xn| |x n 1 xn| |f '(ξn)| m PHƢƠNG PHÁP SỐ-Bài 2 13
14. CÁC PHƢƠNG PHÁP LẶP (11) 4 PHƢƠNG PHÁP NEWTON Cho f(x) khả vi liên tục và một điểm x0, trong công thức dây cung thay độ nghiêng của dây cung bằng độ nghiêng của tiếp tuyến tại xn, ta nhận đƣợc công thức lặp f(x ) x x n n 1 n ' f (xn ) Chọn điểm xuất phát x ? 0 f(x0)f’’(x0) > 0 khi f’’ giữ dấu trên [a, b] PHƢƠNG PHÁP SỐ-Bài 2 14
15. CÁC PHƢƠNG PHÁP LẶP (12) VÍ DỤ: Giải PT x3–x–1=0 trên đoạn [1; 1.5] với độ chính xác ε = 0.0005 bằng phƣơng pháp Newton f(xn) / n x f(x ) f'(x ) f'(x ) Sai số Phương pháp này n n n n thường có dãy 0 1.0000 -1.0000 2.0000 -0.5000 nghiệm xấp xỉ nằm 1 1.5000 0.8750 5.7500 0.1522 0.5625 về một trong hai phía 2 1.3478 0.1007 4.4499 0.0226 0.0521 của nghiệm đúng, khi f’’ không đổi dấu 3 1.3252 0.0021 4.2685 0.0005 0.0012 4 1.3247 0.00002 Thử lại với x0 = 1.5 ?? PHƢƠNG PHÁP SỐ-Bài 2 15
16. CÁC PHƢƠNG PHÁP LẶP (13) SAI SỐ CỦA PHƢƠNG PHÁP NEWTON • Từ công thức lặp Newton và khai triển Taylor cấp 1 tại xn f ''(ξ ) f(x ) n | x x |2 (*) n 1 2 n 1 n Mặt khác, giả sử 0 < m1 ≤ |f′(x)| và |f′′(x)| ≤ M2 với mọi x є [a, b], ta có ' f(x n 1) f(xn 1) f(ξ ) f (ςn 1)(xn 1 ξ) | f(x ) | | x - ξ | n 1 n 1 m 1 M2 2 Thay (*) vào ta đƣợc | x n 1 - ξ | | x n 1 x n | 2m1 PHƢƠNG PHÁP SỐ-Bài 2 16
17. PHÉP LẶP ĐIỂM CỐ ĐỊNH (1) Khái niệm: Nếu phương trình f(x) = 0 (*) ↔ x = g(x) ( ) thì một nghiệm bất kì của ( ), tức là, bất kì điểm cố định nào của g(x), cũng đều là nghiệm của (*), g(x) được gọi là một hàm lặp để giải phương trình (*) Thuật toán 3.6: Phép lặp điểm cố định. Cho hàm lặp g(x) và một điểm xuất phát x0 Với n = 0, 1, 2, , cho đến khi thỏa mãn độ chính xác cho trước, lặp lại việc tính xn+1 = g(xn) PHƢƠNG PHÁP SỐ-Bài 2 17
18. PHÉP LẶP ĐIỂM CỐ ĐỊNH (2) Định lý: Cho g(x) thỏa mãn Hệ quả: Nếu g(x) là 1.Tồn tại khoảng [a, b] sao cho g(x) khả vi liên tục trong xác định trên [a, b] , g(x) є [a, b]. khoảng mở chứa điểm 2. g(x) khả vi trên [a, b], và tồn cố định ξ và |g’(ξ)| 0: mọi x є [a, b], , |g’(x)| ≤ K phép lặp điểm cố định Khi đó: xn+1 = g(xn) hội tụ bất 1. g(x) có đúng một điểm cố cứ khi nào |x0 – ξ| ≤ ε định ξ є [a, b], Điểm cố định ξ mà 2. Với điểm x0 є [a, b] bất kì, dãy |g’(ξ)| < 1, xn+1 = g(xn) luôn hội tụ về ξ. đƣợc gọi là điểm hút. Phƣơng pháp Newton là trƣờng hợp riêng của phép lặp điểm cố định PHƢƠNG PHÁP SỐ-Bài 2 18
19. PHÉP LẶP ĐIỂM CỐ ĐỊNH (3) Ví dụ: Tìm nghiệm dƣơng nhỏ n xn g(xn) nhất của phƣơng trình f(x) = x3 – x – 1 0 1.0000 1.2599 Giải 1 1.2599 1.3123 Bƣớc 1: . Chọn g(x) = (x + 1)1/3, với x є [1 ; 1.5], g(x) є [1; 1.5], 2 1.3123 1.3224 Bƣớc 2: g(x) khả vi và 3 1.3224 1.3243 g’(x)| = 1/[3(x + 1)2/3] < 1 trên [1; 1.5], phép lặp xn+1 = g(xn) sẽ 4 1.3243 1.3246 hội tụ nếu chọn điểm x0 є [1; 1.5], 1/3 xn+1 = (xn + 1) 5 1.3246 1.3247 PHƢƠNG PHÁP SỐ-Bài 2 19
20. PHÉP LẶP ĐIỂM CỐ ĐỊNH (4) SỰ HỘI TỤ VÀ PHÂN KÌ CỦA PHÉP LẶP ĐiỂM CỐ ĐỊNH g’>1 0<g’<1 g’<-1 -1<g’<0 PHƢƠNG PHÁP SỐ-Bài 2 20
21. PHÉP LẶP ĐIỂM CỐ ĐỊNH (5) • Sự hội tụ Chọn x0 là điểm bất kì trong [a, b], phép lặp điểm cố định tạo ra dãy x1, x2, . . . nằm trong [a, b]. Kí hiệu sai số ở lần lặp thứ n là en = ξ – xn n = 0, 1, 2, Từ ξ = g(ξ) và xn = g(xn–1), ta có en = ξ – xn = g(ξ) – g(xn–1) = g’(ηn) (ξ – xn) = g’(ηn) en–1 với ηn nằm giữa ξ và xn–1 (theo định lí về giá trị trung bình đối với đạo hàm) PHƢƠNG PHÁP SỐ-Bài 2 21
22. PHÉP LẶP ĐIỂM CỐ ĐỊNH (6) • Sự hội tụ Vì lim xn =  lim ηn =  nên limg'(ηn ) g'( ) n Vì g’(x) liên tục nên, en+1 = g’(ξ)en, ở đây lim en = 0. Do vậy, nếu g’(ξ) ≠ 0 thì với n đủ lớn 2 n | en+1 |≈ |g’(ξ)en| ≤ K|en| ≤ K |en–2| ≤ ≤ K |e0| 0 sai số en +1 ở lần lặp thứ (n + 1) phụ thuộc tuyến tính vào sai số en ở lần lặp thứ n. Vậy ta nói rằng x0, x1, x2, hội tụ tuyến tính về ξ PHƢƠNG PHÁP SỐ-Bài 2 22
23. PHÉP LẶP ĐIỂM CỐ ĐỊNH (7) • Sai số ξ – xn+1 = g’(ηn)(ξ – xn) = g’(ηn)(ξ – xn+1 + xn+1 – xn) g' (ηn )(xn 1 x n ) ξ x n 1 1 g' (ηn ) g' (ηn ) | ξ x n 1 | | x n 1 x n | 1 g' (ηn ) K | x x | 1 K n 1 n Kn | x x | 1 K 1 0 PHƢƠNG PHÁP SỐ-Bài 2 23
24. SỰ HỘI TỤ CỦA PHƢƠNG PHÁP NEWTON • Phƣơng pháp Newton là trƣờng hợp riêng của phép lặp điểm cố định, với g(x) = x – f(x)/f’(x) ta có g’(ξ) = 0. • Nếu g(x) khả vi hai lần, từ công thức Taylor en+1 = ξ – xn+1 = g(ξ) – g(xn) 2 = – g’(ξ)( xn – ξ) – g’’(ζn)(xn – ξ ) / 2 với ζn nào đó, nằm giữa ξ và xn, tức là 2 en+1 = g’(ξ)en – g” (ζn)en / 2 Vậy, nếu g’(ξ) = 0 và g”(x) liên tục tại ξ thì 2 en+1 ≈ – g” (ζn)en / 2 đối với n đủ lớn, tức là phƣơng pháp Newton có tốc độ hội tụ bậc 2. PHƢƠNG PHÁP SỐ-Bài 2 24
25. NHẬN XÉT CHUNG 1. Phƣơng pháp chia đôi hội tụ chậm nhƣng chắc chắn hội tụ 2. Phƣơng pháp Newton thƣờng hội tụ nhanh về nghiệm khi biết f’’(x) và chọn điểm xuất phát x0 gần nghiệm đúng, ít nhất khi f(x0)f’’(x0) > 0 (Định lí 3.2) 3. Phƣơng pháp dây cung thƣờng hiệu quả và đƣợc ƣa chuộng hơn. 4. Phép lặp điểm cố định chỉ hội tụ tuyến tính, nên không thực sự cạnh tranh với phƣơng pháp dây cung Bảng so sánh tốc độ hội tụ về nghiệm đúng của các phƣơng pháp lặp tìm nghiệm của phƣơng trình x3 – x – 1 = 0 với 5 chữ số có nghĩa Chia đôi Thử sai Dây cung Newton Lặp điểm cố định n 9 5 5 4 5 PHƢƠNG PHÁP SỐ-Bài 2 25
26. NGHIỆM CỦA ĐA THỨC (1) • ĐỊNH VỊ VÀ PHÂN LOẠI NGHIỆM CỦA ĐA THỨC 1. Định lí cơ bản của đại số Mọi đa thức bậc n với an ≠ 0 đều có đúng n nghiệm, kể cả thực và phức, nếu các nghiệm có bội r thì nó đƣợc tính r lần 2. Luật về dấu của Descartes Số np các nghiệm dƣơng của một đa thức p(x) là nhỏ hơn hay bằng số lần thay đổi về dấu υ của các hệ số p(x). Hơn nữa, hiệu υ – np là một số nguyên chẵn không âm xác định đƣợc số nghiệm thực của một đa thức với các hệ số thực PHƢƠNG PHÁP SỐ-Bài 2 26
27. NGHIỆM CỦA ĐA THỨC (2) • ĐỊNH VỊ VÀ PHÂN LOẠI NGHIỆM CỦA ĐA THỨC 3. Định lí về các biên của các nghiệm của đa thức Nếu p(x) là một đa thức với các hệ số ak nhƣ trong dạng lũy thừa, thì p(x) có ít nhất một nghiệm trong hình tròn |x| ≤ R xác định bởi R = min{p1, pn}, trong đó | a0 | | a0 | p1 n pn n | a1 | | a n | Ví dụ Nếu đa thức là p(x) = x5 – 3.7 x4 + 7.4 x3 – 10.8x2 + 10.8x – 6.8 thì a5 = 1, , a1 = 10.8, a0 = –6.8. | 6.8 | 6.8 p1 5 3.14 pn 5 1.46 R 1.46 |10.8 | 1 PHƢƠNG PHÁP SỐ-Bài 2 27
28. NGHIỆM CỦA ĐA THỨC (3) • ĐỊNH VỊ VÀ PHÂN LOẠI NGHIỆM CỦA ĐA THỨC 4. Định lí Cauchy: Mọi nghiệm của p(x) đều nằm trong miền hình khuyên r ≤ |x| ≤ R Với r và R là các nghiệm dƣơng nhỏ nhất và lớn nhất của 2 đa n n 1 thức P(x) | an | x | an 1 | x | a1 | x | a0 | n n 1 Q(x) | an | x | an 1 | x | a1 | x | a0 | Ví dụ Xét lại đa thức p(x) = x5 – 3.7 x4 + 7.4 x3 – 10.8x2 + 10.8x – 6.8 P(x) = x5 – 3.7 x4 – 7.4 x3 – 10.8x2 – 10.8x – 6.8 Q(x) = x5 + 3.7 x4 + 7.4 x3 + 10.8x2 + 10.8x – 6.8 Các nghiệm dƣơng của chúng là R = 5.6 , r = 0.63 tƣơng ứng. Vậy mọi nghiệm của p(x) đều phải thỏa mãn 0.63 < |x| ≤ 5.6 PHƢƠNG PHÁP SỐ-Bài 2 28
29. NGHIỆM CỦA ĐA THỨC (4) • BIỂU DIỄN ĐA THỨC Dạng lũy thừa 2 n p(x) = a0 + a1x+ a2x + anx Dạng Newton (tránh làm mất các chữ số có nghĩa khi tính) p(x) = a0 + a1 (x – c1) + a2 (x – c1) (x – c2) + a3 (x – c1) (x – c2)(x – c3) + + an(x – c1)(x – c2) (x – cn) Cho trƣớc giá trị z tính p(z): lặp lại các phép gán sau đây ' ' ' a n a n ; a i a i (z ci 1)a i 1, i n 1, n 2, , 0 ' a0 p(z) và p(x) p(z) (x - z)q(x), p'(z) q(z) " " " q(x) a1 a2 (x c1 ) a3 (x c1 )(x c2 ) " an (x c1 )(x c2 ) (x cn 1 ) PHƢƠNG PHÁP SỐ-Bài 2 29
30. NGHIỆM CỦA ĐA THỨC (5) PHƢƠNG PHÁP NEWTON TÌM NGHIỆM THỰC CỦA ĐATHỨC Cho n + 1 hệ số ao, . . . , an của đa thức p(x) dạng lũy thừa và nghiệm lặp ban đầu x0. Khi đó xm+1 = xm - p(xm)/q(xm) Với m = 0, 1, , cho đến khi thỏa mãn, lặp các công việc: z = xm ' " ' an an ,an an Với k = n – 1, , 1, lặp: ' ' ak ak zak 1 " ' " ak ak zak 1 ' ' a0 a0 za1 ' " xm 1 xm a0/a1 PHƢƠNG PHÁP SỐ-Bài 2 30
31. NGHIỆM CỦA ĐA THỨC (6) • VÍ DỤ Tìm các nghiệm của phƣơng trình đa thức p(x) = x3 + x – 3 = 0. Đa thức này có một nghiệm thực và hai nghiệm phức. Do p(1) = – 1 và p(2) = 7, nghiệm thực є (1, 2) chọn x0 = 1.1 PHƢƠNG PHÁP SỐ-Bài 2 31
32. NGHIỆM CỦA ĐA THỨC (7) • VÍ DỤ Để tìm các nghiệm phức, chúng ta áp dụng các công thức nghiệm của phƣơng trình bậc hai cho phƣơng trình đa thức 2 ' ' 2 x a2x a1 x 1.21341x 2.47236 0 Điều này dẫn đến các kết quả 2 a' (a' 4a' )1/2 x 2 2 1 2 -1.21341 2.90122i 0.60671 1.45061i 2 PHƢƠNG PHÁP SỐ-Bài 2 32
33. SỬ DỤNG GOAL SEEK TRONG EXCEL (1) PHƢƠNG PHÁP SỐ-Bài 2 33
34. BÀI TẬP 1. Bài tập tính toán: 3.1-3, 3.1-5, 3.1-8, 3.2-9, 3.3-2, 3.5-3, 3.5-4, 3.6-1 2. Bài tập lập trình: 3.2-6, 3.2-7, 3.3-3 PHƢƠNG PHÁP SỐ-Bài 2 35