Universidad de San Buenaventura
Laboratorio 2do corte - Métodos Numéricos
Integrantes
• Ana María Molano Ballesteros - 30000125688
• Andrés Felipe Ramírez Arias - 30000125515
• Miguel Ángel Rosas Rodríguez - 30000125288
Ejercicios
1. Encuentre los polinomios de Taylor con n=0,1,2,3 para la función en x = 0. Luego, verifique sus
resultados en Matlab con el polinomio escrito como una serie. Adicionalmente realizar el grafico de la función y sus
respectivos polinomios de aproximación para n = 0,1,2,3.
syms x
f = exp(x) * sin(x);
a = 0;
n = 3; % Máximo orden del polinomio de Taylor
% Inicializar el polinomio de Taylor con el término de orden 0
P(x) = subs(f, x, a);
% Configurar la gráfica
fplot(f, [-1, 1], 'k', 'LineWidth', 2); % Función original en negro
hold on;
grid on;
colors = ['r', 'b', 'g', 'm']; % Colores para cada polinomio de Taylor
legends = {'Función original'};
% Construcción de los polinomios de Taylor para n = 0,1,2,3
for n = 0:n_max
P(x) = subs(f, a); % Inicializar con f(0)
for k = 1:n
% Calcular la k-ésima derivada
df = diff(f, x, k);
% Evaluar la derivada en x = a
df_a = subs(df, x, a);
% Agregar el término al polinomio de Taylor
P(x) = P(x) + (df_a/factorial(k))*(x - a)^k;
end
P(x) = expand(P(x))
% Graficar cada polinomio de Taylor
fplot(P(x), [-1, 1], colors(n+1), 'LineWidth', 1.5);
legends{end+1} = sprintf('T_%d(x)', n);
end
P(x) =
0
P(x) =
x
P(x) =
P(x) =
El polinomio resultante es:
Pn(x)=
% Configuración de la gráfica
xlabel('x');
ylabel('f(x) y T_n(x)');
legend(legends, 'Location', 'Best');
title('Polinomios de Taylor de sqrt(3x+1) en x=0');
hold off;
2. Encuentre los polinomios de Taylor con n=0,1,2,3 para la función en x = 0. Luego, verifique sus
resultados en Matlab con el polinomio escrito como una serie. Adicionalmente realizar el grafico de la función y sus
respectivos polinomios de aproximación para n = 0,1,2,3
syms x
f = sqrt(3*x + 1);
a = 0;
n = 3; % Máximo orden del polinomio de Taylor
% Inicializar el polinomio de Taylor con el término de orden 0
P(x) = subs(f, x, a);
% Configurar la gráfica
fplot(f, [-1, 1], 'k', 'LineWidth', 2); % Función original en negro
hold on;
grid on;
colors = ['r', 'b', 'g', 'm']; % Colores para cada polinomio de Taylor
legends = {'Función original'};
% Construcción de los polinomios de Taylor para n = 0,1,2,3
for n = 0:n
P(x) = subs(f, a); % Inicializar con f(0)
for k = 1:n
% Calcular la k-ésima derivada
df = diff(f, x, k);
% Evaluar la derivada en x = a
df_a = subs(df, x, a);
% Agregar el término al polinomio de Taylor
P(x) = P(x) + (df_a/factorial(k))*(x - a)^k;
end
P(x) = expand(P(x))
% Graficar cada polinomio de Taylor
fplot(P(x), [-1, 1], colors(n+1), 'LineWidth', 1.5);
legends{end+1} = sprintf('T_%d(x)', n);
end
P0(x) =
1
P1(x) =
P2(x) =
P3(x) =
El polinomio resultante es:
Pn(x)=
% Configuración de la gráfica
xlabel('x');
ylabel('f(x) y T_n(x)');
legend(legends, 'Location', 'Best');
title('Polinomios de Taylor de sqrt(3x+1) en x=0');
hold off;