function [ vc3 ] = nasob_vc_d( vc1,vc2 )
% function [ vc3 ] = nasob_vc_d( vc1,vc2 )
% nasobeni velkych cisel reprezentovanych 
% polem dvojic cifer uloznych v double

d1=length(vc1);
d2=length(vc2);

l=6;
base=10^l;
if d1==1 & vc1==0
    vc3=vc1;
    return;
end
if d2==1 & vc2==0
    vc3=vc2;
    return;
end

vc3=0;

for k1=d1:-1:1
    pc=zeros(1,d2); % pomocne cislo pro mezivysledky
    prenos=0;
    for k2=d2:-1:1
        a=vc1(k1)*vc2(k2)+prenos;
        b=mod(a,base);
        prenos=(a-b)/base;
        pc(k2)=b;
    end
    if prenos>0
        pc=[prenos, pc];
    end
    pc=[pc zeros(1,d1-k1)];
    vc3=secti_vc_d(vc3,pc);
end


end

