clear all;
close all;

fd = fopen('d:\sumon\sp.txt');
c = textscan(fd, '%d %d %c %c %f %f');
node = [c{1} c{2}];
var_type = [c{3}];
dat_type = [c{4}];
val_first = [c(5)];
val_second = [c(6)];

num_entry = length(var_type);
num_node = max(max(node));

for k = 1:num_entry
    if(dat_type(k)==';p')
        val(k) = polar2cartesian(val_first(k),val_second(k));
    else
        val(k) = val_first(k)+j*val_second(k);
    end
end

Ybus = zeros(num_node,num_node);
I = zeros(num_node, 1);

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