blob: 2997bd6c5d492a047124fa790d206dfb12a30de7 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
|
function [L012,L02,L01]=iteratedinterpol (x,f,n)
X=poly(0,"X");
L01=(1/(x(2)-x(1)))*det([f(1) x(1)-X;f(2) x(2)-X]);
L02=(1/(x(3)-x(1)))*det([f(1) x(1)-X;f(3) x(3)-X]);
L012=(1/(x(3)-x(2)))*det([L01 x(2)-X;L02 x(3)-X]);
endfunction
|
