diff options
| author | priyanka | 2015-06-24 15:03:17 +0530 |
|---|---|---|
| committer | priyanka | 2015-06-24 15:03:17 +0530 |
| commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
| tree | ab291cffc65280e58ac82470ba63fbcca7805165 /1673/CH3/EX3.23 | |
| download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip | |
initial commit / add all books
Diffstat (limited to '1673/CH3/EX3.23')
| -rwxr-xr-x | 1673/CH3/EX3.23/3_23.sce | 21 |
1 files changed, 21 insertions, 0 deletions
diff --git a/1673/CH3/EX3.23/3_23.sce b/1673/CH3/EX3.23/3_23.sce new file mode 100755 index 000000000..522722787 --- /dev/null +++ b/1673/CH3/EX3.23/3_23.sce @@ -0,0 +1,21 @@ +//example 3.22
+//newton's divided formula
+//page 114
+clc;clear;close
+x=[-1 0 3 6 7];
+y=[3 -6 39 822 1611];
+for i=1:4
+ d1(i)=(y(i+1)-y(i))/(x(i+1)-x(i));
+end
+for i=1:3
+ d2(i)=(d1(i+1)-d1(i))/(x(i+2)-x(i));
+end
+for i=1:2
+ d3(i)=(d2(i+1)-d2(i))/(x(i+3)-x(i));
+end
+for i=1:1
+ d4(i)=(d3(i+1)-d3(i))/(x(i+4)-x(i));
+end
+X=poly(0,'X')
+f_x=y(1)+(X-x(1))*(d1(1))+(X-x(2))*(X-x(1))*d2(1)+(X-x(1))*(X-x(2))*(X-x(3))*d3(1)+(X-x(1))*(X-x(2))*(X-x(3))*(X-x(4))*d4(1)
+disp(f_x,'the polynomial equation is =')
\ No newline at end of file |