blob: 05dca459d6d4d153cde9ee5b47e70226a4f93441 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
//Example 6.15
clc
clear
x = [1 2 5];
y = [1 4 10];
n = length(x);
Y = 0;
X = poly(0, "X");
//X = 5;
for i = 1:n
t = x;
t(i) = [];
p = 1;
for j = 1:length(t)
p = p * (X-t(j))/(x(i)-t(j));
end
Y = Y + p*y(i);
end
Y5 = horner(Y,3);
disp(Y5,"f(3) = ")
|
