Question: Write a scilab program to integrate x^2 + 2x from 1 to 5 by using Simpson's 1/3 and 3/8 Rule

Simpson's 1/3 rule :

Formula : Integral of f(x) limit a to b =

h/3[f(a)+4[f(x1)+f(x3)....f(xn-1)]+2[f(x2)+f(x4)....f(xn-2)+f(b)]

Here in the question 

F(x)=x^2+2x

H=b-a/n

N= no.of value of x 

a=1

b=5