summaryrefslogtreecommitdiff
path: root/964/CH21/EX21.7/21_7.sce
blob: 82bc6b812f2a3d7b5747a15ed6930262fab32035 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
clc;
clear;
function y=f(x)
    y=(0.2+25*x-200*x^2+675*x^3-900*x^4+400*x^5)
endfunction
tval=1.640533;
x=[0 0.12 0.22 0.32 0.36 0.4 0.44 0.54 0.64 0.7 0.8];
for i=1:11
    func(i)=f(x(i));
end
l=0;
for i=1:10
    l=l+(x(i+1)-x(i))*(func(i)+func(i+1))/2;
end
disp(l,"l=")
Et=tval-l;//error
et=Et*100/tval;//percent relative error
disp(Et,"The Error Et=")
disp("%",et,"The percent relative error et=")