function [A,b] = priklad(h)

a = 1;
b = 1;
[uzly,prvky,ID] = diskretizuj_obdelnik(a,b,h);
nuzly = size(uzly,2);
nprvky = size(prvky,2);

%f = ones(1,nprvky);
f = zeros(1,nprvky);
for i=1:nprvky
    massPoint = (1/3) * sum(uzly(:,prvky(:,i)),2);
    x = massPoint(1); y = massPoint(2);
    f(i) = 1*sin(pi*x)*sin(pi*y) + 1*sin(4*pi*x)*sin(4*pi*y) + 10*sin(20*pi*x)*sin(20*pi*y);
end
[A,b] = mkp_2d(uzly,prvky,f,ID);
u = A\b;
vykresli_reseni(uzly,prvky,u);

