Scribd
Upload
77 views36 pages

Codes

This document contains data and calculations for performing an optimal load flow analysis of the IEEE 30 bus test system using the fast decoupled method. It includes bus data, line data, calculations to form the bus admittance matrix and loss coefficient matrices, and iterations to solve for bus voltages and angles within a specified tolerance. The results include changes in real and reactive power at each bus.

Uploaded by

Pravat Satpathy
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
77 views36 pages

Codes

This document contains data and calculations for performing an optimal load flow analysis of the IEEE 30 bus test system using the fast decoupled method. It includes bus data, line data, calculations to form the bus admittance matrix and loss coefficient matrices, and iterations to solve for bus voltages and angles within a specified tolerance. The results include changes in real and reactive power at each bus.

Uploaded by

Pravat Satpathy
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 36

clc

fprintf('*****************************************************************');
% OPTIMAL LOAD FLOW OF IEEE 30 BUS SYSTEM USING FAST DECOUPLED METHOD
% Busdata
%1=slack bus :: 2=generating bus :: 0=load bus
% Bus Bus Voltage Angle --Load-- ---Generator--- Injected
% No code Mag. Degree MW Mvar MW Mvar Qmin Qmax Mvar
busdata =[1 1 1.06 0 30.38 17.78 40 -40 0 0 0
2 2 1.045 0 0 0 232 0 -40 50 0
3 2 1.01 0 131.88 26.6 0 0 0 40 0
4 0 1 0 66.92 0 0 0 0 0 0
5 0 1 0 10.64 2.24 0 0 0 0 0
6 2 1.07 0 15.68 10.5 0 0 -6 24 0
7 0 1 0 0 0 0 0 0 0 0
8 2 1.09 0 0 0 0 0 -6 24 0
9 0 1 0 41.3 23.24 0 0 0 0 0
10 0 1 0 12.6 8.12 0 0 0 0 0
11 0 1 0 4.9 2.52 0 0 0 0 0
12 0 1 0 8.54 2.24 0 0 0 0 0
13 0 1 0 18.9 8.12 0 0 0 0 0
14 0 1 0 20.86 7 0 0 0 0 0];
% Formation of bus Admittance matrix
% Linedata
% Bus Bus R X 1/2 B 1 for Line Code or
% nl nr pu pu pu tap setting value
linedata=[1 2 0.01938 0.05917 0.0264 1
2 3 0.04699 0.19797 0.0219 1
2 4 0.05811 0.17632 0.0187 1
1 5 0.05403 0.22304 0.0246 1
2 5 0.05695 0.17388 0.017 1
3 4 0.06701 0.17103 0.0173 1
4 5 0.01335 0.04211 0.0064 1
5 6 0 0.25202 0 0.932
4 7 0 0.20912 0 0.978
7 8 0 0.17615 0 1
4 9 0 0.55618 0 0.969
7 9 0 0.11001 0 1
9 10 0.03181 0.0845 0 1
6 11 0.09498 0.1989 0 1
6 12 0.12291 0.25581 0 1
6 13 0.06615 0.13027 0 1
9 14 0.12711 0.27038 0 1
10 11 0.08205 0.19207 0 1
12 13 0.22092 0.19988 0 1
13 12 0.17093 0.34802 0 1];
j=sqrt(-1);
% i = sqrt(-1);
nl = linedata(:,1); nr = linedata(:,2); R = linedata(:,3);
X = linedata(:,4); Bc = j*linedata(:,5); a = linedata(:, 6);
nbr=length(linedata(:,1)); nbus = max(max(nl), max(nr));
Z = R + j*X; y= ones(nbr,1)./Z; %branch admittance
Ybus = zeros(nbus,nbus);
for k = 1:nbr % formation of off diagonal elements
if nl(k) > 0 && nr(k) > 0
Ybus(nl(k),nr(k))=Ybus(nl(k),nr(k))-y(k)/a(k);
Ybus(nr(k),nl(k))=Ybus(nl(k),nr(k));
end
end
% formation of the diagonal elements
for n=1:nbus
for k=1:nbr
if nl(k)==n
Ybus(n,n) = Ybus(n,n)+y(k)/(a(k)^2) + Bc(k);
elseif nr(k)==n
Ybus(n,n) = Ybus(n,n)+y(k)/a(k) +Bc(k);
end
end
end
fprintf(' bus admittance matix ');
display(Ybus); % Displaying the YBUS
ns=0; yload=0; deltad=[]; %Vm=0;
basemva = 100;
accuracy = 0.001; maxiter = 100;
nbus = length(busdata(:,1));
kb=[];Vm=[]; delta=[]; Pd=[]; Qd=[]; Pg=[]; Qg=[]; Qmin=[]; Qmax=[];
Pk=[]; P=[]; Qk=[]; Q=[]; S=[]; V=[];
for k=1:nbus
n=busdata(k,1);
kb(n)=busdata(k,2); Vm(n)=busdata(k,3); delta(n)=busdata(k,4);
Pd(n)=busdata(k,5); Qd(n)=busdata(k,6); Pg(n)=busdata(k,7); Qg(n) =
busdata(k,8);
Qmin(n)=busdata(k,9); Qmax(n)=busdata(k,10);
Qsh(n)=busdata(k,11);
if Vm(n) <= 0 Vm(n) = 1.0; V(n) = 1 + j*0;
else delta(n) = pi/180*delta(n);
V(n) = Vm(n)*(cos(delta(n)) + j*sin(delta(n)));
P(n)=(Pg(n)-Pd(n))/basemva;
Q(n)=(Qg(n)-Qd(n)+Qsh(n))/basemva;
S(n) = P(n) + j*Q(n);
end
if kb(n) == 1,
ns = ns+1;
nss(n) = ns;
end
end
Ym = abs(Ybus); t = angle(Ybus);
ii=0;
%formation of LOSS COEFICIENT MATRIX
for ib=1:nbus
if kb(ib) == 0 || kb(ib) == 2
ii = ii+1;
jj=0;
for jb=1:nbus
if kb(jb) == 0 || kb(jb) == 2
jj = jj+1;
B2(ii,jj)=imag(Ybus(ib,jb));
B1(ii,jj)=busdata(jb,3)*B2(ii,jj);
end

