function z = VC_minus(x,y)
% function z = VC_minus(x,y)
% Odecitani velkych cisel
% x,y,z - pole uint64
% Cisla jsou ulozena po 9 cifrach, 
% nejmensi mocniny jsou na zacatku pole.
% Predpoklada se, ze x>y

h=1e9;    % hranice pro 9 cifer
prenos=0; % prenos do dalsi polozky pri odecitani od mensioho cisla
nx=length(x);
ny=length(y);
if nx<ny  % doplnime kratsi cislo nulami
    for k=nx+1:ny
        x(k)=0;
    end
else
    for k=ny+1:nx
        y(k)=0;
    end
end
N=max([nx,ny]);
z=uint64(zeros(1,N));
for k=1:N
    yy=y(k)+prenos;
    if yy>x(k)
        v=h+x(k)-yy;
        prenos=1;
    else
        v=x(k)-yy;
        prenos=0;
    end
    z(k)=v;
end
k=N;
while z(k)==0 && k>1
    z(k)=[];
    k=k-1;
end
        
end

