Variable 'k' is not fully defined on some execution paths.

function symfault1(zdata, Zbus)
persistent Ink Inkm Ign

nl = zdata(:,1); nr = zdata(:,2); R = zdata(:,3);
X = zdata(:,4);
nc = length(zdata(1,:));
  if nc > 4
  BC = zdata(:,5);
  elseif nc ==4, BC = zeros(length(zdata(:,1)), 1);
  end
ZB = R + j*X;
nbr=length(zdata(:,1)); nbus = max(max(nl), max(nr));

 V0 = ones(nbus, 1) + j*zeros(nbus, 1);
 Vf = ones(nbus, 1) + j*zeros(nbus, 1);

fprintf('Three-phase balanced fault analysis \n')
ff = 999;
while ff > 0
nf =3;
     while nf <= 0 | nf > nbus
     fprintf('Faulted bus No. must be between 1 & %g \n', nbus)
     nf = 3;
     end
fprintf('\nEnter Fault Impedance Zf = R + j*X in ')
Zf = 0+j*0.16;
fprintf('  \n')
fprintf('Balanced three-phase fault at bus No. %g\n', nf)

If = V0(nf)/(Zf + Zbus(nf, nf));
Ifm = abs(If); Ifmang=angle(If)*180/pi;
fprintf('Total fault current = %8.4f per unit \n\n', Ifm)
%fprintf(' p.u. \n\n', Ifm)
fprintf('Bus Voltages during fault in per unit \n\n')
fprintf('     Bus     Voltage       Angle\n')
fprintf('     No.     Magnitude     degrees\n')

for n = 1:nbus
    if n==nf
    Vf(nf) = V0(nf)*Zf/(Zf + Zbus(nf,nf)); Vfm = abs(Vf(nf)); angv=angle(Vf(nf))*180/pi;
    else, Vf(n) = V0(n) - V0(n)*Zbus(n,nf)/(Zf + Zbus(nf,nf));
    Vfm = abs(Vf(n)); angv=angle(Vf(n))*180/pi;
    end
    fprintf('   %4g',  n), fprintf('%13.4f', Vfm),fprintf('%13.4f\n', angv)

end

fprintf('  \n')

fprintf('Line currents for fault at bus No.  %g\n\n', nf)
fprintf('     From      To     Current     Angle\n')
fprintf('     Bus       Bus    Magnitude   degrees\n')


for n= 1:nbus
   %Ign=0;
    for I = 1:nbr
      if (nl(I) == n | nr(I) == n)
         if (nl(I) ==n)   
             k = nr(I);
         elseif( nr(I) == n)
             k = nl(I);
         end
           if (k==0)
            Ink = (V0(n) - Vf(n))/ZB(I);
            Inkm = abs(Ink); th=angle(Ink);
               %if th <= 0
                if real(Ink) > 0
                fprintf('      G   '), fprintf('%7g',n), fprintf('%12.4f', Inkm)
                fprintf('%12.4f\n', th*180/pi)
                elseif real(Ink) ==0 & imag(Ink) < 0
                fprintf('      G   '), fprintf('%7g',n), fprintf('%12.4f', Inkm)
                fprintf('%12.4f\n', th*180/pi)
               else, end
           Ign=Ink;
           elseif (k ~= 0)
            Ink = (Vf(n) - Vf(k))/ZB(I)+BC(I)*Vf(n);
            %Ink = (Vf(n) - Vf(k))/ZB(I);
            Inkm = abs(Ink); th=angle(Ink);
            %Ign=Ign+Ink;
            %if th <= 0
             if real(Ink) > 0
                  fprintf('%7g', n), fprintf('%10g', k),
                  fprintf('%12.4f', Inkm), fprintf('%12.4f\n', th*180/pi)
                  elseif real(Ink) ==0 & imag(Ink) < 0
                  fprintf('%7g', n), fprintf('%10g', k),
                  fprintf('%12.4f', Inkm), fprintf('%12.4f\n', th*180/pi)
                  else, end
              else, end
      else, end
   end

 if n==nf
 fprintf('%7g',n), fprintf('         F'),  fprintf('%12.4f', Ifm)
 fprintf('%12.4f\n', Ifmang)
 else, end
end
   resp=0;
   while strcmp(resp, 'n')~=1 & strcmp(resp, 'N')~=1 & strcmp(resp, 'y')~=1 & strcmp(resp, 'Y')~=1
   resp = 1;
   if strcmp(resp, 'n')~=1 & strcmp(resp, 'N')~=1 & strcmp(resp, 'y')~=1 & strcmp(resp, 'Y')~=1
   fprintf('\n Incorrect reply, try again \n\n'), end
   end
   if resp == 'y' | resp == 'Y'
   nf = 999;
   else ff = 0; end
end   % end for while

Answer 1

[DGM]

If neither nl(I) nor nr(I) equal n, then k is undefined.

if (nl(I) == n)
    k = nr(I);
elseif (nr(I) == n)
    k = nl(I);
% else
    % k = some default value
end

I'm guessing that's what it's complaining about.

[Himanshu Kaushish]

Thank you very much