%Método de Newton - Raphson aplicando la derivada al inicio clc clear pkg load symbolic syms x disp("Método de Newton - Raphson") disp ("Tabla de valores de iteración por Newton Raphson") format short N=input("Ingrese el número aproximado de interaciones "); %N=100; tole=input("Ingrese la tolerancia "); %tol=10^(-6); f=input("Ingrese la ecuación a calcular \n"); ezplot(f) grid on x0=input("Ingrese el valor aproximado inicial "); %x0=1 ; fa=input("Ingrese la ecuación en forma anónima \n"); df1=diff(f); disp("La primera derivada es "); disp(df1) df1a=input("Ingrese la primera derivada en forma anónima\n"); df2=diff(df1); disp("La segunda derivada es ") disp(df2) df2a=input("Ingrese la segunda derivada en forma anónima\n"); disp(" i x abs x abs g") for I=1 : N x=x0-(f(x0))/(df1(x0)); dist = abs(x - x0); dg=abs(1-((df1a(x)^2-(fa(x)*df2a(x)))/df1a(x)^2)); fprintf(" %i %f %f %f\n",I,x,dist,dg) if abs(x-x0)