function lab_q
global points

a = 0;
b = 1;

c = (a + b)/2;
h = (b - a)/2;
s = [.942882415695480 sqrt(2/3) .641853342345781 1/sqrt(5) .236383199662150];
points = [a c-h*s c c+h*fliplr(s) b];

disp(['*** Integralens värde med tol = 1e-6  är ', ...
           num2str(quadl(@f_lab, a, b)), '.'])
disp(['*** Integralens värde med tol = 1e-14 är ', ...
           num2str(quadl(@f_lab, a, b, 1e-14)), '.'])
disp('Här följer en en plot av integranden:')

figure
x = linspace(0, 1, 500);
plot(x, f_lab(x))
xlabel('x')
ylabel('y')
title('f(x)')
grid on

% -----------------------------------------

function y = f_lab(x)
global points

y = x;
for j = 1:length(x)
  y(j) = 1e5 * min(abs(x(j) - points))^2;
end
