% primer 1
f = @(r) 1+0*r;
R = 5;
n = 10;
[r,u] = upogib_opne(f,R,n);

% >> u
% u =
%   -6.383328470634976
%   -6.258328470634976
%   -6.049995137301643
%   -5.724995137301643
%   -5.278566565873071
%   -4.709122121428625
%   -4.015940303246806
%   -3.198632610939113
%   -2.256965944272447
%   -1.190789473684211
%                    0


% primer 2
f = @(r) 1-r.^2;
R = 2;
n = 20;
[r,u] = upogib_opne(f,R,n);

% u(1)
% ans =
%   -0.009941434047129


% primer 3
f = @(r) sin(2*pi*r);
R = 3;
n = 50;
[r,u] = upogib_opne(f,R,n);

% u(1)
% ans =
%   -0.038213752522659