end
end
end
for ib=1:nbus
if kb(ib) == 1
ii=1;
jj=1;
for jb=2:nbus
B22(ii,jj)=imag(Ybus(ib,jb));
B11(ii,jj)=B22(ii,jj)*busdata(jb,3);
jj=jj+1;
end
end
end
fprintf('*****************************************************************');
display(B1);
fprintf('*****************************************************************');
display(B2);
X1=inv(B1);
X2=inv(B2);
maxerror = 1; converge = 1;
%%%% added for parallel lines
mline=ones(nbr,1);
for k=1:nbr
for m=k+1:nbr
if((nl(k)==nl(m)) && (nr(k)==nr(m)));
mline(m)=2;
elseif ((nl(k)==nr(m)) && (nr(k)==nl(m)));
mline(m)=2;
end
end
end
%%% end of statements for parallel lines
% Start of iterations
iter=0;
while maxerror > accuracy && iter <= maxiter % Test for max. power mismatch
iter = iter+1;
id=0; iv=0;
for n=1:nbus
% nn=n-nss(n);
J11=0; J33=0;
for ii=1:nbr
if mline(ii)==1 % Added to include parallel lines
if nl(ii) == n || nr(ii) == n
if nl(ii) == n, l = nr(ii);end
if nr(ii) == n, l = nl(ii);end
j11=J11+ Vm(n)*Vm(l)*Ym(n,l)*sin(t(n,l)- delta(n) +
delta(l));
J33=J33+ Vm(n)*Vm(l)*Ym(n,l)*cos(t(n,l)- delta(n) +
delta(l));
end
end
end
Pk = Vm(n)^2*Ym(n,n)*cos(t(n,n))+J33;
Qk = -Vm(n)^2*Ym(n,n)*sin(t(n,n))-J11;
if kb(n) == 1
P(n)=Pk; % Slack bus P
Q(n) = Qk;
end
if kb(n) == 2
Q(n)=Qk;
Qgc = Q(n)*basemva - Qd(n) + Qsh(n);
if Qmax(n) ~= 0
if iter <= 20 % Between the 1th & 6th iterations
if iter >= 10 % the Mvar of generator buses are
if Qgc < Qmin(n) % tested. If not within limits Vm(n)
Vm(n) = Vm(n) + 0.005;% is changed in steps of 0.05 pu to
% Qgc = Qmin(n);
elseif Qgc > Qmax(n) % bring the generator Mvar within
Vm(n) = Vm(n) - 0.005; % the specified limits.
% Qgc = Qmax(n);
end
end
end
end
end
if kb(n) ~= 1
id = id+1;
DP(id) = P(n)-Pk; %change in real power
DPV(id) = (P(n)-Pk)/Vm(n);
end
if kb(n) == 0 || kb(n)==2
iv=iv+1;
DQ(iv) = Q(n)-Qk; %change in reactive power
DQV(iv) = (Q(n)-Qk)/Vm(n);
end
end
Dd=-X1*(DPV)';
DV=-X2*(DQV)';
id=0;iv=0;
maxerror=max(max(abs(DP)),max(abs(DQ)));
end
fprintf('*****************************************************************');
DP=DPV.*busdata(2:nbus,3)';
DQ=DQV.*busdata(2:nbus,3)';

fprintf('del P(except slack bus= ');


display(DP);
for s1=1:length(B11)
yi(s1)=0;
yii(s1)=0;
for s=1:length(B11)
yi(s1)=yi(s1)+B11(1,s)*X1(s1,s);
yii(s1)=yii(s1)+B22(1,s)*X2(s1,s);
end

end
fprintf('*****************************************************************');
display(yi);
fprintf('*****************************************************************');
display(yii);

DP1=B11*X1*DP';
DQ1=B22*X2*DQ';
for s=1:length(B11)
jj=0;
for s1=1:nbus
if s1 == 1
f1(s,s1)=1/yi(s)*B1(s,s);
f2(s,s1)=1/yii(s)*B2(s,s);
elseif (s+1) == s1
m=0;m1=0;
for s2=1:length(B11)
if (s2 < s) || (s2 > s)
m=m+B1(s,s2)*X1(s2,s);
m1=m1+B2(s,s2)*X2(s2,s);
end
end
f1(s,s1)=-m/B1(s,s);
f2(s,s1)=-m1/B2(s,s);
else
m=0;m1=0;
for s2=1:length(B11)
if ((s2 < s) || (s2 > s))
m=m+B1(s,s2)*X1(s2,s);
m1=m1+B2(s,s2)*X2(s2,s);
end
end
jj=jj+1;
f1(s,s1)=-yi(s1-1)/(yi(s)*B1(s,s))-1/B1(s,s)*m;
f2(s,s1)=-yii(s1-1)/(yii(s)*B2(s,s))-1/B2(s,s)*m1;
end
end
end
fprintf('*****************************************************************');
display(f1);
fprintf('*****************************************************************');
display(f2);

vk=busdata(:,3);
dk=busdata(:,4);

