function x=diagsol(A,b)
  // Mat01032 - Calculo Numerico A - UFRGS
  // solucao de Ax=b quando A eh diagonal e nao-singular
  x=zeros(b); n=size(A,1);
  if n~=size(A,2) 
    error('matriz A deve ser quadrada');
  end
  if n~=length(b)
    error('A e b tem dimensoes incompativeis');
  end
  for i=1:n
    x(i)=b(i)/A(i,i);
  end
  // numero de operacoes ponto flutuante: n
endfunction
