function RP = fun(x,a)
% funkcija vrne vrednosti na robu kvadrata x krat x
% vrednosti na stranicah zlozi v vrstice 
n = length(x);
%RP = [ zeros(1,n); 1-x.^2; 2-2*x.^2; zeros(1,n)];
RP=[a^2-x.^2;a^2-x.^2;a^2-x.^2;zeros(1,n)];
%RP=zeros(4,n);
