1 function mat=next_mat(f)
\r
2 %mat=reshape(a(k:k+mm*nn),mm,nn);
\r
3 [sz,count]=fscanf(f,'%g',[1 2]);
\r
8 error('not enough terms in size')
\r
10 %fprintf('reading matrix size %g %g\n',sz)
\r
11 if any(any(sz~=round(sz) | sz < 0)),
\r
14 [mat,count]=fscanf(f,'%g',sz);
\r
15 if(count<sz(1)*sz(2)),
\r
16 error('not enough terms in matrix')
\r