¡Ayuda a compilar un programa matlab! Interpolación lagrangiana ¡Pintura del fenómeno Runge!

% Entrada: nodo de interpolación x, y (los puntos son diferentes entre sí, de lo contrario se modificarán como parámetros)

% punto de espera x0

% salida: x0, y0, resultado de interpolación

n = longitud (x% número de nodos de interpolación, xey deben ser); Coherente.

m = length(x0); porcentaje de puntos a encontrar

y0 = zero(1,m); % asigna espacio para la salida

For j = 1:m

Para i=1:n

% Fórmula de interpolación lagrangiana: li (x) = yi * prod (x-xj)/prod (xi-xj), I ~ = j;

XJ = x;

XJ(I)=[];

y0(j)= y0(j)+ y(I) )* prod(x0(j)-XJ)/prod(x(I)-XJ);

Fin

Fin

%% %%% %%%%%%%%%

Método de llamada:

x =[-5:5];

y = 1. /(1+x.^2);

x0 =[-5:0.01:5];%Elija un tamaño de paso apropiado según el rendimiento de su computadora.

[x0, y0] = Lagrangiano_2(x, y, x0);

Espera

plot(x0, 1. /(1+x0); .^2),'r');%función original y = 1/(1+x^2) imagen

plot(x0, y0);%función de interpolación imagen