Newton raphson en cas y matlab

Codificación por medio de la calculadora TI-nspire cx cas
Define LibPub nr()=
Prgm
:Request "dame la funcion evaluada en x",fun1
:Request "dame el valor inicial",x0
:Request "dameel error maximo",error
:For i,1,10
:x1:=x0-((fun1|x=x0)/((fun1,x)|x=x0))
:errorc:=abs(x1-x0)
:If errorc≤error Then
:Disp "tu resultado es",x1,"en la iteracion",i
:Disp"error=",errorc
:Exit
:Else
:x0:=x1
:EndIf
:EndFor
:EndPrgm
Corrido del programa en calculadora TI-nspire cx cas
dame la funcion evaluada en x((1)/(0.012))*((4*x+4*x^(2))/(2))*(((((4*x+4*x^(2))/(2)))/(2+2*√(x^(2)+2*x^(2)))))^(((2)/(3)))*(0.5)^(((1)/(2)))-4
dame el valor inicial 1
dame el error maximo 0.0001
tu resultado es 0.1262 en la iteracion 6.0000error= 2.8197−5
Codificación MATLAB R2010a
syms x,
fun1=input('Introduce la función evaluada en x ');
x0=input('Introduce el valor de x0 ');
error=input('Introduce el error máximo ');funcion=inline(fun1);
der=diff(fun1);
derivada=inline(der);
for i=1:10;
fx=feval(funcion,x0);
dfx=feval(derivada,x0);
x1=x0-(fx/dfx);
error1=(abs(x1-x0));if(error1<=error)
break
else x0=x1;
end
end
fprintf(' el valor de x= %6.5f con un error de %6.5f\n',x0,error1)

Corrido del programa en MATLAB R2010a
Introduce lafunción evaluada en x (((1)/(0.012))*((4*x+4*x^(2))/(2))*(((((4*x+4*x^(2))/(2)))/(2+2*sqrt(x^(2)+2*x^(2)))))^(((2)/(3)))*(0.5)^(((1)/(2)))-4)
Introduce el valor de x0 1
Introduce el errormáximo 0.0001
el valor de x= 0.12626 con un error de 0.00003

Conclusión
Se utilizo el método de Newton Rhapson por que fue el método más sencillo de programar y de utilizarpara cualquier usuario, además de que nos arroja un resultado muy certero en todos los programas en los que se codifico el sistema de solución.
Pérez Gámez Francisco Antonio