Simpsons Rule

[Jeff_Emanuel]

@lmiq

function simps(f::Function,a,b,n)
    h=(b-a)/(2*n)
    x=a:h:b
    return simps(x,a,b)
end 

It appears you dropped f