for k = 1:nbr
if nl(k) > 0 && nr(k) > 0
g(k)=1/(R(k)*a(k));
b(k)=1/(X(k)*a(k));
Ir(k)=(vk(nl(k))-vk(nr(k)))*g(k);
Ij(k)=-(vk(nl(k))-vk(nr(k)))*b(k)+Bc((nl(k)))*vk((nl(k)));
I(k)=Ir(k)+j*Ij(k);
DIDK(k)=Ir(k)/I(k)*(-b(k)*vk(nl(k))-
Bc((nl(k)))*vk((nl(k)))+Ij(k)/I(k)*g(k)*vk(nl(k)));
DIDVK(k)=Ir(k)/I(k)*g(k)+Ij(k)/I(k)*-b(k)-Bc((nl(k)));
DIDM(k)=Ir(k)/I(k)*b(k)*vk(nr(k))+Ij(k)/I(k)*-g(k)*vk(nr(k));
DIDVM(k)=Ir(k)/I(k)*-g(k)+Ij(k)/I(k)*-b(k)*cos(dk(nr(k)));
end
end

for k=1:nbr
for ii=1:length(B11)
if nl(k)>1
f11(k,ii)=DIDK(k)*(f1(nl(k)-1,1)*yi(ii)+f1(nl(k)-
1,ii+1))+DIDM(k)*(f1(nr(k)-1,1)*yi(ii)+f1(nr(k)-1,ii+1));
f22(k,ii)=DIDVK(k)*(f2(nl(k)-1,1)*yii(ii)+f2(nl(k)-
1,ii+1))+DIDVM(k)*(f2(nr(k)-1,1)*yii(ii)+f2(nr(k)-1,ii+1));
else
f11(k,ii)=-(DIDK(k)*(f1(nr(k)-1,1)*yi(ii)+f1(nr(k)-
1,ii+1))+DIDM(k)*(f1(nl(k),1)*yi(ii)+f1(nl(k),ii+1)));
f22(k,ii)=-
(DIDVK(k)*(f2(nr(k),1)*yii(ii)+f2(nr(k),ii+1))+DIDVM(k)*(f2(nl(k),1)*yii(ii)+
f2(nl(k),ii+1)));
end
end
end
f11
f22
P=busdata(2:nbus,5);
Q=busdata(2:nbus,6);
beta=Q./P;
vkk=vk(2:nbus);
for k=1:nbr
for k1=1:13
F(k,k1)=(f11(k,k1)+beta(k1)*f22(k,k1))/vkk(k1);
end
end
fprintf('*****************************************************************');
fprintf('senaitivity matix ');
display(F);
DI=F*DP';
fprintf('*****************************************************************');
fprintf('change in current ');
display(DI);

%=====================================================
%selection of paticipating generators
for i=1:4
C(i)=input('enter incremental fuel cost(in Rs/MW) of generators: ');
end
%calculation for f_up_cut and f_lo_cut for generator1
%buscode=busdata(:,2);
for i=1:20
for j=1:13
if F(i,j)>0
F1(i,j)=C(1)*F(i,j);
else if F(i,j)<0
F1(i,j)=F(i,j)/C(1);
else F1(i,j)=0;
end
end
end
end
MF_up1=input('multiplying factor for increase generation of 1: ');
MF_lo1=input('multiplying factor for decrease generation of 1: ');
for i=1:20
for j=1:13
if F1(i,j)>0
f_up_cut=max(F1(i,j));
else if F1(i,j)<0
f_lo_cut=min(F1(i,j));
end
end
end
end
f_lo_cut1=MF_lo1*f_lo_cut
f_up_cut1=MF_up1*f_up_cut
%calculation for f_up_cut and f_lo_cut for generator2
for i=1:20
for j=1:13
if F(i,j)>0
F2(i,j)=C(2)*F(i,j);
else if F(i,j)<0
F2(i,j)=F(i,j)/C(2);
else
F2(i,j)=0;
end
end
end
end
MF_up2=input('multiplying factor for increase generation of 2: ');
MF_lo2=input('multiplying factor for decrease generation of 2: ');
for i=1:20
for j=1:13
if F2(i,j)>0
f_up_cut=max(F2(i,j));
else if F2(i,j)<0
f_lo_cut=min(F2(i,j));
end
end
end
end
f_lo_cut2=MF_lo2*f_lo_cut
f_up_cut2=MF_up2*f_up_cut
%calculation for f_up_cut and f_lo_cut for generator3
for i=1:20
for j=1:13
if F(i,j)>0
F3(i,j)=C(3)*F(i,j);
else if F(i,j)<0
F3(i,j)=F(i,j)/C(3);
else F3(i,j)=0;
end
end
end
end
MF_up3=input('multiplying factor for increase generation of 3: ');
MF_lo3=input('multiplying factor for decrease generation of 3: ');
for i=1:20
for j=1:13
if F3(i,j)>0
f_up_cut=max(F3(i,j));
else if F3(i,j)<0
f_lo_cut=min(F3(i,j));
end
end
end
end
f_lo_cut3=MF_lo3*f_lo_cut
f_up_cut3=MF_up3*f_up_cut
%calculation for f_up_cut and f_lo_cut for generator4
for i=1:20
for j=1:13
if F(i,j)>0
F4(i,j)=C(4)*F(i,j);
else if F(i,j)<0
F4(i,j)=F(i,j)/C(4);
else F4(i,j)=0;
end
end
end
end
MF_up4=input('multiplying factor for increase generation of 4: ');
MF_lo4=input('multiplying factor for decrease generation of 4: ');
for i=1:20
for j=1:13
if F4(i,j)>0
f_up_cut=max(F4(i,j));
else if F4(i,j)<0
f_lo_cut=min(F4(i,j));
end
end
end
end
f_lo_cut4=MF_lo4*f_lo_cut
f_up_cut4=MF_up4*f_up_cut

%participating loads for load shedding


