% Datos de entrada
x = [1, 2, 3, 4, 5];
y = [1, 4, 9, 16, 25];
% Punto a evaluar
x_eval = 8;
% Número de puntos
n = length(x);
% Tabla de diferencias divididas
div_diff = zeros(n, n);
div_diff(:, 1) = y';
% Calcular las diferencias divididas
for j = 2:n
for i = 1:(n-j+1)
div_diff(i, j) = (div_diff(i+1, j-1) - div_diff(i, j-1)) / (x(i+j-1)
- x(i));
end
end
% Coeficientes del polinomio de Newton
coeff = div_diff(1, :);
% Evaluar el polinomio de Newton en x_eval
P = coeff(1);
product_term = 1;
for k = 2:n
product_term = product_term * (x_eval - x(k-1));
P = P + coeff(k) * product_term;
end
% Mostrar el resultado
fprintf('El valor de P(%d) es %.2f\n', x_eval, P);
El valor de P(8) es 64.00
Published with MATLAB® R2024a
1
HECHO POR: MARCELO EMILIO AZPILCUETA RUBINA