stacksize(10000000);
N=5000;
timer();
s=0.0;
for i=1:N
  s=s+i;
end;
t1=timer()
s = ones(1,N)*(1:N)';
t2 = timer()

