SoloCodigo

quisiera sabe rqu edebo hacer para que fujncione este algortimo?
me da el sig. error:
??? Input argument 'A' is undefined.

Error in ==> C:MATLAB6p5workGaussSeidel.m
On line 16  ==> n = size(A,1);

function x = GaussSeidel(A,b,varargin)

% Implementacion del metodo Gauss-Seidel para la solucion de sistemas
% de ecuaciones, tomando como aproximacion inicial x0.
%
%    x = GaussSeidel(A,b,x0=zeros,eps=0.001,nrm=Inf)
%
% Regresa x, la solucion del sistema Ax=b. El criterio de terminacion es
% que norm(x-xAnt,nrm)/norm(x,nrm)<eps.

% 21 agosto 2007
% Manuel Valenzuela

MAXCICLOS = 1000;
MAXPIVOTE = 1;
n = size(A,1);
x = zeros(n,1);
eps = 0.001;
nrm = Inf;
if length(varargin)>=1
   x = varargin{1};
end
if length(varargin)>=2
   eps = varargin{2};
end
if length(varargin)>=3
   nrm = varargin{3};
end
A = [A b];

if MAXPIVOTE==0
   for i=1:n-1
      % Encontrar renglon del maximo pivote
      k = find(abs(A(:,i))==max(abs(A(i:n,i))),1,'last');

      if k~=i
         % Intercambiar renglones
         rPivote = A(k,:);
         A(k,:) = A(i,:);
         A(i,:) = rPivote;
      end
   end
end

b = A(:,n+1);
A = A(:,1:n);

% Se obtienen el vector c y la matriz T
c = b./diag(A);
T = zeros(n);
for i=1:n
   T(i,:) = A(i,:)/A(i,i);
end
T = -T+eye(n);

xAnt = x;
j = 0;
while 1
   j = j+1;
   if j>MAXCICLOS
       warning('numero maximo de ciclos excedido')
       break
   end
   for i=1:n
      x(i) = T(i,:)*x + c(i);
   end
   if (norm(x-xAnt,nrm)/norm(x,nrm)<eps)
      break
   end
   xAnt = x;
end