use matlab to slove this question.

Write a program to evaluate I = integral^b_a f(x) dx using the Romberg algorithm. Compute seven rows in the Romberg array. (a) integral^1_0 e^-x^2 dx (b) integral^4_-4 1/1 + x^2 dx