chp9ex8.m

solve power flows by gauss seidel method ans newton Rapson

%        Bus   Bus     R        X       1/2B
%        No.   No.    p.u.     p.u.     p.u.
linedata=[1     2     0.00     0.06     0.0000
          2     3     0.08     0.30     0.0004
          2     5     0.04     0.15     0.0002
          2     6     0.12     0.45     0.0005
          3     4     0.10     0.40     0.0005
          3     6     0.04     0.40     0.0005
          4     6     0.15     0.60     0.0008
          4     9     0.18     0.70     0.0009
          4    10     0.00     0.08     0.0000
          5     7     0.05     0.43     0.0003
          6     8     0.06     0.48     0.0000
          7     8     0.06     0.35     0.0004
          7    11     0.00     0.10     0.0000
          8     9     0.052    0.48     0.0000];
%         Gen.  Ra     Xd'
gendata=[ 1     0     0.20
          10    0     0.15
          11    0     0.25];
Zbus=zbuildpi(linedata, gendata)
symfault(linedata, Zbus)