fast decoupled power flow method

[shoo]

I made Octave code for fast decoupled power flow method. could you find its
problem?
% Fast Decoupled Method:
V1= 1.02; V2 = 1.0; V3 = 1.0;
d1 = 0; d2 = 0; d3=0;
Ps2= 0.5;  Ps3 = -1.2;
Qs3= -0.5;
YB = [1.2113 - j*13.094  -0.2212 + j*3.3186   -0.9901 + j*9.9010
     -0.3309 + j*3.3096    1.2113 - j*13.0946  -0.9901 + j*9.9010
     -0.9901 + j*9.9010   -0.9901 + j*9.9010    1.9802 - j*19.7019];
Y = abs(YB); t=angle(YB);
B = [-13.3333   10.0000
      10.0000  -20.0000];
Binv = inv(B)
iter=0;
pwracur = 0.0001;    % Power accuracy
DC = 10;              % Power mismatch limit
while max(abs(DC)) > pwracur
iter = iter +1;
P2= V2*V1*Y(2,1)*cos(t(2,1)-d2+d1)+V2^2*Y(2,2)*cos(t(2,2))+ ...
    V2*V3*Y(2,3)*cos(t(2,3)-d2+d3);
P3= V3*V1*Y(3,1)*cos(t(3,1)-d3+d1)+V3^2*Y(3,3)*cos(t(3,3))+ ...
    V3*V2*Y(3,2)*cos(t(3,2)-d3+d2);
Q3=-V3*V1*Y(3,1)*sin(t(3,1)-d3+d1)-V3^2*Y(3,3)*sin(t(3,3))- ...
    V3*V2*Y(3,2)*sin(t(3,2)-d3+d2);
DP2 = Ps2 - P2; DP2V = DP2/V2;
DP3 = Ps3 - P3; DP3V = DP3/V3;
DQ3 = Qs3 - Q3; DQ3V = DQ3/V3;
DC =[DP2; DP3; DQ3];
Dd = -Binv*[DP2V;DP3V];
DV = -1/B(1,1)*DQ3V;
d2 =d2+Dd(1);
d3 =d3+Dd(2);
V3= V3+DV;
angle2 =180/pi*d2;
angle3 =180/pi*d3;
R = [iter  d2  d3  V2  DP2  DP3 DQ3];
disp(R)
end
Q2=-V2*V1*Y(2,1)*sin(t(2,1)-d2+d1)-V2^2*Y(2,2)*sin(t(2,2))- ...
    V2*V3*Y(2,3)*sin(t(2,3)-d2+d3);
P1= V1^2*Y(1,1)*cos(t(1,1))+V1*V2*Y(1,2)*cos(t(1,2)-d1+d2)+ ...
    V1*V3*Y(1,3)*cos(t(1,3)-d1+d3);
Q1=-V1^2*Y(1,1)*sin(t(1,1))-V1*V2*Y(1,2)*sin(t(1,2)-d1+d2)- ...
    V1*V3*Y(1,3)*sin(t(1,3)-d1+d3);
S1=P1+j*Q1
Q2




--
View this message in context: 
http://octave.1599824.n4.nabble.com/fast-decoupled-power-flow-method-tp4661645.html
Sent from the Octave - General mailing list archive at Nabble.com.