function [x, Ex]=xfel(A,b,Ea,Eb);
% [x, Ex]=xfel(A,b,Ea,Eb) beräknar störningen Ex i lösningen till det
%linjära ekvationssystemet A x = b. Störningen i koefficientmatrisen
% A är Ea (matris) och störningen i högerledet b är Eb(vektor).
Ex=zeros(size(b));
nm=size(A); n=nm(1);
Wz=A\[eye(n) b];
x=Wz(:,n+1);
for i=1:n,
  for j=1:n,
   if Ea(i,j) ~= 0,
    s=Ea(i,j)*x(j)/(1+Ea(i,j)*Wz(j,i));
    ex=abs(s*Wz(:,i));
    Ex=Ex+ex;
   end %if Ea(i,j) <> 0
  end %j
  if Eb(i) ~= 0,
    ex=abs(Eb(i)*Wz(:,i));
    Ex=Ex+ex;
  end %if Eb(i) <> 0
end %i



