diff --git a/DA_THPPT_5.pdf b/DA_THPPT_5.pdf new file mode 100644 index 0000000..ad2a5b3 Binary files /dev/null and b/DA_THPPT_5.pdf differ diff --git a/FINAL/DoAnMatlab.mlapp b/FINAL/DoAnMatlab.mlapp index 1de05bf..7f5e92d 100644 Binary files a/FINAL/DoAnMatlab.mlapp and b/FINAL/DoAnMatlab.mlapp differ diff --git a/FINAL/DoAnTimNghiem/chiadoi.m b/FINAL/DoAnTimNghiem/chiadoi.m index b280c3b..5189411 100644 --- a/FINAL/DoAnTimNghiem/chiadoi.m +++ b/FINAL/DoAnTimNghiem/chiadoi.m @@ -9,6 +9,6 @@ else b = c; end - n = n + 1; + n = n; end end \ No newline at end of file diff --git a/FINAL/DoAnTimNghiem/tieptuyen.m b/FINAL/DoAnTimNghiem/tieptuyen.m index 20dbd35..4785a02 100644 --- a/FINAL/DoAnTimNghiem/tieptuyen.m +++ b/FINAL/DoAnTimNghiem/tieptuyen.m @@ -1,10 +1,10 @@ -function [x3, n3] = tieptuyen(fx, a, b ,saiso) - n3 = 0; - x3 = (a+b)/2; - syms x - df = matlabFunction(diff(fx(x))); - while abs((x3-fx(x3)/df(x3))-x3) > saiso - x3 = x3 - fx(x3)/df(x3); - n3 = n3 + 1 ; - end -end \ No newline at end of file +function [x1, n] = tieptuyen(fx, a, b, saiso) + syms x; + df = matlabFunction(diff(fx(x))); % dao ham fx + x1 = (a + b)/2; + n = 0; + while abs(fx(x1)) > saiso +x1 = x1 - fx(x1)/df(x1); % newton +n = n + 1; +end +end \ No newline at end of file