vnames = strvcat('dpfo','w', 'zam','vpm');
y = vsetky; % use all eight states
nlag = 5;
pterm = 0;
result = johansen(y,pterm,nlag);
prt(result,vnames);