for i=1:9
C1(i)=input('enter incremental cost(in Rs/MW) of load shedding: ');
end
F1L=F*C1(1);
F2L=F*C1(2);
F3L=F*C1(3);
F4L=F*C1(4);
F5L=F*C1(5);
F6L=F*C1(2);
F7L=F*C1(7);
F8L=F*C1(8);
F9L=F*C1(9);
Fl=[min(min(F1L)) min(min(F2L)) min(min(F3L)) min(min(F4L)) min(min(F5L))
min(min(F6L)) min(min(F7L)) min(min(F8L)) min(min(F9L))];
for i=1:9
t=0;
for j=(i+1):9
if Fl(i)>Fl(j)
t=Fl(i);
Fl(i)=Fl(j);
Fl(j)=t;
end
end
end
Fl
MF_lod_shdng=input('multiplying factor for load shedding: ');
f_lo_lod_shdng=MF_lod_shdng*min(Fl)
%Simultanous generator rescheduling and load shedding

MF_up_gen=[MF_up1 MF_up2 MF_up3];


MF_lo_gen=[MF_lo1 MF_lo2 MF_lo3];
G=0;
for i=1:3
if MF_lo_gen(i)>0 || MF_up_gen(i)>0
G=G+1;
end
end
D=0;
for i=1:9
if Fl(i)<(f_lo_lod_shdng/2)
D=D+1;
end
end

for i=1:20
for j=1:(G+D)
Fgd(i,j)=F(i,j);
end
end
Ct=[C C1];
for i=1:G+D
for j=1:G+D
if(G==0) && (i==j)
W(i,j)=C1(i);
else if i==j
W(i,j)=Ct(i);
end
end
end
end
F_new=Fgd*inv(W);
for i=1:G+D
DP_new1(i)=DP(i);
end
DP_new=W*DP_new1';
if nbr<(G+D)
DP_new=F_new'*inv((F_new*F_new'))*DI;
else
DP_new=inv((F_new'*F_new))*F_new'*DI;
end
DP_new
output
Ybus =

Columns 1 through 4

6.0250 -19.4471i -4.9991 +15.2631i 0 0

-4.9991 +15.2631i 9.5213 -30.2707i -1.1350 + 4.7819i -1.6860 + 5.1158i

0 -1.1350 + 4.7819i 3.1210 - 9.8115i -1.9860 + 5.0688i

0 -1.6860 + 5.1158i -1.9860 + 5.0688i 10.5130 -38.6352i

-1.0259 + 4.2350i -1.7011 + 5.1939i 0 -6.8410 +21.5786i

0 0 0 0

0 0 0 0 + 4.8895i

0 0 0 0

0 0 0 0 + 1.8555i

0 0 0 0

0 0 0 0

0 0 0 0

0 0 0 0

0 0 0 0

Columns 5 through 8

-1.0259 + 4.2350i 0 0 0
-1.7011 + 5.1939i 0 0 0

0 0 0 0

-6.8410 +21.5786i 0 0 + 4.8895i 0

9.5680 -35.5275i 0 + 4.2574i 0 0

0 + 4.2574i 6.5799 -17.6302i 0 0

0 0 0 -19.6566i 0 + 5.6770i

0 0 0 + 5.6770i 0 - 5.6770i

0 0 0 + 9.0901i 0

0 0 0 0

0 -1.9550 + 4.0941i 0 0

0 -1.5260 + 3.1760i 0 0

0 -3.0989 + 6.1028i 0 0

0 0 0 0

Columns 9 through 12

0 0 0 0

0 0 0 0

0 0 0 0

0 + 1.8555i 0 0 0

0 0 0 0

0 0 -1.9550 + 4.0941i -1.5260 + 3.1760i

0 + 9.0901i 0 0 0

0 0 0 0

5.3261 -24.3400i -3.9020 +10.3654i 0 0


-3.9020 +10.3654i 5.7829 -14.7683i -1.8809 + 4.4029i 0

0 -1.8809 + 4.4029i 3.8359 - 8.4970i 0

0 0 0 5.1520 - 7.7429i

0 0 0 -3.6260 + 4.5669i

-1.4240 + 3.0291i 0 0 0

Columns 13 through 14

0 0

0 0

0 0

0 0

0 0

-3.0989 + 6.1028i 0

0 0

0 0

0 -1.4240 + 3.0291i

0 0

0 0

-3.6260 + 4.5669i 0

6.7249 -10.6697i 0

0 1.4240 - 3.0291i

*****************************************************************

B1 =
Columns 1 through 8

-31.6329 4.8297 5.1158 5.1939 0 0 0 0

4.9970 -9.9096 5.0688 0 0 0 0 0

5.3461 5.1195 -38.6352 21.5786 0 4.8895 0 1.8555

5.4277 0 21.5786 -35.5275 4.5555 0 0 0

0 0 0 4.2574 -18.8644 0 0 0

0 0 4.8895 0 0 -19.6566 6.1879 9.0901

0 0 0 0 0 5.6770 -6.1879 0

0 0 1.8555 0 0 9.0901 0 -24.3400

0 0 0 0 0 0 0 10.3654

0 0 0 0 4.3807 0 0 0

0 0 0 0 3.3983 0 0 0

0 0 0 0 6.5299 0 0 0

0 0 0 0 0 0 0 3.0291

Columns 9 through 13

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

0 4.0941 3.1760 6.1028 0

0 0 0 0 0
0 0 0 0 0

10.3654 0 0 0 3.0291

-14.7683 4.4029 0 0 0

4.4029 -8.4970 0 0 0

0 0 -7.7429 4.5669 0

0 0 4.5669 -10.6697 0

0 0 0 0 -3.0291

*****************************************************************

B2 =

Columns 1 through 8

-30.2707 4.7819 5.1158 5.1939 0 0 0 0

4.7819 -9.8115 5.0688 0 0 0 0 0

5.1158 5.0688 -38.6352 21.5786 0 4.8895 0 1.8555

5.1939 0 21.5786 -35.5275 4.2574 0 0 0

0 0 0 4.2574 -17.6302 0 0 0

0 0 4.8895 0 0 -19.6566 5.6770 9.0901

0 0 0 0 0 5.6770 -5.6770 0

0 0 1.8555 0 0 9.0901 0 -24.3400

0 0 0 0 0 0 0 10.3654

0 0 0 0 4.0941 0 0 0

0 0 0 0 3.1760 0 0 0

0 0 0 0 6.1028 0 0 0
0 0 0 0 0 0 0 3.0291

Columns 9 through 13

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

0 0 0 0 0

0 4.0941 3.1760 6.1028 0

0 0 0 0 0

0 0 0 0 0

10.3654 0 0 0 3.0291

-14.7683 4.4029 0 0 0

4.4029 -8.4970 0 0 0

0 0 -7.7429 4.5669 0

0 0 4.5669 -10.6697 0

0 0 0 0 -3.0291

*****************************************************************del P(except slack bus=

DP =

Columns 1 through 8

2.8085 -1.2644 -0.7754 -0.0619 -0.2709 -0.0000 -0.0000 -0.4130


Columns 9 through 13

-0.1260 -0.0197 -0.0625 -0.1425 -0.2086

*****************************************************************

yi =

Columns 1 through 8

-0.9841 -1.0087 -1.0019 -0.9978 -0.9334 -1.0014 -0.9187 -1.0011

Columns 9 through 13

-1.0007 -0.9998 -0.9988 -0.9988 -1.0011

*****************************************************************

yii =

Columns 1 through 8

-0.9911 -0.9857 -0.9730 -0.9715 -0.9719 -0.9728 -0.9728 -0.9727

Columns 9 through 13

-0.9726 -0.9722 -0.9719 -0.9719 -0.9727


*****************************************************************

f1 =

Columns 1 through 8

32.1439 -0.0204 0.0120 0.0118 0.0117 0.0096 0.0118 0.0091

9.8237 0.0378 -0.0607 0.0396 0.0392 0.0327 0.0395 0.0312

38.5604 -0.0504 -0.0498 -0.0759 -0.0501 -0.0517 -0.0500 -0.0521

35.6063 -0.0358 -0.0351 -0.0353 -0.0635 -0.0372 -0.0353 -0.0376

20.2099 -0.1407 -0.1393 -0.1397 -0.1400 -0.1966 -0.1397 -0.1444

19.6290 -0.1342 -0.1330 -0.1333 -0.1335 -0.1368 -0.1842 -0.1375

6.7354 -0.0426 -0.0382 -0.0394 -0.0402 -0.0515 -0.0395 -0.2157

24.3129 -0.1715 -0.1705 -0.1708 -0.1710 -0.1736 -0.1708 -0.1742

14.7578 -0.1648 -0.1631 -0.1636 -0.1639 -0.1682 -0.1636 -0.1692

8.4989 -0.1035 -0.1006 -0.1014 -0.1019 -0.1094 -0.1014 -0.1112

7.7525 -0.1835 -0.1803 -0.1812 -0.1817 -0.1900 -0.1812 -0.1919

10.6829 -0.2064 -0.2041 -0.2047 -0.2051 -0.2112 -0.2048 -0.2125

3.0257 0.0715 0.0796 0.0774 0.0760 0.0548 0.0772 0.0500

Columns 9 through 14

0.0118 0.0118 0.0117 0.0117 0.0117 0.0118

0.0395 0.0394 0.0394 0.0393 0.0393 0.0395

-0.0500 -0.0500 -0.0500 -0.0501 -0.0501 -0.0500


-0.0353 -0.0353 -0.0353 -0.0354 -0.0354 -0.0353

-0.1398 -0.1398 -0.1398 -0.1399 -0.1399 -0.1398

-0.1334 -0.1334 -0.1334 -0.1335 -0.1335 -0.1334

-0.0396 -0.0396 -0.0398 -0.0400 -0.0400 -0.0396

-0.2119 -0.1709 -0.1709 -0.1709 -0.1709 -0.1708

-0.1637 -0.2314 -0.1638 -0.1638 -0.1638 -0.1637

-0.1015 -0.1015 -0.2193 -0.1018 -0.1018 -0.1015

-0.1813 -0.1813 -0.1814 -0.3107 -0.1816 -0.1813

-0.2048 -0.2048 -0.2049 -0.2050 -0.2988 -0.2048

0.0771 0.0770 0.0767 0.0764 0.0764 -0.2530

*****************************************************************

f2 =

Columns 1 through 8

30.5434 -0.0213 0.0116 0.0111 0.0111 0.0111 0.0111 0.0111

9.9538 0.0412 -0.0613 0.0393 0.0392 0.0392 0.0393 0.0393

39.7068 -0.0495 -0.0496 -0.0759 -0.0500 -0.0500 -0.0500 -0.0500

36.5682 -0.0348 -0.0350 -0.0354 -0.0635 -0.0354 -0.0354 -0.0354

18.1402 -0.1525 -0.1529 -0.1536 -0.1537 -0.2104 -0.1536 -0.1536

20.2057 -0.1324 -0.1327 -0.1333 -0.1334 -0.1334 -0.1842 -0.1333

5.8356 -0.0556 -0.0566 -0.0589 -0.0592 -0.0591 -0.0589 -0.2351

25.0226 -0.1701 -0.1703 -0.1708 -0.1709 -0.1709 -0.1708 -0.1708

15.1847 -0.1624 -0.1628 -0.1637 -0.1638 -0.1637 -0.1637 -0.1637


8.7396 -0.0994 -0.1000 -0.1015 -0.1017 -0.1017 -0.1016 -0.1016

7.9669 -0.1790 -0.1797 -0.1814 -0.1816 -0.1816 -0.1814 -0.1814

10.9783 -0.2032 -0.2037 -0.2049 -0.2051 -0.2050 -0.2049 -0.2049

3.1140 0.0834 0.0815 0.0772 0.0767 0.0768 0.0772 0.0772

Columns 9 through 14

0.0111 0.0111 0.0111 0.0111 0.0111 0.0111

0.0393 0.0393 0.0393 0.0392 0.0392 0.0393

-0.0500 -0.0500 -0.0500 -0.0500 -0.0500 -0.0500

-0.0354 -0.0354 -0.0354 -0.0354 -0.0354 -0.0354

-0.1536 -0.1536 -0.1536 -0.1537 -0.1537 -0.1536

-0.1333 -0.1333 -0.1334 -0.1334 -0.1334 -0.1333

-0.0590 -0.0590 -0.0590 -0.0591 -0.0591 -0.0590

-0.2119 -0.1708 -0.1709 -0.1709 -0.1709 -0.1708

-0.1637 -0.2314 -0.1637 -0.1637 -0.1637 -0.1637

-0.1016 -0.1016 -0.2193 -0.1017 -0.1017 -0.1016

-0.1815 -0.1815 -0.1815 -0.3107 -0.1816 -0.1815

-0.2049 -0.2050 -0.2050 -0.2050 -0.2988 -0.2049

0.0771 0.0771 0.0770 0.0768 0.0768 -0.2530

f11 =

1.0e+003 *
Columns 1 through 4

0.0616 - 0.1696i 0.0630 - 0.1736i 0.0626 - 0.1725i 0.0624 - 0.1718i

0.2241 + 0.0816i 0.2284 + 0.0833i 0.2279 + 0.0829i 0.2269 + 0.0826i

-0.0798 + 0.0382i -0.0821 + 0.0390i -0.0818 + 0.0387i -0.0812 + 0.0386i

-0.0371 - 0.0492i -0.0383 - 0.0505i -0.0380 - 0.0502i -0.0381 - 0.0501i

-0.0491 + 0.0476i -0.0507 + 0.0487i -0.0503 + 0.0483i -0.0504 + 0.0481i

-0.3438 - 0.0641i -0.3513 - 0.0651i -0.3504 - 0.0653i -0.3486 - 0.0650i

0 - 2.6274i 0 - 2.6931i 0 - 2.6750i 0 - 2.6660i

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN NaN NaN NaN

0.0775 + 0.0806i 0.0794 + 0.0826i 0.0789 + 0.0821i 0.0786 + 0.0817i

0.0651 + 0.0656i 0.0667 + 0.0673i 0.0663 + 0.0668i 0.0660 + 0.0665i

0.0802 + 0.1096i 0.0822 + 0.1123i 0.0817 + 0.1116i 0.0813 + 0.1111i

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

Columns 5 through 8
0.0583 - 0.1607i 0.0626 - 0.1724i 0.0574 - 0.1582i 0.0626 - 0.1723i

0.2123 + 0.0773i 0.2278 + 0.0829i 0.2089 + 0.0760i 0.2277 + 0.0829i

-0.0760 + 0.0361i -0.0815 + 0.0387i -0.0748 + 0.0355i -0.0815 + 0.0387i

-0.0355 - 0.0468i -0.0380 - 0.0502i -0.0349 - 0.0460i -0.0380 - 0.0501i

-0.0469 + 0.0450i -0.0503 + 0.0483i -0.0462 + 0.0443i -0.0503 + 0.0483i

-0.3261 - 0.0608i -0.3499 - 0.0652i -0.3210 - 0.0598i -0.3498 - 0.0652i

0 - 2.4924i 0 - 2.6735i 0 - 2.4532i 0 - 2.6728i

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN NaN NaN NaN

0.0739 + 0.0768i 0.0789 + 0.0820i 0.0723 + 0.0753i 0.0788 + 0.0820i

0.0620 + 0.0625i 0.0662 + 0.0668i 0.0607 + 0.0613i 0.0662 + 0.0668i

0.0765 + 0.1044i 0.0816 + 0.1115i 0.0748 + 0.1023i 0.0816 + 0.1115i

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

Columns 9 through 12

0.0625 - 0.1723i 0.0625 - 0.1721i 0.0624 - 0.1719i 0.0624 - 0.1719i


0.2276 + 0.0828i 0.2274 + 0.0828i 0.2272 + 0.0827i 0.2272 + 0.0827i

-0.0814 + 0.0387i -0.0813 + 0.0387i -0.0813 + 0.0386i -0.0813 + 0.0386i

-0.0380 - 0.0501i -0.0380 - 0.0501i -0.0379 - 0.0500i -0.0379 - 0.0500i

-0.0503 + 0.0483i -0.0502 + 0.0482i -0.0502 + 0.0482i -0.0502 + 0.0482i

-0.3496 - 0.0652i -0.3493 - 0.0651i -0.3490 - 0.0650i -0.3490 - 0.0650i

0 - 2.6717i 0 - 2.6692i 0 - 2.6665i 0 - 2.6665i

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN NaN NaN NaN

0.0788 + 0.0820i 0.0777 + 0.0814i 0.0787 + 0.0818i 0.0787 + 0.0818i

0.0662 + 0.0667i 0.0661 + 0.0667i 0.0652 + 0.0662i 0.0661 + 0.0666i

0.0816 + 0.1114i 0.0815 + 0.1113i 0.0814 + 0.1112i 0.0802 + 0.1106i

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

Column 13

0.0626 - 0.1723i

0.2277 + 0.0829i

-0.0815 + 0.0387i
-0.0380 - 0.0501i

-0.0503 + 0.0483i

-0.3498 - 0.0652i

0 - 2.6728i

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN

0.0788 + 0.0820i

0.0662 + 0.0668i

0.0816 + 0.1115i

NaN

NaN

NaN

NaN

f22 =

1.0e+003 *

Columns 1 through 4
-0.7681 - 0.2867i -0.7570 - 0.2829i -0.7526 - 0.2810i -0.7515 - 0.2805i

-0.4723 - 0.1086i -0.4670 - 0.1073i -0.4629 - 0.1064i -0.4622 - 0.1063i

0.0208 + 0.0198i 0.0213 + 0.0199i 0.0214 + 0.0197i 0.0210 + 0.0196i

-0.1653 - 0.0472i -0.1639 - 0.0468i -0.1617 - 0.0462i -0.1615 - 0.0461i

-0.0222 + 0.0046i -0.0215 + 0.0048i -0.0212 + 0.0047i -0.0208 + 0.0048i

0.2976 + 0.1754i 0.2943 + 0.1738i 0.2925 + 0.1724i 0.2918 + 0.1720i

0 - 1.7962i 0 - 1.7865i 0 - 1.7641i 0 - 1.7615i

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN NaN NaN NaN

-0.1364 - 0.0630i -0.1357 - 0.0627i -0.1339 - 0.0619i -0.1337 - 0.0618i

-0.1100 - 0.0509i -0.1094 - 0.0506i -0.1080 - 0.0500i -0.1079 - 0.0499i

-0.1787 - 0.0880i -0.1777 - 0.0876i -0.1754 - 0.0864i -0.1752 - 0.0863i

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

Columns 5 through 8

-0.7518 - 0.2806i -0.7525 - 0.2809i -0.7525 - 0.2809i -0.7524 - 0.2809i

-0.4624 - 0.1063i -0.4628 - 0.1064i -0.4628 - 0.1064i -0.4628 - 0.1064i


0.0210 + 0.0196i 0.0210 + 0.0196i 0.0210 + 0.0196i 0.0210 + 0.0196i

-0.1605 - 0.0459i -0.1617 - 0.0462i -0.1617 - 0.0462i -0.1617 - 0.0462i

-0.0212 + 0.0047i -0.0212 + 0.0047i -0.0212 + 0.0047i -0.0212 + 0.0047i

0.2919 + 0.1720i 0.2922 + 0.1722i 0.2922 + 0.1722i 0.2921 + 0.1722i

0 - 1.7615i 0 - 1.7632i 0 - 1.7632i 0 - 1.7630i

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN NaN NaN NaN

-0.1344 - 0.0621i -0.1339 - 0.0619i -0.1339 - 0.0619i -0.1339 - 0.0619i

-0.1084 - 0.0501i -0.1080 - 0.0499i -0.1080 - 0.0499i -0.1080 - 0.0499i

-0.1761 - 0.0868i -0.1754 - 0.0864i -0.1754 - 0.0864i -0.1754 - 0.0864i

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

Columns 9 through 12

-0.7523 - 0.2808i -0.7520 - 0.2807i -0.7518 - 0.2806i -0.7518 - 0.2806i

-0.4627 - 0.1064i -0.4626 - 0.1063i -0.4624 - 0.1063i -0.4624 - 0.1063i

0.0210 + 0.0196i 0.0210 + 0.0196i 0.0210 + 0.0196i 0.0210 + 0.0196i

-0.1616 - 0.0461i -0.1616 - 0.0461i -0.1615 - 0.0461i -0.1615 - 0.0461i


-0.0212 + 0.0047i -0.0212 + 0.0047i -0.0212 + 0.0047i -0.0212 + 0.0047i

0.2921 + 0.1721i 0.2920 + 0.1721i 0.2919 + 0.1720i 0.2919 + 0.1720i

0 - 1.7627i 0 - 1.7621i 0 - 1.7615i 0 - 1.7615i

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN NaN NaN NaN

-0.1339 - 0.0618i -0.1330 - 0.0614i -0.1338 - 0.0618i -0.1338 - 0.0618i

-0.1080 - 0.0499i -0.1079 - 0.0499i -0.1072 - 0.0495i -0.1079 - 0.0499i

-0.1753 - 0.0864i -0.1753 - 0.0864i -0.1752 - 0.0863i -0.1744 - 0.0859i

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

Column 13

-0.7524 - 0.2809i

-0.4628 - 0.1064i

0.0210 + 0.0196i

-0.1617 - 0.0462i

-0.0212 + 0.0047i

0.2921 + 0.1722i
0 - 1.7630i

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN

-0.1339 - 0.0619i

-0.1080 - 0.0499i

-0.1754 - 0.0864i

NaN

NaN

NaN

NaN

*****************************************************************senaitivity matix

F=

1.0e+003 *

Columns 1 through 4

NaN + NaNi -0.0888 - 0.2284i 0.0626 - 0.1725i -0.0959 - 0.2308i

NaN + NaNi 0.1329 + 0.0611i 0.2279 + 0.0829i 0.1296 + 0.0602i

NaN + NaNi -0.0770 + 0.0426i -0.0818 + 0.0387i -0.0768 + 0.0427i


NaN + NaNi -0.0706 - 0.0594i -0.0380 - 0.0502i -0.0721 - 0.0598i

NaN + NaNi -0.0545 + 0.0491i -0.0503 + 0.0483i -0.0548 + 0.0491i

NaN + NaNi -0.2890 - 0.0297i -0.3504 - 0.0653i -0.2872 - 0.0288i

NaN + NaNi 0 - 3.0232i 0 - 2.6750i 0 - 3.0368i

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN NaN NaN NaN

NaN + NaNi 0.0516 + 0.0693i 0.0789 + 0.0821i 0.0504 + 0.0687i

NaN + NaNi 0.0442 + 0.0565i 0.0663 + 0.0668i 0.0433 + 0.0560i

NaN + NaNi 0.0459 + 0.0937i 0.0817 + 0.1116i 0.0445 + 0.0929i

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

Columns 5 through 8

-0.4160 - 0.3258i NaN + NaNi NaN + NaNi -0.3608 - 0.3304i

-0.0910 + 0.0057i NaN + NaNi NaN + NaNi -0.0327 + 0.0230i

-0.0579 + 0.0460i NaN + NaNi NaN + NaNi -0.0696 + 0.0497i

-0.1336 - 0.0724i NaN + NaNi NaN + NaNi -0.1290 - 0.0761i

-0.0571 + 0.0450i NaN + NaNi NaN + NaNi -0.0622 + 0.0509i


-0.1221 + 0.0509i NaN + NaNi NaN + NaNi -0.1854 + 0.0317i

0 - 3.4317i NaN + NaNi NaN + NaNi 0 - 3.6648i

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN NaN NaN NaN

-0.0151 + 0.0329i NaN + NaNi NaN + NaNi 0.0035 + 0.0472i

-0.0099 + 0.0270i NaN + NaNi NaN + NaNi 0.0055 + 0.0387i

-0.0387 + 0.0433i NaN + NaNi NaN + NaNi -0.0171 + 0.0629i

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

Columns 9 through 12

-0.4223 - 0.3532i -0.3243 - 0.3165i -0.1348 - 0.2455i -0.2606 - 0.2925i

-0.0706 + 0.0143i -0.0105 + 0.0281i 0.1059 + 0.0548i 0.0285 + 0.0370i

-0.0679 + 0.0513i -0.0705 + 0.0487i -0.0758 + 0.0438i -0.0722 + 0.0470i

-0.1422 - 0.0799i -0.1211 - 0.0738i -0.0803 - 0.0621i -0.1073 - 0.0698i

-0.0639 + 0.0513i -0.0611 + 0.0506i -0.0557 + 0.0494i -0.0593 + 0.0502i

-0.1614 + 0.0458i -0.1991 + 0.0234i -0.2724 - 0.0199i -0.2235 + 0.0089i

0 - 3.8077i 0 - 3.5754i 0 - 3.1285i 0 - 3.4233i


NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN + NaNi NaN + NaNi NaN + NaNi NaN + NaNi

NaN NaN NaN NaN

-0.0075 + 0.0421i 0.0093 + 0.0499i 0.0436 + 0.0656i 0.0212 + 0.0553i

-0.0034 + 0.0346i 0.0106 + 0.0410i 0.0371 + 0.0532i 0.0197 + 0.0452i

-0.0314 + 0.0558i -0.0087 + 0.0669i 0.0355 + 0.0886i 0.0053 + 0.0737i

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

NaN NaN NaN NaN

Column 13

-0.1899 - 0.2666i

0.0724 + 0.0472i

-0.0744 + 0.0453i

-0.0923 - 0.0656i

-0.0574 + 0.0499i

-0.2517 - 0.0074i

0 - 3.2644i

NaN + NaNi

NaN + NaNi
NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN

0.0339 + 0.0612i

0.0300 + 0.0500i

0.0228 + 0.0825i

NaN

NaN

NaN

NaN

*****************************************************************change in current

DI =

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN + NaNi
NaN + NaNi

NaN + NaNi

NaN

NaN + NaNi

NaN + NaNi

NaN + NaNi

NaN

NaN

NaN

NaN

enter incremental fuel cost(in Rs/MW) of generators: 10

enter incremental fuel cost(in Rs/MW) of generators: 10

enter incremental fuel cost(in Rs/MW) of generators: 10

enter incremental fuel cost(in Rs/MW) of generators: 10

multiplying factor for increase generation of 1: 10

multiplying factor for decrease generation of 1: 10

f_lo_cut1 =

-8.6518 +66.9180i

f_up_cut1 =
2.2758e+003 +8.2489e+003i

multiplying factor for increase generation of 2: 10

multiplying factor for decrease generation of 2: 10

f_lo_cut2 =

-8.6518 +66.9180i

f_up_cut2 =

2.2758e+003 +8.2489e+003i

multiplying factor for increase generation of 3: 10

multiplying factor for decrease generation of 3: 10

f_lo_cut3 =

-8.6518 +66.9180i

f_up_cut3 =

2.2758e+003 +8.2489e+003i
multiplying factor for increase generation of 4: 10

multiplying factor for decrease generation of 4: 10

f_lo_cut4 =

-8.6518 +66.9180i

f_up_cut4 =

2.2758e+003 +8.2489e+003i

enter incremental cost(in Rs/MW) of load shedding: 10

enter incremental cost(in Rs/MW) of load shedding: 10

enter incremental cost(in Rs/MW) of load shedding: 10

enter incremental cost(in Rs/MW) of load shedding: 10

enter incremental cost(in Rs/MW) of load shedding: 10

enter incremental cost(in Rs/MW) of load shedding: 10

enter incremental cost(in Rs/MW) of load shedding: 10

enter incremental cost(in Rs/MW) of load shedding: 10

enter incremental cost(in Rs/MW) of load shedding: 10

Fl =
1.0e+002 *

Columns 1 through 4

-0.9889 + 2.7029i -0.9889 + 2.7029i -0.9889 + 2.7029i -0.9889 + 2.7029i

Columns 5 through 8

-0.9889 + 2.7029i -0.9889 + 2.7029i -0.9889 + 2.7029i -0.9889 + 2.7029i

Column 9

-0.9889 + 2.7029i

multiplying factor for load shedding: 10

f_lo_lod_shdng =

-9.8891e+002 +2.7029e+003i

Warning: Matrix is singular, close to singular or badly scaled.

Results may be inaccurate. RCOND = NaN.

> In f_decoupled1 at 485

DP_new =
NaN + NaNi

NaN + NaNi

NaN + NaNi

>>

You might also like