***************************** ** Oppgave 4 b >> format long >> g = inline('x - (x.^3 - 3*x - 2) ./ (3.*x.^2 - 3)'); >> p0 = 2.1; >> p1 = g(p0) p1 = 2.00606060606061 >> p2 = g(p1) p2 = 2.00002433978338 >> p3 = g(p2) p3 = 2.00000000039494 >> p4 = g(p3) p4 = 2 ********************************** ** Oppgave 22 a >> g = inline('x*(x^2+15)/(3*x^2+5)'); >> p0 = 2; >> p1 = g(p0) p1 = 2.23529411764706 >> p2 = g(p1) p2 = 2.23606797747661 >> p3 = g(p2) p3 = 2.23606797749979 >> >> sqrt(5) ans = 2.23606797749979 >> ******************************* ** Oppgave 22 b >> g = inline('(x^3+2*x^2+4*x+2)/(2*x^2+4*x+3)'); >> p0 = -2.4; >> p1 = g(p0) p1 = -2.01300813008130 >> p2 = g(p1) p2 = -2.00000072111960 >> p3 = g(p2) p3 = -2.00000000000000 >> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% Oppgave 23 d >> g=inline('x - (sin(2*x^3)*x^2)/(6*x^4-2*x*sin(2*x^3))'); >> p0 = 1; >> p1 = g(p0) p1 = 0.78253783237921 >> p2 = g(p1) p2 = 0.26558132223138 >> p3 = g(p2) p3 = 1.862855151226195e-04