function [p,cas]=pascal_t(n)

c1=cputime;
N=(n+1)*(n+2)/2; % celkovy rozsah
p=zeros(1,N);
p(1)=1;          % poc. hodnota
ki=1;            % koncovy index spocitane casti
pom=1;           % poc hodnota v pomocnem poli
for k=1:n
    pom=[0,pom]+[pom,0]; % vypocet dalsioho radku
    p(ki+1:ki+k+1)=pom;  % umisteni na misto
    ki=ki+k+1;           % oprava ki
end    
c2=cputime;
cas=c2-c1;