function [Uzly,Segmenty] = diskretizace_magnetostatika_MHP(sirkaOmegar,vyskaOmegar,h)

nr = 2 * ( floor(vyskaOmegar/h) + floor(sirkaOmegar/h) );

Uzly = zeros(2, nr );
Segmenty{1} = [ [1:nr] ; [2:nr,1] ];                      % Gammar

%
% Uzly Gama^r
%
x = [sirkaOmegar/2;-vyskaOmegar/2];
% smer nahoru
n = floor(vyskaOmegar/h);
Uzly(:,1) = x;
idxUzly = 2;
for i=1:n
    x(2) = x(2) + h;
    Uzly(:,idxUzly) = x;
    idxUzly = idxUzly+1;
end
% smer doleva
n = floor(sirkaOmegar/h);
for i=1:n
    x(1) = x(1) - h;
    Uzly(:,idxUzly) = x;
    idxUzly = idxUzly+1;
end
% smer dolu
n = floor(vyskaOmegar/h);
for i=1:n
    x(2) = x(2) - h;
    Uzly(:,idxUzly) = x;
    idxUzly = idxUzly+1;
end
% smer doprava
n = floor(sirkaOmegar/h);
for i=1:n-1
    x(1) = x(1) + h;
    Uzly(:,idxUzly) = x;
    idxUzly = idxUzly+1;
end
