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 /191/CH5/EX5.7 | |
| 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 '191/CH5/EX5.7')
| -rwxr-xr-x | 191/CH5/EX5.7/Example5_7.sce | 40 | ||||
| -rwxr-xr-x | 191/CH5/EX5.7/Result5_7.txt | 29 |
2 files changed, 69 insertions, 0 deletions
diff --git a/191/CH5/EX5.7/Example5_7.sce b/191/CH5/EX5.7/Example5_7.sce new file mode 100755 index 000000000..223cce35e --- /dev/null +++ b/191/CH5/EX5.7/Example5_7.sce @@ -0,0 +1,40 @@ +//Illustration of Central Difference Formula +close(); +clear; +clc; +xi = 0:0.2:1.2; +fi = sin(xi); +x0 = 0; +h = 0.2; +format('v',8); +// First order difference +delta1_fi = diff(fi); +// Second order difference +delta2_fi = diff(delta1_fi); +// Third order difference +delta3_fi = diff(delta2_fi); +// Fourth order difference +delta4_fi = diff(delta3_fi); +//Fifth order difference +delta5_fi = diff(delta4_fi); +//Sixth order difference +delta6_fi = diff(delta5_fi); +disp(fi , 'Values of f(x) : ') +disp(delta1_fi , 'First Order Difference :') +disp(delta2_fi , 'Second Order Difference :') +disp(delta3_fi , 'Third Order Difference :') +disp(delta4_fi , 'Fourth Order Difference :') +disp(delta5_fi , 'Fifth Order Difference :') +disp(delta6_fi , 'Sixth Order Difference :') +//Calculating p2(0.67) +xm = 0.6; +x = 0.67; +s = (x-xm)/0.2; +p2 = fi(4) + {s*(delta1_fi(3)+delta1_fi(4))/2} + s*s*(delta2_fi(3))/2; +disp(p2 , 'Value of p2(0.67) : '); +//Calculating p4(0.67) +p4 = p2 + s*(s*s-1)*(delta3_fi(3)+delta3_fi(2))/12 + s*s*(s*s-1)*delta4_fi(2)/24; +disp(p4 , 'Value of p4(0.67) : '); +//Exact value of sin(0.67) is 0.62099 so error in estimation +err = 0.62099-0.62098; +disp(err , 'Error in estimation : ');
\ No newline at end of file diff --git a/191/CH5/EX5.7/Result5_7.txt b/191/CH5/EX5.7/Result5_7.txt new file mode 100755 index 000000000..689644837 --- /dev/null +++ b/191/CH5/EX5.7/Result5_7.txt @@ -0,0 +1,29 @@ + First Order Difference :
+
+ 0.19867 0.19075 0.17522 0.15271 0.12411 0.09057
+ Second Order Difference :
+
+ - 0.00792 - 0.01552 - 0.02251 - 0.02860 - 0.03355
+ Third Order Difference :
+
+ - 0.00760 - 0.00699 - 0.00609 - 0.00495
+ Fourth Order Difference :
+
+ 0.00062 0.00090 0.00114
+ Fifth Order Difference :
+
+ 0.00028 0.00024
+
+ Sixth Order Difference :
+
+ - 0.00004
+ Value of p2(0.67) :
+
+ 0.62065
+ Value of p4(0.67) :
+
+ 0.62098
+
+ Error in estimation :
+
+ 0.00001
\ No newline at end of file |