chp10ex8.m

Those who are working on power system optimisatio, for them here is all the load flow programs from

zdata1 = [ 0     1     0.00    0.20
           0    10     0.00    0.15
           0    11     0.00    0.25
           1     2     0.00    0.06
           2     3     0.00    0.30
           2     5     0.00    0.15
           2     6     0.00    0.45
           3     4     0.00    0.40
           3     6     0.00    0.40
           4     6     0.00    0.60
           4     9     0.00    0.70
           4    10     0.00    0.08
           5     7     0.00    0.43
           6     8     0.00    0.48
           7     8     0.00    0.35
           7    11     0.00    0.10
           8     9    0.000    0.48];

zdata0 = [ 0     1     0.00    0.06+3*0.05
           0    10     0.00    0.04+3*0.05
           0    11     0.00    0.08
           0     2     0.00    0.06
           0     7     0.00    0.10+3*.08
           1     2      inf     inf
           2     3     0.00    0.60
           2     5     0.00    0.30
           2     6     0.00    0.90
           3     4     0.00    0.80
           3     6     0.00    0.80
           4     6     0.00    1.00
           4     9     0.00    1.10
           4    10     0.00    0.08
           5     7     0.00    0.80
           6     8     0.00    0.95
           7     8     0.00    0.70
           7    11      inf     inf
           8     9     0.00   0.90];
zdata2=zdata1;
Zbus0 = zbuild(zdata0)
Zbus1 = zbuild(zdata1)
Zbus2 = Zbus1;
lgfault(zdata0, Zbus0, zdata1, Zbus1, zdata2, Zbus2)
llfault(zdata1, Zbus1,zdata2, Zbus2)
dlgfault(zdata0, Zbus0, zdata1, Zbus1, zdata2, Zbus2)