50/CH5/EX5.13/ex_5_13.sce


1
2
3
4
5
6
7
8
9
10
11
12

// example 5.13
// caption: simpson 3-8 rule


// let integration of f(x)=1/(1+x) in the range [0,1] by simpson 3-8 rule is equal to I

x=0:1/3:1;
deff('[y]=f(x)','y=1/(1+x)')

[I] = simpson38(x,f)