1332/CH12/EX12.2/12_2.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20

//Example 12.2
//Lagarangian Method
//Page no. 373
clc;close;clear;

xk=[-1,0,2,5];
yk=[10,7,7,22];

P=0;
x=poly(0,"x");
for k=0:3
    p=yk(k+1)
    for j=0:3
        if(j~=k)
            p=p*((x-xk(j+1))/(xk(k+1)-xk(j+1)))
        end
    end
    P=P+p;
end
disp(P,'P=')