b34sexec matrix; n=100000; a=array(n:); a=rn(a); b=rn(a); unitroot=cusum(a); broot=cusum(b); call graph(unitroot broot); call olsq(unitroot broot :print); call stataio(unitroot, broot :writefile 'test.dct'); b34srun;