f x = x * x
n = 10000
integral f startX endX =
sum $ map (\x -> squareLen * f (x * squareLen)) [1..n]
where squareLen = (endX - startX) / n