Как векторизовать вложенные циклы в октаве?

Question

Мне нужно ускорить свою функцию, и я ищу эффективный способ вычисления этих двух циклов ниже. Любая помощь будет оценена.

  for i=1:n
    for j=1:m
      ig = i + (j-1)*n;
      coord(ig,1) = i;
      coord(ig,2) = j;
    end
  end

  for j=1:m
    for i=1:n-1
      k = (j-1)*(n-1) + i;
      conec(k,1) = (j-1)*n + i;
      conec(k,2) = (j-1)*n + i+1;
      C(k,1) = CH;
    end
  end

   for k=1:nc
     p = conec(k,1);
     q = conec(k,2);
     aux = [C(k) -C(k) ; -C(k) C(k)];
     A([p q],[p q]) =  A([p q],[p q]) + aux;
   end