Đề tài Nghiên cứu một số đặc tính của enzyme β-Fructofuranosidase (β-ffase) từ nấm mốc aspergillus

TẠP CHÍ PHÁT TRIỂN KH&CN, TẬP 9, SỐ 6-2006 Trang 49 NGHIÊN CỨU MỘT SỐ ĐẶC TÍNH CỦA ENZYME β-FRUCTOFURANOSIDASE (β-Ffase) TỪ NẤM MỐC Aspergillus Mai Huỳnh Đoan Anh(1), Phạm Thị Ánh Hồng(1), Nguyễn Như Nhứt(2) (1)Trường Đại học Khoa học Tự Nhiên, ĐHQG - HCM (2) Công ty TNHH Gia Tường (Bài nhận ngày 08 tháng 11 năm 2005, hoàn chỉnh sửa chữa ngày 07 tháng 08 năm 2006) TÓM TẮT: Nhằm góp phần vào bước đầu xây dựng mô hình sản xuất enzyme Fructo- Oligosaccharide (FOS) sau này, chúng tôi tiến hành các giai đoạn thí nghiệm sau: 1. Khảo sát khả năng tổng hợp β-Ffase của một số chủng vi nấm Aspergillus được phân lập từ một số sản phẩm nông nghiệp ở Việt Nam và chọn chủng có khả năng tổng hợp cao nhất. 2. Khảo sát một số điều kiện để tổng hợp và thu nhận enzyme β-Ffase. 3. Khảo sát một số yếu tố ảnh hưởng lên hoạt tính thủy phân sucrose của enzyme β-Ffase như: pH, nhiệt độ, nồng độ cơ chất, động học thủy phân cơ chất. 1. MỞ ĐẦU β-Ffase (EC 3.2.1.26) là một enzyme thuộc họ GH 68 còn được gọi với nhiều tên khác như invertase; invertin; acid invertase; saccharase; glucosucrase; β-fructosidase; sucrase; β -h- fructosidase; fructosylinvertase. β-Ffase được tổng hợp bởi nhiều loài thực vật và vi sinh vật khác nhau với nhiều chức năng sinh học quan trọng như tham gia vào quá trình biến dưỡng và vận chuyển hexose, làm chậm quá trình lão ở thực vật. β-Ffase là một trong những enzyme có ứng dụng quan trọng trong sản xuất các fructo- oligosaccharide (FOS), một thành phần tối cần thiết và hữu ích trong các sản phẩm thực phẩm chức năng. FOS là một trong hơn 12 loại oligosaccharide đã được thương mại hóa. Ở Hàn Quốc, Bỉ, Pháp, Mỹ và nhất là Nhật Bản, FOS hấp dẫn người tiêu dùng trong những năm gần đây bởi FOS mang nhiều đặc tính chức năng có lợi cho sức khỏe con người như giảm cholesterol và mỡ trong máu, phòng ngừa bệnh tiểu đường và bệnh xơ cứng động mạch, kích thích hoạt động của hệ tiêu hóa, chống béo phì, không gây sâu răng….(1) Hiện nay, một lượng đáng kể FOS trên thế giới được sản xuất bằng enzyme β-Ffase thu nhận từ nấm mốc Aspergillus, đặc biệt là các chủng A. niger, A.japonicus, A.sydowi. Tuy nhiên, ở nước ta vẫn chưa có các nghiên cứu ứng dụng nào của β-Ffase vào việc sản xuất FOS. Mục tiêu của đề tài nhằm sàng lọc các chủng Aspergillus có khả năng tổng hợp β-Ffase cao và một số điều kiện để thu nhận cũng như vài đặc tính của enzyme này. 2. VẬT LIỆU - PHƯƠNG PHÁP THỰC NGHIỆM: 2.1. Vật liệu: Chúng tôi sử dụng 8 chủng Aspergillus phân lập từ các nguồn khác nhau. STT Tên chủng Ký hiệu Nguồn phân lập 1 Aspergillus niger A1 Phân lập từ thức ăn hỗn hợp gà cò 2 Aspergillus niger A2 Phân lập từ thức ăn hỗn hợp vịt cò 3 Aspergillus niger A3 Phân lập từ đậu phộng ở Củ Chi 4 Aspergillus niger A4 Phân lập từ thức ăn hỗn hợp heo Faco Science & Technology Development, Vol 9, No.6- 2006 Trang 50 5 Aspergillus sydowi A5 Phân lập từ tinh bột bắp 6 Aspergillus sydowi A6 Phân lập từ bột sò 7 Aspergillus foetidus A7 Phân lập từ tinh bột bắp 8 Aspergillus foetidus A8 Phân lập từ thức ăn hỗn hợp heo 2.2. Phương pháp 2.2.1. Phương pháp lên men: nuôi cấy lắc trong môi trường lỏng Hidaka (với 2% sucrose, 1,2% cao nấm men và 0,2% CMC; tất cả được pha trong đệm McIlvaine pH 5,0)(1). Thời gian nuôi là 4 ngày, tốc độ 200 vòng/phút ở 30oC với mật độ giống ban đầu 105 bào tử/100ml môi trường Hidaka. 2.2.2. Phương pháp xác định hoạt tính enzyme β-Ffase dựa trên nguyên tắc enzyme β- Ffase có khả năng cắt liên kết glucoside của sucrose tạo thành sản phẩm glucose và fructose. Do đó, dựa vào hàm lượng glucose sinh ra để xác định hoạt tính enzyme β-Ffase. Hỗn hợp phản ứng gồm 1 ml dung dịch đệm McIlvaine pH 5,0, 1 ml dung dịch sucrose 25% và 1 ml dịch β- Ffase. Phản ứng được thực hiện ở 40oC trong 10 phút và được dừng bằng cách đun sôi cách thủy trong 5 phút(4). Sau đó, xác định hàm lượng đường glucose trong hỗn hợp bằng phương pháp xác định đường khử theo Miler(7). Một đơn vị hoạt tính là lượng enzyme cần thiết để giải phóng 1 γ glucose trong một phút. 3. KẾT QUẢ VÀ BIỆN LUẬN 3.1. Khảo sát chọn chủng Aspergillus cho enzyme β-ffase có hoạt tính cao nhất Nuôi cấy lắc 8 chủng Aspergillus trong môi trường lỏng Hidaka với 105 bào tử/100ml môi trường ở nhiệt độ phòng 29oC đến 32oC, trong 4 ngày. Sau đó thu dịch enzyme ngoại bào thô và tiến hành xác định hoạt tính của dịch enzyme β-Ffase ngoại bào. Dựa vào kết quả hoạt tính của enzyme để chọn chủng Aspergillus cho hoạt tính enzyme cao nhất. Bảng 1.Khả năng tổng hợp β-Ffase của 8 chủng nấm mốc trên môi trường Hidaka Chủng Aspergillus ODTB thử không ODTB thử thật ∆OD575 Nồng độ Glucose (γ/ml) Hoạt tính (UI/mlDD.E) A1 0,406 0,852 0,466 310,7 0,518 A2 0,454 0,636 0,182 121,3 0,202 A3 0,466 0,681 0,215 143,3 0,239 A4 0,483 0,756 0,274 182,7 0,304 A5 0,079 0,169 0,090 60,0 0,100 A6 0,168 0,437 0,269 179,3 0,299 A7 0,082 0,257 0,175 116,7 0,194 A8 0,075 0,369 0,294 196,0 0,327 TẠP CHÍ PHÁT TRIỂN KH&CN, TẬP 9, SỐ 6-2006 Trang 51 0 0,1 0,2 0,3 0,4 0,5 0,6 A1 A2 A3 A4 A5 A6 A7 A8 Chuûng Aspergillus H oa ït t ín h (U I/m ld dE ) Hình 1. Khả năng tổng hợp β-Ffase của 8 chủng nấm mốc trên môi trường Hidaka Kết quả trên cho thấy khả năng tổng hợp β-Ffase của các chủng nấm mốc rất khác nhau. Trong số 8 chủng được khảo sát thì chủng A1 có khả năng tổng hợp β-Ffase cao nhất trên môi trường Hidaka nên được chọn tiếp tục cho các nghiên cứu tiếp theo. 3.2. Tổng hợp và thu nhận β-Ffase 3.2.1. Ảnh hưởng của pH môi trường lên sự tổng hợp β-Ffase của chủng chọn được Nuôi cấy chủng Aspergillus niger A1 trong môi trường Hidaka có pH thay đổi từ pH2,2- pH8,0. Kết quả kiểm tra hoạt tính được trình bày trong bảng 2. Bảng 2. Ảnh hưởng của pH môi trường lên khả năng tổng hợp β-Ffase của chủng A1 pH ODTB thử không ODTB thử thật ∆OD575 Nồng độ Glucose (γ/ml) Hoạt tính (UI/mlDD.E) 2,2 0,008 0,012 0,004 2,67 0,004 3,0 0,01 0,045 0,035 23,33 0,038 4,0 0,244 0,574 0,334 222,67 0,371 5,0 0,301 0,790 0,489 326,00 0,543 5,6 0,123 0,846 0,723 482,00 0,803 6,0 0,105 0,785 0,680 453,33 0,756 6,5 0,115 0,752 0,637 424,67 0,707 7,0 0,098 0,506 0,408 272,00 0,453 7,5 0,07 0,237 0,167 111,33 0,185 8,0 0,06 0,106 0,046 30,67 0,051 Science & Technology Development, Vol 9, No.6- 2006 Trang 52 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0 0,5 1 1,5 2 2,5 3 3,5 4 4,5 5 5,5 6 6,5 7 7,5 8 8,5 pH H oa ït t ín h (U I/ m ld dE ) Hình 2. Ảnh hưởng của pH môi trường lên khả năng tổng hợp β-Ffase của chủng A1 Các kết quả đạt được cho thấy chủng A1 có khả năng tổng hợp β-Ffase trong một dãy pH từ 5,0 đến 7,0; tuy nhiên, pH 5,6 là pH tối ưu nhất. 3.2.2. Thời gian thu nhận β-Ffase Nuôi cấy lắc chủng A1 Hidaka có pHop = 5,6 sao cho 105 bào tử/100ml. Sau mỗi 6 giờ thu dịch enzyme và kiểm tra hoạt tính một lần. Bảng 3. Sự biến thiên hoạt tính β-Ffase theo thời gian nuôi cấy Thời gian (giờ) Độ pha loãng (lần) ODTB thử không ODTB thử thật ∆OD575 Nồng độ Glucose (γ/ml) Hoạt tính (UI/mlddE) 0 50 0 0 0 0 0 6 50 0,010 0,014 0,004 2,784 0,232 12 50 0,013 0,036 0,023 15,68 1,307 18 50 0,013 0,071 0,058 38,92 3,246 24 50 0,014 0,100 0,086 57,41 4,784 30 50 0,013 0,115 0,102 68,09 5,674 36 50 0,012 0,128 0,116 77,65 6,471 42 50 0,015 0,153 0,138 91,89 7,658 48 50 0,014 0,159 0,145 96,67 8,056 54 50 0,013 0,157 0,144 96,12 8,010 60 50 0,014 0,152 0,138 91,85 7,654 66 50 0,015 0,147 0,132 87,94 7,328 72 50 0,013 0,135 0,122 81,74 6,812 TẠP CHÍ PHÁT TRIỂN KH&CN, TẬP 9, SỐ 6-2006 Trang 53 0 1 2 3 4 5 6 7 8 9 0 6 12 18 24 30 36 42 48 54 60 66 72 78 Thôøi gian (giôø) H oa ït t ín h (U I/ m ld dE ) Hình 3. Sự biến thiên hoạt tính β-Ffase theo thời gian nuôi cấy Từ 6 giờ sau khi nuôi cấy, hoạt tính β-Ffase tăng liên tục và tăng chậm ở khoảng 42 đến 48 giờ, sau đó, hoạt tính β-Ffase giảm dần. Do đó, có thể thu nhận enzyme sau 48 giờ nuôi cấy. Thời gian ổn định lâu (từ 48 đến 54 giờ) là ưu điểm cho quá trình thu nhận tiếp theo. 3.2.3. Ảnh hưởng của tỷ lệ giống ban đầu lên sự tổng hợp β-Ffase Nuôi cấy chủng A1 theo các điều kiện tối ưu với tỷ lệ giống thay đổi từ 101 đến 107 bào tử/100ml môi trường. Kết quả thu được như sau: Bảng 4. Ảnh hưởng của tỷ lệ giống Tỷ lệ giống 10n bào tử/100mlMT Độpha loãng (lần) ODTB thử không ODTB thử thật ∆OD575 Nồng độ Glucose (γ/ml) Hoạt tính (UI/mldd. E) 101 50 0,017 0,053 0,036 24,19 2,016 102 50 0,016 0,076 0,060 40.09 3,341 103 50 0,017 0,114 0,097 65,07 5,423 104 50 0,018 0,142 0,124 82,69 6,891 105 50 0,018 0,163 0,145 96.72 8,060 106 50 0,016 0,128 0,112 74,59 6,216 107 50 0,018 0,072 0,054 36,19 3,016 0 2 4 6 8 10 101 102 103 104 105 106 107 Baøo töû/100ml moâi tröôøng H oa ït t ín h (U I/ m ld d. E ) Hình 4. Ảnh hưởng của tỷ lệ giống Science & Technology Development, Vol 9, No.6- 2006 Trang 54 Thông qua các số liệu thí nghiệm ở bảng 4 và đồ thị 3 cho thấy khi tỷ lệ giống bổ sung vào môi trường lên men càng tăng thì enzyme β-Ffase thu được có hoạt tính càng cao. Nhưng chỉ đạt cực đại ở tỷ lệ 105 bào tử/100ml môi trường (với hoạt tính là 8,060 UI/mldd.E). Khi tỷ lệ giống vượt quá ngưỡng 105 trong điều kiện thành phần dinh dưỡng trong môi trường không được bổ sung nên kéo theo hoạt tính enzyme không tăng. 3.2.4. Tinh sạch sơ bộ β-Ffase Trước khi tiến hành xác định các yếu tố ảnh hưởng lên hoạt tính xúc tác, β-Ffase được tinh sạch sơ bộ bằng ethanol, acetone và sulfate ammonium. Kết quả thí nghiệm được trình bày trong bảng 5. Bảng 5. Ảnh hưởng của các tác nhân tủa lên khả năng tinh sạch β-Ffase Chế phẩm enzyme Tỷ lệ tác nhân tủa/dd. β-Ffase Hoạt tính chung (UI/gCP.E) Hoạt tính riêng (UI/mg protein) CP.E tủa bằng ethanol 1 : 2 736,04 49,61 CP.E tủa bằng aceton 1 : 1,5 724,77 48,43 CP.E tủa bằng (NH4)2SO4 60% 418,78 36,61 Kết quả trên cho thấy enzyme β-Ffase do chủng A1 tổng hợp thích hợp để tinh sạch bằng các hai dung môi hữu cơ ethanol và aceton hơn là muối sulfate amonium. Chế phẩm được tủa bằng ethanol và aceton co hoạt tính gần bằng nhau, cao nhất là chế phẩm được tủa bằng ethanol (49,61 UI/mg protein). 3.3. Một số yếu tố ảnh hưởng lên hoạt tính thủy phân sucrose của β-Ffase 3.3.1. Ảnh hưởng của pH Bảng 6. Ảnh hưởng pH lên phản ứng thủy phân của enzyme β-Ffase pH ODTB thử Không ODTB thử thật ∆OD575 Nồng độ Glucose (γ/ml) Độ pha loãng (lần) Hoạt tính (UI/gCP.E) 4.0 0,005 0,010 0,005 3,33 50 34,68 4.6 0,009 0,027 0,018 12,01 50 125,10 5.0 0,015 0,076 0,061 40,66 50 423,54 5.6 0,015 0,129 0,114 76,00 50 791,67 6.0 0,016 0,125 0,109 72,67 50 756,98 6.6 0,017 0,105 0,088 58,67 50 611,14 7.0 0,016 0,082 0,066 44,00 50 405,88 TẠP CHÍ PHÁT TRIỂN KH&CN, TẬP 9, SỐ 6-2006 Trang 55 0 100 200 300 400 500 600 700 800 900 3 3,5 4 4,5 5 5,5 6 6,5 7 7,5 pH H oa ït t ín h (U I/g C P. E ) Hình 5. Ảnh hưởng pH lên phản ứng thủy phân của enzyme β-Ffase Enzyme β-Ffase do chủng A1 tổng hợp có dãy pH hoạt động trong vùng acid và tối ưu nhất ở pH 5,6. Kết quả này phù hợp với kết quả của Sarote Sirisansaneeyakull và cộng tác viên (2000) khi nghiên cứu β-Ffase của Aspergillus niger ATCC 20611 (pHop 5,5) và các nghiên cứu trước đây trên enzyme β-Ffase do nấm men và thực vật tổng hợp. 3.3.2. Ảnh hưởng của nhiệt độ Bảng 7. Ảnh hưởng nhiệt độ lên phản ứng thủy phân của enzyme β-Ffase Nhiệt độ (oC) ODTB thử không ODTB thử thật ∆OD575 Nồng độ Glucose (γ/ml) Độ pha loãng (lần) Hoạt tính (UI/gCP.E) 30 0,018 0,074 0,056 37,33 50 388,85 35 0,018 0,082 0,064 42,67 50 444,48 40 0,018 0,131 0,113 75,33 50 784,68 45 0,018 0,153 0,135 90,00 50 937,50 50 0,018 0,181 0,163 108,67 50 1131,98 55 0,018 0,173 0,155 103,33 50 1076,35 60 0,018 0,165 0,147 98,00 50 1020,83 65 0,018 0,148 0,130 86,67 50 902,81 70 0,018 0,126 0,108 72,00 50 750,00 75 0,018 0,108 0,090 60,00 50 625,00 80 0,018 0,063 0,045 30,00 50 312,50 Science & Technology Development, Vol 9, No.6- 2006 Trang 56 0 200 400 600 800 1000 1200 25 30 35 40 45 50 55 60 65 70 75 80 85 Nhieät ñoä (oC) H oa ït t ín h (U I/g C P. E ) Hình 6. Ảnh hưởng nhiệt độ lên phản ứng thủy phân của enzyme β-Ffase Kết quả bảng 7 và hình 6 cho thấy enzyme này có dãy nhiệt độ hoạt động rộng (từ 30 đến 75oC); nhiệt độ tốt nhất cho hoạt tính thủy phân sucrose là 50oC với hoạt tính là 1131,98 UI/g chế phẩm. So với β-Ffase do Aspergillus niger ATCC 20611(1, 2) tổng hợp (top 40oC) thì enzyme này có mức nhiệt độ tối ưu cao hơn. 3.3.3. Ảnh hưởng của nồng cơ chất Bảng 8. Ảnh hưởng nồng độ sucrose lên phản ứng thủy phân của enzyme β-Ffase Nồng độ Sucrose (%) ODTB Thử không ODTB Thử thật ∆OD575 Nồng độ Glucose (γ/ml) Độ pha loãng (lần) Hoạt tính (UI/gCP.E) 5 0,014 0,085 0,071 47,33 50 493,02 10 0,014 0,096 0,082 54,67 50 569,48 15 0,016 0,153 0,137 91,33 50 951,35 20 0,015 0,182 0,167 111,33 50 1159,69 25 0,017 0,178 0,161 107,33 50 1118,02 30 0,016 0,145 0,129 86,00 50 895,83 0 200 400 600 800 1000 1200 1400 0 5 10 15 20 25 30 35 40 Noàng ñoä sucrose (%) H oa ït t ín h (U I/ g C P. E ) Hình 7. Ảnh hưởng nồng độ cơ chất lên phản ứng thủy phân của enzyme β-Ffase TẠP CHÍ PHÁT TRIỂN KH&CN, TẬP 9, SỐ 6-2006 Trang 57 Các số liệu cho thấy hoạt tính β-Ffase cao nhất tại nồng độ sucrose là 20% (1159.69 UI/gCP.E). Khi nồng độ sucrose cao hơn 20% thì hoạt tính β-Ffase giảm dần vì nồng độ sucrose càng cao làm hỗn hợp phản ứng sánh đặc hơn, đồng thời, có thể cùng lúc 2 phân tử sucrose gắn vào trung tâm hoạt động của β-Ffase làm enzyme này bị bất hoạt (Arnd Sturm, 1999). 3.3.4. Động học thủy phân Bảng 9. Động học thủy phân của β-Ffase Thời gian (phút) ODTB thử không ODTB thử thật ∆OD575 Nồng độ Glucose (γ/ml) 5 0,019 0,185 0,166 110,67 10 0,019 0,190 0,171 113,67 20 0,019 0,204 0,185 123,33 30 0,019 0,213 0,194 129,87 40 0,019 0,215 0,196 131,12 50 0,019 0,216 0,197 131,49 60 0,019 0,217 0,198 131,58 70 0,019 0,217 0,198 131,58 105 110 115 120 125 130 135 0 10 20 30 40 50 60 70 80 Thôøi gian (phuùt) N oàn g ño ä G lu co se (g am a/ m l) Hình 8. Động học thủy phân sucrose của enzyme β-Ffase Kết quả cho thấy enzyme β-Ffase do chủng A1 hoạt động tốt nhất kể từ thời điểm 30 phút sau khi tiếp xúc với cơ chất. 4. KẾT LUẬN Thông qua các kết quả thực hiện thí nghiệm, chúng tôi rút ra được một số kết luận sau: • So sánh khả năng tổng hợp β-Ffase trong 8 chủng Aspergillus dùng trong thí nghiệm thì chủng A1 (Aspergillus niger) là chủng có khả năng tổng hợp enzyme β-Ffase cao nhất trên môi trường Hidaka. • Điều kiện thích hợp để tổng hợp và thu nhận β-Ffase từ chủng A1 như sau: - pHop: 5,6. - Tỷ lệ giống: 105 bào tử/100ml môi trường Hidaka. - Thời gian lên men: 48 giờ. - Tinh sạch bằng ethanol với tỷ lệ 1 dịch chiết enzyme : 1 ethanol. • Điều kiện tối ưu để enzyme β-Ffase hoạt động thủy giải sucrose: Science & Technology Development, Vol 9, No.6- 2006 Trang 58 - pHop:5.6 - Nhiệt độ: 50oC - Nồng độ sucrose: 20% - Thời gian phản ứng: 30 phút. RESEARCH SOME CHARACTERISTICS OF ENZYME β-FRUCTOFURANOSIDASE (β-Ffase) FROM Aspergillus Mai Huynh Đoan Anh(1), Pham Thi Anh Hong(1), Nguyen Nhu Nhut (2) (1) University of Natural Sciences, VNU- HCM (2) Gia Tuong Company ABSTRACT: To construct the manufacturing process of fructo-oligosacharide enzime (FOS) in the future, we carried out the following reseach stages : 1. Studying some types of Aspergillus to synthesize enzyme β - Ffase. We seperated the stypes from agricultural products in Viet Nam and chose the types which have the ability to synthesize enzyme β - Ffase with high productivity. 2. Studying some conditions for culturing Aspergillus to have high producttivity of enzyme β - Ffase. 3. Studying some factors influencing hydrolytic ability enzyme β - Ffase to produce sucrose such as : pH, temperature, concentration of substrates, hydrolytic kinetics of substrates. TÀI LIỆU THAM KHẢO [1]. Carol Brannon, Prebiotic: Feeding friendly bacteria, Today’s Dietitian., 2003. [2]. Mirza Ahsen Baig, Kiran Shafiq, Sikander Ali, and Shazia Mirza, Optimization of Cultural Condition for the Biosynthesis of β- fructofuranosidase by Saccharomyces cerevisiae, Pakistan Journal of Biological Science, Vol.6, p. 952-954. ,2003. [3]. Mirza Ahsen Baig, Kiran Shafiq, Sikander Ali, and Shazia Mirza, Effect of Nitrogen and Phosphate Sources on the Biosynthesis of β- fructofuranosidase, Online Journal of Biological Science, Vol. 3, p. 519-595., 2003. [4]. Sarote Sirisansaneeyakul1, Sanya Jitbanjongkit1, Nararit Prasomsart1 and Pairojana Luangpituksa2, Production of β-Fructofuranosidase from Aspergillus niger, ATCC 20611. Kasetsart J. (Nat. Sci.), Vol. 34, p. 378–386., 2000. [5]. Sarote Sirisansaneeyakul, Sittiwat Lertsiri., Enzymatic production of fructo- oligosaccharides from sucrose, Kasetsart J (Nat. Sci.), Vol. 34, p. 262-269., 2000. [6]. Tirso Pons, Daniil G. Naumoff, Carlos Martonez-Fleites,1 and Lazaro Hernandez., Three Acidic Residues Are at the Active Site of a Propeller Architecture in Glycoside Hydrolase Families 32, 43, 62, and 68, Proteins: Structure, Funtion, and Bioinformatics, Vol. 54, p. 424-432., 2004. [7]. httpp://www.glue.umd.edu/~nsw/ench485/lab9d.htm#Method