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