clear all;
close all;

fid=fopen('c:\test.txt');
c=textscan(fid,'%d%d%c%c%f%f');
node=[c{1} c{2}];
var_type=[c{3}];
dat_type=[c{4}];
val=[c{5}+j*c{6}];
num_entry=size(var_type);
num_node=max(max(node));
Ybus=zeros(num_node'num_node);
I=zeros(num_node,1);
val
var_type

for k=1:num_entry
    if(var_type(k)=='I')
        if(node(k,1)~=0)
            I(node(k,1))=I(node(k,1)-val(k));
        end
        if(node(k,2)~=0)
             I(node(k,2))=I(node(k,2)+val(k));
        end
    end
    if(var_type(k)=='Y')
        if(node(k,1)~=0)
            Ybus(node(k,1),node(k,1))==Ybus(node(k,1),node(k,1))+val(k);
              if(node(k,2)~=0)
                  Ybus(node(k,1),(node(k,2))==-val(k);
              end
        end
         if(node(k,2)~=0)
             Ybus(node(k,2),node(k,2))==Ybus(node(k,2),node(k,2))+val(k);
                 if(node(k,1)~=0)
                   Ybus(node(k,2),(node(k,1))=-val(k);
              end
         end
    end
end
I
Zbus= inv(Ybus);
v=Zbus*I;
V