Matlab - Problema programa interpolación Lagrange

   
Vista:

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

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
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder