Problema programa interpolación Lagrange
Publicado por Jose Ramon (1 intervención) el 02/05/2012 22:41:26
Buenas noche, al poner en marcha el programa el resultado del polinomio de interpolación (p) me da como resultado siempre infinito, y no encuentro el fallo
Gracias de antemano
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
clear
syms x p l Y
k=input('introduce el numero de nodos:');
M(k,k+1)=0;
for L=1:1:k
fprintf('Introduce el numero%d:', L)
f=input('')
for H=1:1:(k+1)
M(L,1)=1;
M(L,H)=f^(H-1);
end
end
M
for L=1:1:(k+1)
X(L)=x^(L-1);
end
X
d=M(:,2);
Punto=d';
for L=1:1:k
fprintf('dame el punto %d a dibujar', Punto(L))
Y(L)=input('');
end
Y
l=X/M;p=l*Y
disp(expand(p));
ezplot(p,[Y(1),Y(k)])
hold on, plot(Puntuak,Y,'s')
hold off
Gracias de antemano
Valora esta pregunta
0