blob: b7c5b84e561dd2f0716c407384de3769f621712d (
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
30
31
32
33
34
35
|
//Example 6.13
clc
clear
x = 1974:2:1982;
y = [40 43 48 52 57];
n = length(x);
del = %nan*ones(n,5);
del(:,1) = y';
for j = 2:5
for i = 1:n-j+1
del(i+j-1,j) = del(i+j-1,j-1) - del(i+j-2,j-1);
end
end
X = 1979;
h = x(2) - x(1);
p = (X-x(n)) / h;
xn = x(n);
yn = y(n);
delyn = del(n,:);
Y = 0;
for i = 0:length(delyn)-1
t = 1;
for j = 1:i
t = t * (p+j-1);
end
Y = Y + t*delyn(i+1)/factorial(i);
end
Y = round(Y*10^4)/10^4;
disp(Y,"Estimated sales for the year 1979: ")
|
