res=load('laine.out');
[n1,n2]=size (res);
n=n1/n2;

for i =1:n
clf   
surfc(res(1+(i-1)*n2:i*n2,:))
pause (0.1)
end