// object mat tyle a=tlist(['mat','dim'],...)
function x=mat(m,n)
   s='tlist(["'mat'","'dim'"],[m,n]';
   d=m*n;
   for i=1:d, s=s+',0', end;
   s=s+')';
   x=evstr(s);
endfunction


function x=m_e(a,m,n)
   d=a('dim');
   s=(m-1)*d(2)+n+2;
   x=a(s);
endfunction
 

function x=m_d(a,m,n,b)
   d=a('dim');
   s=(m-1)*d(2)+n+2;
   x=a;
   x(s)=b
endfunction


function x=%mat_a_mat(a,b)
   c=a('dim');
   d=c(1)*c(2);
   x=a;
   for i=3:2+d, x(i)=a(i)+b(i), end;
endfunction


function x=%mat_s_mat(a,b)
   c=a('dim');
   d=c(1)*c(2);
   x=a;
   for i=3:2+d, x(i)=a(i)-b(i), end;
endfunction


function x=%mat_m_mat(a,b)
   c=a('dim');
   d=b('dim');
   x=mat(c(1),d(2));
   for i=1:c(1), 
      for j=1:d(2),
         ii=(i-1)*d(2)+j+2; 
         x(ii)=a(2+(i-1)*c(2)+1)*b(2+(j-1)*d(1)+1); 
         for k=2:c(2), 
            x(ii)=x(ii)+a(2+(i-1)*c(2)+k)*b(2+(j-1)*d(1)+k);
         end;
      end;
   end;
endfunction
   

function x=%mat_s(a)
   c=a('dim');
   cc=c(1)*c(2);
   x=a;
   for i=1:cc,
      x(i+2)=(-1)*a(i+2);
   end;
endfunction


function x=%s_m_mat(s,a)
   c=a('dim');
   x=a;
   for i=1:c(1),
      for j=1:c(2),
          ii=(i-1)*c(2)+j+2; 
          x(ii)=s*a(ii);
      end;
   end;
endfunction


function x=%mat_m_s(a,s)
   c=a('dim');
   x=a;
   for i=1:c(1),
      for j=1:c(2),
          ii=(i-1)*c(2)+j+2;
          x(ii)=s*a(ii);
      end;
   end;
endfunction
   


function x=%intval_m_mat(s,a)
   c=a('dim');
   x=a;
   for i=1:c(1),
      for j=1:c(2),
          ii=(i-1)*c(2)+j+2; 
          x(ii)=s*a(ii);
      end;
   end;
endfunction


function x=%mat_m_intval(a,s)
   c=a('dim');
   x=a;
   for i=1:c(1),
      for j=1:c(2),
          ii=(i-1)*c(2)+j+2;
          x(ii)=s*a(ii);
      end;
   end;
endfunction


function x=%mat_r_s(a,s)
   c=a('dim');
   x=a;
   for i=1:c(1),
      for j=1:c(2),
          ii=(i-1)*c(2)+j+2;
          x(ii)=a(ii)/s;
      end;
   end;
endfunction