clear all;
close all;
fid = fopen('c:\test.txt');
c = textscan(fid, '%d %d %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)
            if(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,2)~=0)
                    Ybus(node(k,2),node(k,1))=-val(k);
                end
            end
        end
    end
    I