blob: 02a9092e48b3a86b58c19f8a84ddbea1a0e051de (
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
|
//Example 6.24
clc
clear
x = 1:3:10;
F = [500426 329240 175212 40365];
n = length(x);
del = %nan*ones(n,4);
del(:,1) = F';
for j = 2:4
for i = 1:n-j+1
del(i,j) = del(i+1,j-1) - del(i,j-1);
end
end
del0 = del(1,:);
X = 2;
x0 = x(1);
h = x(2) - x(1);
p = (X-x0) / h;
F2 = 0;
for i = 0:length(del0)-1
t = 1;
for j = 1:i
t = t * (p-j+1);
end
F2 = F2 + t*del0(i+1)/factorial(i);
end
f2 = F(1) - F2;
disp(f2,"f(2) = ")
|
