clear
basemva = 100; accuracy = 0.001; maxiter = 50;
busdata =[...
1 2 0.955 0 51 27 0 0 0 0 0
2 0 1 0 20 9 0 0 0 0 0
3 0 1 0 39 10 0 0 0 0 0
4 2 0.998 0 39 12 0 0 0 0 0
5 0 1 0 0 0 0 0 0 0 0
6 2 0.99 0 52 22 0 0 0 0 0
7 0 1 0 19 2 0 0 0 0 0
8 2 1.015 0 28 0 0 0 0 0 0
9 0 1 0 0 0 0 0 0 0 0
10 2 1.05 0 0 0 450 0 0 0 0
11 0 1 0 70 23 0 0 0 0 0
12 2 0.99 0 47 10 85 0 0 0 0
13 0 1 0 34 16 0 0 0 0 0
14 0 1 0 14 1 0 0 0 0 0
15 2 0.97 0 90 30 0 0 0 0 0
16 0 1 0 25 10 0 0 0 0 0
17 0 1 0 11 3 0 0 0 0 0
18 2 0.973 0 60 34 0 0 0 0 0
19 2 0.962 0 45 25 0 0 0 0 0
20 0 1 0 18 3 0 0 0 0 0
21 0 1 0 14 8 0 0 0 0 0
22 0 1 0 10 5 0 0 0 0 0
23 0 1 0 7 3 0 0 0 0 0
24 2 0.992 0 13 0 0 0 0 0 0
25 2 1.05 0 0 0 220 0 0 0 0
26 2 1.015 0 0 0 314 0 0 0 0
27 2 0.968 0 71 13 0 0 0 0 0
28 0 1 0 17 7 0 0 0 0 0
29 0 1 0 24 4 0 0 0 0 0
30 0 1 0 0 0 0 0 0 0 0
31 2 0.967 0 43 27 7 0 0 0 0
32 2 0.963 0 59 23 0 0 0 0 0
33 0 1 0 23 9 0 0 0 0 0
34 2 0.984 0 59 26 0 0 0 0 0
35 0 1 0 33 9 0 0 0 0 0
36 2 0.98 0 31 17 0 0 0 0 0
37 0 1 0 0 0 0 0 0 0 0
38 0 1 0 0 0 0 0 0 0 0
39 0 1 0 27 11 0 0 0 0 0
40 2 0.97 0 66 23 0 0 0 0 0
41 0 1 0 37 10 0 0 0 0 0
42 2 0.985 0 96 23 0 0 0 0 0
43 0 1 0 18 7 0 0 0 0 0
44 0 1 0 16 8 0 0 0 0 0
45 0 1 0 53 22 0 0 0 0 0
46 2 1.005 0 28 10 19 0 0 0 0
47 0 1 0 34 0 0 0 0 0 0
48 0 1 0 20 11 0 0 0 0 0
49 2 1.025 0 87 30 204 0 0 0 0
50 0 1 0 17 4 0 0 0 0 0
51 0 1 0 17 8 0 0 0 0 0
52 0 1 0 18 5 0 0 0 0 0
53 0 1 0 23 11 0 0 0 0 0
54 2 0.955 0 113 32 48 0 0 0 0
55 2 0.952 0 63 22 0 0 0 0 0
56 2 0.954 0 84 18 0 0 0 0 0
57 0 1 0 12 3 0 0 0 0 0
58 0 1 0 12 3 0 0 0 0 0
59 2 0.985
