1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
//Example 5.17
//Error Propagation
//Page no. 140
clc;close;clear;
printf(' x\t y\t\tdy\td2y\t d3y\t d4y\t d5y\n')
printf('------------------------------------------------------------------------------------------')
x=[1,1;1.1,1.5191;1.2,2.0736;1.3,2.6611;1.4,3.2816;1.5,3.9375;1.6,4.6363;1.7,5.3771;1.8,6.1776;1.9,7.0471;2,8]
for i=3:7
for j=1:13-i
x(j,i)=x(j+1,i-1)-x(j,i-1)
end
end
disp(x)
for i=1:11
if abs(x(i,7))<10^-5 then
continue
else
break
end
end
printf("\n\Therefore the error is in the value corresponding to %g i.e. %g",x(i+5,1),x(i+5,2))
|
