blob: 8d8e9ef963e36269a9e5fa12136af4714bcafc87 (
plain)
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
|
//Example 5.36
//Lagrange's Interpolation Method
//Page no. 176
clc;close;clear;
x=[7,8,9,10]
y=[3,1,1,9]
x0=9.5
printf('\tx\ty=f(x)\n-----------------------\n')
for i=1:4
printf('x%i\t%i\t %i\n',i-1,x(i),y(i))
end
p=1;p1=1;i=1;
for k=1:4
for j=1:4
if k~=j then
p=p*(x0-x(j))
p1=p1*(x(k)-x(j))
end
end
L(k)=p/p1
p=1;p1=1;
end
p=0;
for i=1:4
printf('\n L%i (x) = %g\n',i-1,L(i))
p=p+L(i)*y(i)
end
disp(p,"P(9.5) = ")
|
