
problem1 = @(x) sqrt(x);
problem1a = @(y) 2*y^2;
problem2 = @(x) 1./(1+x.^2);
problem3 = @(x) 1/pi*cos(8*sin(x)-x);
problem4 = @(x) 1/(1+.5*sin(2*pi*x));
problem5 = @(x) exp(x.^2*(1-x)^2);

A2 = Romberg(problem1,0,1,2/3,12);
A = Romberg(problem1a,0,1,2/3,12);
B = Romberg(problem2,0,1,pi/4,12);
C = Romberg(problem3,0,pi,0.234636347,12);
D = Romberg(problem4,0,1,1.15470054,12);
E = Romberg(problem5,0,1,1.03414105,12);
