function [spojeni,matice,vektor] = vicepromennych(pocpoz,b)         %funkce pro model s vice vysvetlujicimi promennymi bude vracet tri parametry a jsou ji predany dva
                                                                    %nasleduje obdobny vypis jako ve funkci jednapromenna.m
spojeni = ddeinit('excel','linearniregrese.xls');
if spojeni==0
           Vypis=str2mat(... 
        'Není otevøen soubor linearniregrese.xls pro naètení dat!');
        set(findobj('Tag','List'),'String',Vypis); 
        error('Není otevøen soubor linearniregrese.xls pro naètení dat!');    
end

X=ddereq(spojeni, 'r32c12:r131c15');
y=ddereq(spojeni, 'r32c11:r131c11');

c=size(X);
d=size(y);
if c(2)~=4 || d(2)~=1
               Vypis=str2mat(... 
        'Špatnì zadané vstupní hodnoty v souboru linearniregrese.xls! Používejte desetinnuo teèku!');
        set(findobj('Tag','List'),'String',Vypis); 
        error('Špatnì zadané vstupní hodnoty v souboru linearniregrese.xls! Používejte desetinnuo teèku!'); 
end

if c(1)<100
        X(c(1)+1 : 100,:)=zeros;
end
if length(y)<100
        y(length(y)+1 : 100,:)=zeros;
end
if pocpoz<100
    X(pocpoz+1 : 100,:)=[];
    y(pocpoz+1 : 100,:)=[];
end
x=4;                                                                %naplneni matice X prislusnymi sloupci matice Z, jednodussi oproti funkci jednapromenna.m, protoze mame primo nactenou matice dat ruznych vysvetlujicich a pouze je priradime podle zadaneho modelu
for i=4:-1:1
   if b(i)==1
     X(:,i)=[];
     x=x-1;
   end
end

X=[ones(length(X),1) X];                                

[n,k]=size(X);
if x~=k-1
           Vypis=str2mat(... 
        'Poèet vysvìtlujících promìnných musí být stejný jako poèet parametrù!');
        set(findobj('Tag','List'),'String',Vypis); 
        error('Poèet vysvìtlujících promìnných musí být stejný jako poèet parametrù!'); 
end

matice = X;
vektor =y;