Hello Jose,
Thank you for rapid response, i'm quite in a rush to finish this study
case for my thesis research.
- i have a radial topology of wind turbines ( 100 wind turbines with 8 MW
each) like i said i represent all my turbines nodes like a PQ node
with P=8Mw et Q=0MVAR. and i have 3 AC offshores plateformes (ie 3
transformers ) related to my wind farm (each transformer is related to one
cluster in my wind farm) after that i have a long transmission line (
distance about 60km) and finally i have a terrestrial network which is
composed of one final transformer and the output node of the transformer
is my slack bus ( my only generator).
I put you my case file so you can more understand what i mean.
from what I understood from your answer is that I have to add PV nodes and
i have to add also Q for each turbine ?
thank you again and i hope you can help me.
Best regards
> Hello Asma,
>
> Your description is not very clear, so let me ask you a couple of things
> first:
>
> - You seem to be saying that your generators should have P=8 MW.
> That's
> fine, but then you total load should be balanced according to that (is
> it?). Otherwise the imbalance, which has to be picked up by the slack
> bus,
> could render the case infeasible.
> - You also seem to be saying that all your generators are PQ type. A
> case with no PV buses is quite prone to numerical ill-conditioning
> (unless
> the total number of buses is low). Worse yet, if all gens have Q=0
> then
> all reactive power will have to be provided by the slack, and then
> voltages
> will quickly drop as you move away from the slack bus. And if the
> network
> topology is dominantly radial, things are even worse yet.
>
> Maybe we can help you faster if you could share your case file here.
>
> --
> Jose L. Marin
> Grupo AIA
>
>
>
> 2018-05-18 17:31 GMT+02:00 Asma DABBABI <[email protected]>:
>
>> the library Matpower works perfectly with small wind farms (30 wind
>> turbines with
>> 3 mw each, the turbines are represented as PQnode with P=3MW and Q=0MVAR
>> and i have one slack bus in the terrestrial network)
>> but when I tried to make an example of a larger farm such as Borsselle
>> Park (100 wind turbines of 8 mw each) I had a divergence of the system.
>>
>> "Newton's method power flow did not converge in 10 iterations ",
>> I looked on the internet and I found that in this case we have to do a
>> continuation power flow that gradually increases the loading/generation
>> like this:
>>
>> {define_constants;
>> mpcbase = loadcase ('casefile');
>> mpcbase.bus (:, PD) = 0;
>> mpcbase.bus (:, QD) = 0;
>> mpcbase.gen (:, PG) = 0;
>> mpctarget = loadcase ('casefile');
>> results = runcpf (mpcbase, mpctarget);
>> results.cpf.max_lam}
>>
>> If the resulting value of results is greater than 1, it indicates that
>> the
>> load for the box is greater than that of the loading and unloading
>> system
>> at least by a factor of results.cpf.max_lam to get a convergent power
>> flow
>> solution.
>> I tried this and actually I found that it shows that results.cpf.max_lam
>> is less than 1 so i must reduce the loads to 1.35 mw for each turbine.
>> in my case i don't need to reduce the loads, i must found the losses of
>> the network with 8mw for each turbine.
>> I've tried many solutions that i found in matpower page but until now I
>> have not found the right solution to calculate with 8 mw.
>>
>> So can you help me please on that, I 'll be so grateful if you answer
>> me.
>>
>> Best regards,
>>
>>
>>
>>
>
function [baseMVA, bus, gen, branch] = caseborselle
%case 5 nodes Power flow data for 5 bus, 2 generator case.
% Please see 'help caseformat' for details on the case file format.
%
% case file can be used together with dc case files "case5_stagg_....m"
%
% Network data from ...
% G.W. Stagg, A.H. El-Abiad, "Computer methods in power system analysis",
% McGraw-Hill, 1968.
%
% MATPOWER case file data provided by Jef Beerten.
\%% MATPOWER Case Format : Version 1
\%%----- Power Flow Data -----\%%
\%% system MVA base
baseMVA = 100;
\%% bus data
% bus_i type Pd Qd Gs Bs area Vm Va
baseKV zone Vmax Vmin
bus =[4 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
5 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
6 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
7 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
8 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
9 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
10 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
11 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
12 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
13 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
14 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
15 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
16 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
17 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
18 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
19 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
20 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
21 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
22 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
23 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
24 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
25 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
26 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
27 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
28 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
29 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
30 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
31 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
32 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
33 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
34 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
35 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
36 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
37 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
38 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
39 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
40 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
41 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
42 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
43 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
44 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
45 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
46 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
47 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
48 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
49 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
50 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
51 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
52 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
53 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
54 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
55 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
56 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
57 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
58 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
59 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
60 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
61 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
62 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
63 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
64 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
65 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
66 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
67 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
68 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
69 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
70 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
71 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
72 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
73 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
74 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
75 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
76 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
77 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
78 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
79 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
80 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
81 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
82 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
83 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
84 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
85 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
86 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
87 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
88 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
89 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
90 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
91 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
92 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
93 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
94 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
95 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
96 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
97 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
98 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
99 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
100 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
101 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
102 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
103 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
1 1 0 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
2 1 0 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
3 1 0 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;
1001 1 0 0 0 0 1 1 0 132 1 1.10000000000000 0.900000000000000;
1002 1 0 0 0 0 1 1 0 132 1 1.10000000000000 0.900000000000000;
1003 1 0 0 0 0 1 1 0 132 1 1.10000000000000 0.900000000000000;
3001 1 0 0 0 0 1 1 0 132 1 1.10000000000000 0.900000000000000;
3002 3 0 0 0 0 1 1 0 230 1 1.10000000000000 0.900000000000000];
%bus = 0 33 1 1.10000000000000 0.900000000000000;16 1 3.89377814200271 0 0 0 1
1 0 33 1 1.10000000000000 0.900000000000000;17 1 3.89377814200271 0 0 0 1 1 0
33 1 1.10000000000000 0.900000000000000;18 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;19 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;20 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;21 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;22 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;23 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;24 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;25 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;26 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;27 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;28 1 3.89377814200271 0 0 0 1 1 0 33 1
1.1000000000!
0000 0.90
0000000000000;29 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;30 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;31 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;32 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;33 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;34 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;35 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;36 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;37 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;38 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;39 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;40 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;41 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;42 1 3.89377!
814200271
0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;43 1 3.89377814200271 0 0
0 1 1 0 33 1 1.10000000000000 0.900000000000000;44 1 3.89377814200271 0 0 0 1 1
0 33 1 1.10000000000000 0.900000000000000;45 1 3.89377814200271 0 0 0 1 1 0 33
1 1.10000000000000 0.900000000000000;46 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;47 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;48 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;49 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;50 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;51 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;52 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;53 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;54 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;55 1 3.89377814200271 0 0 0 1 1 0 33 1
1.100000!
00000000
0.900000000000000;56 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;57 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;58 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;59 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;60 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;61 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;62 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;63 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;64 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;65 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;66 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;67 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;68 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;69 1 3.8!
937781420
0271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;70 1 3.89377814200271
0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;71 1 3.89377814200271 0 0 0
1 1 0 33 1 1.10000000000000 0.900000000000000;72 1 3.89377814200271 0 0 0 1 1 0
33 1 1.10000000000000 0.900000000000000;73 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;74 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;75 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;76 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;77 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;78 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;79 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;80 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;81 1 3.89377814200271 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;82 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10!
000000000
000 0.900000000000000;83 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;84 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;85 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;86 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;87 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;88 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;89 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;90 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;91 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;92 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;93 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;94 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;95 1 3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;96 1!
3.893778
14200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;97 1
3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;98 1
3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;99 1
3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;100 1
3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;101 1
3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;102 1
3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;103 1
3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;104 1
3.89377814200271 0 0 0 1 1 0 33 1 1.10000000000000 0.900000000000000;1 1 0 0 0
0 1 1 0 33 1 1.10000000000000 0.900000000000000;2 1 0 0 0 0 1 1 0 33 1
1.10000000000000 0.900000000000000;3 1 0 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;4 1 0 0 0 0 1 1 0 33 1 1.10000000000000
0.900000000000000;1001 1 0 0 0 0 1 1 0 132 1 1.10000000000000
0.900000000000000;1002 1 0 0 0 0 1 1 0 132 1 1.1000000!
0000000 0
.900000000000000;1003 1 0 0 0 0 1 1 0 132 1 1.10000000000000
0.900000000000000;1004 1 0 0 0 0 1 1 0 132 1 1.10000000000000
0.900000000000000;3001 1 0 0 0 0 1 1 0 132 1 1.10000000000000
0.900000000000000;3002 3 0 0 0 0 1 1 0 230 1 1.10000000000000
0.900000000000000];
bus(1:100,3)=8;
\%% generator data
% bus Pg Qg Qmax Qmin Vg mBase status
Pmax Pmin
gen = [3002 0 0 200 -200 1 100 1 1000 100];
\%% branch data commençant par 1 en ordre
% fbus tbus r x b rateA rateB rateC
ratio angle status
%branch = [26 1 0.0172855936562986 0.0737162362582345 8.42343855970121e-05 0 0
0 0 0 1;35 1 0.0243547816607643 0.103863533681322 0.000118683228955594 0 0 0 0
0 1;37 1 0.00617786852890702 0.0263461715636985 3.01053565286129e-05 0 0 0 0 0
1;48 1 0.00370216740417866 0.0157882831484080 1.80410232283021e-05 0 0 0 0 0
1;57 1 0.0118059272896418 0.0503476214144171 5.75314363752619e-05 0 0 0 0 0
1;75 1 0.0318041822635946 0.135632287809149 0.000154984885428745 0 0 0 0 0 1;95
1 0.0393009627115836 0.167603098281064 0.000191517428513369 0 0 0 0 0 1;39 2
0.00410718597111310 0.0175155275209628 2.00147182496861e-05 0 0 0 0 0 1;92 2
0.0397570510418288 0.169548135042629 0.000193739991477655 0 0 0 0 0 1;101 2
0.0392660726500572 0.167454305935850 0.000191347405836998 0 0 0 0 0 1;8 3
0.0152644272511473 0.0650967590680456 7.43850443645110e-05 0 0 0 0 0 1;22 4
0.0181135748789532 0.0772472494614949 8.82692189364361e-05 0 0 0 0 0 1;32 4
0.0200580430735703 0.0855396390479402 9.77448023004368e-05 0 0 !
0 0 0 1;3
4 4 0.00976266960402693 0.0416339336300028 4.75744421761401e-05 0 0 0 0 0 1;62
4 0.0308038947833482 0.131366456407221 0.000150110386872646 0 0 0 0 0 1;6 5
0.00912638238744698 0.0389204197225177 4.44737524447050e-05 0 0 0 0 0 1;16 6
0.00538399440034720 0.0229606116585059 2.62367303888197e-05 0 0 0 0 0 1;17 7
0.00556250990428545 0.0237219098427470 2.71066538691915e-05 0 0 0 0 0 1;19 9
0.00490501493335039 0.0209179532312604 2.39026166801269e-05 0 0 0 0 0 1;20 9
0.00921902020705440 0.0393154834693818 4.49251855845113e-05 0 0 0 0 0 1;11 10
0.00976268208386926 0.0416339868516059 4.75745029916303e-05 0 0 0 0 0 1;20 10
0.00499802272865703 0.0213145947784938 2.43558527476387e-05 0 0 0 0 0 1;12 11
0.00961404427872044 0.0410001052633237 4.68501764546220e-05 0 0 0 0 0 1;13 12
0.00960528119464459 0.0409627341675470 4.68074730903210e-05 0 0 0 0 0 1;24 13
0.00975266362101406 0.0415912620606926 4.75256820439664e-05 0 0 0 0 0 1;23 14
0.0123099686965883 0.0524971591264103 5.99876793642972e-0!
5 0 0 0 0
0 1;24 14 0.00477365454099451 0.0203577529910822 2.32624846625697e-05 0 0 0 0
0 1;25 15 0.00486685321703487 0.0207552085693257 2.37166509105339e-05 0 0 0 0 0
1;26 16 0.00534996781921228 0.0228155017164422 2.60709155367031e-05 0 0 0 0 0
1;27 17 0.00548373439278381 0.0233859633700552 2.67227731102560e-05 0 0 0 0 0
1;27 18 0.0121606422830955 0.0518603408952573 5.92599971715949e-05 0 0 0 0 0
1;28 18 0.00482882573552007 0.0205930363658486 2.35313392802311e-05 0 0 0 0 0
1;29 19 0.00542304205266040 0.0231271344878911 2.64270134110294e-05 0 0 0 0 0
1;30 21 0.00899060813334644 0.0383413961036550 4.38121112479109e-05 0 0 0 0 0
1;31 21 0.00495765037258990 0.0211424226103101 2.41591143146963e-05 0 0 0 0 0
1;33 23 0.00479546404119598 0.0204507619036769 2.33687644864611e-05 0 0 0 0 0
1;26 25 0.00969518467342325 0.0413461370297143 4.72455815202965e-05 0 0 0 0 0
1;37 28 0.0118528434397935 0.0505477002822694 5.77600633556933e-05 0 0 0 0 0
1;39 29 0.00485969713995452 0.0207246907242793 2.36!
817786482
250e-05 0 0 0 0 0 1;41 31 0.00495299366745768 0.0211225636004006
2.41364217359213e-05 0 0 0 0 0 1;42 32 0.00525647592980326 0.0224167956988799
2.56153204313651e-05 0 0 0 0 0 1;43 33 0.00546519462231459 0.0233068985645728
2.66324269985958e-05 0 0 0 0 0 1;44 34 0.00508083533829873 0.0216677578817629
2.47594066799800e-05 0 0 0 0 0 1;45 35 0.00473904133109186 0.0202101413088008
2.30938111111914e-05 0 0 0 0 0 1;46 36 0.00533691081518869 0.0227598187464172
2.60072874812522e-05 0 0 0 0 0 1;47 37 0.00560820282914347 0.0239167721373918
2.73293199533351e-05 0 0 0 0 0 1;48 38 0.00535737153202996 0.0228470756301942
2.61069945521384e-05 0 0 0 0 0 1;49 39 0.00557191200112050 0.0237620060757950
2.71524712050593e-05 0 0 0 0 0 1;50 40 0.00548022412965790 0.0233709934829305
2.67056672553126e-05 0 0 0 0 0 1;52 41 0.00919319218281090 0.0392053371374088
4.47993231006074e-05 0 0 0 0 0 1;53 43 0.00503658551663129 0.0214790498527950
2.45437729785841e-05 0 0 0 0 0 1;54 44 0.00547533167466701 0.0233!
501290928
978 2.66818258444584e-05 0 0 0 0 0 1;55 45 0.00488615739105893
0.0208375332543451 2.38107220353507e-05 0 0 0 0 0 1;57 46 0.0121866412009225
0.0519712160209315 5.93866924366168e-05 0 0 0 0 0 1;56 47 0.00970701353612190
0.0413965824615980 4.73032247231591e-05 0 0 0 0 0 1;59 49 0.00476351635453962
0.0203145176262577 2.32130802901151e-05 0 0 0 0 0 1;60 50 0.00559013835758562
0.0238397342942192 2.72412900196767e-05 0 0 0 0 0 1;52 51 0.00910037851010055
0.0388095234464502 4.43470330114608e-05 0 0 0 0 0 1;61 51 0.00547321020917285
0.0233410818796698 2.66714877359723e-05 0 0 0 0 0 1;63 53 0.00528756049852834
0.0225493591188980 2.57667985697599e-05 0 0 0 0 0 1;65 55 0.00558430814001594
0.0238148707168178 2.72128788002159e-05 0 0 0 0 0 1;66 56 0.00525842797697354
0.0224251204098069 2.56248329478174e-05 0 0 0 0 0 1;67 57 0.00527890492596403
0.0225124465172953 2.57246191195517e-05 0 0 0 0 0 1;68 58 0.00477051753406135
0.0203443748943372 2.32471977214971e-05 0 0 0 0 0 1;68 59 0.00904564!
380578758
0.0385761015302067 4.40803054532229e-05 0 0 0 0 0 1;70 60 0.00555041126079265
0.0236703139022284 2.70476959981784e-05 0 0 0 0 0 1;71 61 0.00507105233662308
0.0216260372398703 2.47117331576337e-05 0 0 0 0 0 1;72 62 0.00530792343417993
0.0226361990044843 2.58660291433617e-05 0 0 0 0 0 1;73 63 0.00483380522758819
0.0206142719346719 2.35556048312606e-05 0 0 0 0 0 1;74 64 0.00514500330096878
0.0219414089226499 2.50721034272250e-05 0 0 0 0 0 1;76 66 0.00540951627108247
0.0230694523665708 2.63611009567863e-05 0 0 0 0 0 1;78 67 0.00939639400113749
0.0400719127115956 4.57895454012859e-05 0 0 0 0 0 1;70 69 0.00985454357230561
0.0420257398523563 4.80221529938414e-05 0 0 0 0 0 1;78 69 0.00923629126735668
0.0393891377266199 4.50093491476552e-05 0 0 0 0 0 1;81 71 0.00517699618529603
0.0220778459503010 2.52280078762342e-05 0 0 0 0 0 1;82 72 0.00528202002260934
0.0225257311753082 2.57397992896523e-05 0 0 0 0 0 1;84 73 0.00918044359359009
0.0391509694348210 4.47371979806030e-05 0 0 0 0 0 1!
;84 74 0.
00471647389443944 0.0201139000963154 2.29838377889721e-05 0 0 0 0 0 1;85 75
0.00474833756160434 0.0202497860637566 2.31391125501329e-05 0 0 0 0 0 1;86 76
0.00475035115993407 0.0202583732660913 2.31489250113949e-05 0 0 0 0 0 1;87 77
0.00483335808981308 0.0206123650684963 2.35534258852255e-05 0 0 0 0 0 1;88 77
0.00956213680505146 0.0407787403701859 4.65972262673540e-05 0 0 0 0 0 1;80 79
0.00950027454349908 0.0405149222347640 4.62957654267766e-05 0 0 0 0 0 1;89 79
0.00560094358214886 0.0238858143133707 2.72939449339608e-05 0 0 0 0 0 1;90 80
0.00478171292400633 0.0203921188148873 2.33017539497810e-05 0 0 0 0 0 1;91 81
0.00532137316179225 0.0226935567856513 2.59315709792048e-05 0 0 0 0 0 1;83 82
0.00981674830956982 0.0418645579703410 4.78379729883218e-05 0 0 0 0 0 1;94 83
0.00942511359449478 0.0401943904448607 4.59294987838049e-05 0 0 0 0 0 1;96 86
0.00566373043163182 0.0241535754514852 2.75999113817977e-05 0 0 0 0 0 1;97 87
0.00533756650267599 0.0227626149948242 2.6010482710397!
0e-05 0 0
0 0 0 1;98 88 0.00515332669286922 0.0219769049047957 2.51126641286252e-05 0 0
0 0 0 1;99 89 0.00526167324949368 0.0224389602165601 2.56406474776661e-05 0 0 0
0 0 1;100 90 0.00478462973284901 0.0204045578536639 2.33159678440592e-05 0 0 0
0 0 1;100 91 0.00966819736102520 0.0412310467912031 4.71140697120262e-05 0 0 0
0 0 1;102 93 0.00941354720968026 0.0401450643775822 4.58731346612991e-05 0 0 0
0 0 1;103 93 0.00492175213670045 0.0209893307992507 2.39841787062185e-05 0 0 0
0 0 1;104 94 0.00495102139394362 0.0211141526320984 2.41268106545208e-05 0 0 0
0 0 1;97 96 0.00917564019725510 0.0391304848448142 4.47137905612744e-05 0 0 0 0
0 1;99 98 0.00960373767819124 0.0409561516789264 4.67999513839398e-05 0 0 0 0 0
1;104 103 0.00942042666040345 0.0401744025203686 4.59066588963607e-05 0 0 0 0 0
1;1 1001 0.0112000000000000 0.160000000000000 0 0 0 0 1 0 1;2 1002
0.0224000000000000 0.320000000000000 0 0 0 0 1 0 1;3 1003 0.140000000000000 2 0
0 0 0 1 0 1;4 1004 0.0448000000000000 0.6400000!
00000000
0 0 0 0 1 0 1;1001 3001 0.00309401045964875 0.00644326171506984
0.00982439919822515 0 0 0 0 0 1;1002 3001 0.00342839651484387
0.0214188571398819 0.00362872463974085 0 0 0 0 0 1;1003 1001
0.00245162227766026 0.0100697507372565 0.00135557030348993 0 0 0 0 0 1;1004
1003 0.00428430014227493 0.0139743711039425 0.00174447370041669 0 0 0 0 0
1;3001 3002 0.00700000000000000 0.100000000000000 0 0 0 0 1 0 1];
branch=[54 1 0.0952946876361481 0.406394240638720 0.000464380317118305 0 0 0
0 0 1;
64 1 0.0866953548078592 0.369721479318640 0.000422474928634805 0 0 0 0 0
1;
74 1 0.0871691034430874 0.371741830312415 0.000424783551989567 0 0 0 0 0
1;
84 1 0.0814033758962182 0.347153277411268 0.000396686598706547 0 0 0 0 0
1;
94 1 0.0854513172265461 0.364416149916500 0.000416412612037230 0 0 0 0 0
1;
48 2 0.0750872359963269 0.320217433010530 0.000365907432286755 0 0 0 0 0
1;
4 3 0.0385583103764922 0.164435979113515 0.000187898411174279 0 0 0 0 0 1;
14 3 0.0508967107636467 0.217054388181390 0.000248024640942685 0 0 0 0 0
1;
24 3 0.0629109274563820 0.268290281711525 0.000306571091915930 0 0 0 0 0
1;
34 3 0.0751968524576016 0.320684903964539 0.000366441604005146 0 0 0 0 0
1;
5 4 0.00778846923097465 0.0332147480344888 3.79539710036702e-05 0 0 0 0 0
1;
6 5 0.00799929185778330 0.0341138233497728 3.89813302481531e-05 0 0 0 0 0
1;
7 6 0.00767181905478493 0.0327172810617862 3.73855233057490e-05 0 0 0 0 0
1;
8 7 0.00808921620858120 0.0344973152228663 3.94195403895162e-05 0 0 0 0 0
1;
9 8 0.00794032413156135 0.0338623492653867 3.86939747608541e-05 0 0 0 0 0
1;
10 9 0.00842476892739692 0.0359283151635770 4.10547215505954e-05 0 0 0 0
0 1;
11 10 0.00827613112224812 0.0352944335752949 4.03303950135872e-05 0 0 0 0
0 1;
12 11 0.00826736803817227 0.0352570624795182 4.02876916492862e-05 0 0 0 0
0 1;
13 12 0.00802225598854363 0.0342117563060978 3.90932366994922e-05 0 0 0 0
0 1;
15 14 0.00784895602989983 0.0334727003645465 3.82487290805694e-05 0 0 0 0
0 1;
16 15 0.00828811685153586 0.0353455479812835 4.03888026426545e-05 0 0 0 0
0 1;
17 16 0.00787078583953579 0.0335657958888654 3.83551078743180e-05 0 0 0 0
0 1;
18 17 0.00817276553318498 0.0348536201247499 3.98266844057982e-05 0 0 0 0
0 1;
19 18 0.00763379291610701 0.0325551145849786 3.72002182192214e-05 0 0 0 0
0 1;
20 19 0.00806914921997029 0.0344117374346314 3.93217519947389e-05 0 0 0 0
0 1;
21 20 0.00818669176570296 0.0349130100113165 3.98945482965710e-05 0 0 0 0
0 1;
22 21 0.00804860996008835 0.0343241455957781 3.92216621759468e-05 0 0 0 0
0 1;
23 22 0.00853853085678045 0.0364134649033184 4.16090945398309e-05 0 0 0 0
0 1;
25 24 0.00835727151695092 0.0356404653416854 4.07258000792616e-05 0 0 0 0
0 1;
26 25 0.00852100132302608 0.0363387083587977 4.15236714103184e-05 0 0 0 0
0 1;
27 26 0.00850559927049217 0.0362730246822039 4.14486157045137e-05 0 0 0 0
0 1;
28 27 0.00829942387107845 0.0353937679580197 4.04439028528666e-05 0 0 0 0
0 1;
29 28 0.00796214820207527 0.0339554203146147 3.88003255873006e-05 0 0 0 0
0 1;
30 29 0.00814102080349137 0.0347182413786925 3.96719892385478e-05 0 0 0 0
0 1;
31 30 0.00808772509363267 0.0344909562059867 3.94122740407844e-05 0 0 0 0
0 1;
32 31 0.00847236122219408 0.0361312775215438 4.12866434499008e-05 0 0 0 0
0 1;
33 32 0.00851770340489966 0.0363246440393071 4.15076003332903e-05 0 0 0 0
0 1;
35 34 0.00832381910020646 0.0354978039841887 4.05627831862793e-05 0 0 0 0
0 1;
36 35 0.00823974060712040 0.0351392423875316 4.01530604803167e-05 0 0 0 0
0 1;
37 36 0.00789166794496100 0.0336548498286781 3.84568684129168e-05 0 0 0 0
0 1;
38 37 0.00856373691227392 0.0365209587839939 4.17319261092889e-05 0 0 0 0
0 1;
39 38 0.00796137816677984 0.0339521364179263 3.87965731301222e-05 0 0 0 0
0 1;
40 39 0.00809048800271292 0.0345027389230035 3.94257379664940e-05 0 0 0 0
0 1;
41 40 0.00765612063323099 0.0326503334882657 3.73090233647731e-05 0 0 0 0
0 1;
42 41 0.00848109373375153 0.0361685182375931 4.13291978313347e-05 0 0 0 0
0 1;
43 42 0.00825481874670202 0.0352035447031960 4.02265377266791e-05 0 0 0 0
0 1;
45 44 0.00819130500176282 0.0349326836427847 3.99170290463122e-05 0 0 0 0
0 1;
46 45 0.00810462163057197 0.0345630132688640 3.94946125150132e-05 0 0 0 0
0 1;
47 46 0.00848106900847798 0.0361684127940577 4.13290773426625e-05 0 0 0 0
0 1;
48 47 0.00821774669594069 0.0350454470348811 4.00458820037232e-05 0 0 0 0
0 1;
49 48 0.00806943196316117 0.0344129432227735 3.93231298299142e-05 0 0 0 0
0 1;
50 49 0.00816669721088600 0.0348277411246313 3.97971128787438e-05 0 0 0 0
0 1;
51 50 0.00776246535362823 0.0331038517584214 3.78272515704260e-05 0 0 0 0
0 1;
52 51 0.00792266054275477 0.0337870210289649 3.86078983679836e-05 0 0 0 0
0 1;
53 52 0.00773921402175262 0.0330046940026658 3.77139455600107e-05 0 0 0 0
0 1;
55 54 0.00818591550834969 0.0349096995800065 3.98907655186972e-05 0 0 0 0
0 1;
56 55 0.00831333191488403 0.0354530802769048 4.05116781082679e-05 0 0 0 0
0 1;
57 56 0.00784336532210685 0.0334488582018350 3.82214850156355e-05 0 0 0 0
0 1;
58 57 0.00791208145340866 0.0337419053872842 3.85563454326916e-05 0 0 0 0
0 1;
59 58 0.00792812028227456 0.0338103046636697 3.86345042117349e-05 0 0 0 0
0 1;
60 59 0.00844102746506933 0.0359976514113278 4.11339509921045e-05 0 0 0 0
0 1;
61 60 0.00791856977878127 0.0337695755347838 3.85879636757316e-05 0 0 0 0
0 1;
62 61 0.00772443292474394 0.0329416584563656 3.76419158311086e-05 0 0 0 0
0 1;
63 62 0.00788457541577832 0.0336246029903401 3.84223058257178e-05 0 0 0 0
0 1;
65 64 0.00838835802183494 0.0357730370186574 4.08772876527479e-05 0 0 0 0
0 1;
66 65 0.00829800851842243 0.0353877320374228 4.04370057011833e-05 0 0 0 0
0 1;
67 66 0.00854184266512575 0.0364275884591116 4.16252333052519e-05 0 0 0 0
0 1;
68 67 0.00842679420619334 0.0359369521784932 4.10645909319120e-05 0 0 0 0
0 1;
69 68 0.00851663041583327 0.0363200681643274 4.15023715528065e-05 0 0 0 0
0 1;
70 69 0.00832853386538781 0.0355179105972979 4.05857586973425e-05 0 0 0 0
0 1;
71 70 0.00783209060567815 0.0334007759851332 3.81665422216513e-05 0 0 0 0
0 1;
72 71 0.00823256118806760 0.0351086249981832 4.01180744702970e-05 0 0 0 0
0 1;
73 72 0.00854150289169390 0.0364261394594953 4.16235775561441e-05 0 0 0 0
0 1;
75 74 0.00848267279522405 0.0361752523116948 4.13368927520647e-05 0 0 0 0
0 1;
76 75 0.00787516298314277 0.0335844626791814 3.83764381224851e-05 0 0 0 0
0 1;
77 76 0.00832962459310597 0.0355225621206279 4.05910739200083e-05 0 0 0 0
0 1;
78 77 0.00822505194819867 0.0350766010531973 4.00814812112377e-05 0 0 0 0
0 1;
79 78 0.00816236138702674 0.0348092505467351 3.97759839857417e-05 0 0 0 0
0 1;
80 79 0.00809766354510279 0.0345333397799142 3.94607050851568e-05 0 0 0 0
0 1;
81 80 0.00788864223772347 0.0336419463812175 3.84421238461279e-05 0 0 0 0
0 1;
82 81 0.00847883515309749 0.0361588862823122 4.13181915472870e-05 0 0 0 0
0 1;
83 82 0.00777636088183007 0.0331631107031951 3.78949658360683e-05 0 0 0 0
0 1;
85 84 0.00819922846279961 0.0349664740776049 3.99556408455658e-05 0 0 0 0
0 1;
86 85 0.00854723994917187 0.0364506057458654 4.16515347973821e-05 0 0 0 0
0 1;
87 86 0.00804152894932795 0.0342939478789661 3.91871557191443e-05 0 0 0 0
0 1;
88 87 0.00827539797412397 0.0352913069878367 4.03268223111960e-05 0 0 0 0
0 1;
89 88 0.00784247779758031 0.0334450732599358 3.82171600219626e-05 0 0 0 0
0 1;
90 89 0.00804861336291936 0.0343241601075096 3.92216787582746e-05 0 0 0 0
0 1;
91 90 0.00826679207382245 0.0352546062188348 4.02848849187742e-05 0 0 0 0
0 1;
92 91 0.00770663095574618 0.0328657400830315 3.75551650980554e-05 0 0 0 0
0 1;
93 92 0.00770987288759666 0.0328795656431467 3.75709633484946e-05 0 0 0 0
0 1;
95 94 0.00826026935503344 0.0352267893969855 4.02530990732580e-05 0 0 0 0
0 1;
96 95 0.00783772704078277 0.0334248131567854 3.81940091202395e-05 0 0 0 0
0 1;
97 96 0.00851878050615132 0.0363292374511343 4.15128491528543e-05 0 0 0 0
0 1;
98 97 0.00826582452171891 0.0352504799908975 4.02801699429050e-05 0 0 0 0
0 1;
99 98 0.00777030294704574 0.0331372760017276 3.78654449283869e-05 0 0 0 0
0 1;
100 99 0.00836558417931414 0.0356759155665886 4.07663084945792e-05 0 0 0
0 0 1;
101 100 0.00768329582781305 0.0327662249701556 3.74414507412703e-05 0 0 0
0 0 1;
102 101 0.00783601129551909 0.0334174961802481 3.81856481260498e-05 0 0 0
0 0 1;
103 102 0.00808251350393111 0.0344687308323397 3.93868774553258e-05 0 0 0
0 0 1;
1 1001 0.0140000000000000 0.200000000000000 0 0 0 0 1 0 1;
2 1002 0.0700000000000000 1 0 0 0 0 1 0 1;
3 1003 0.0155555555555556 0.222222222222222 0 0 0 0 1 0 1;
% 1 1001 0.007 0.1 0 0 0 0 1 0 1;
% 2 1002 0.007 0.1 0 0 0 0 1 0 1;
% 3 1003 0.007 0.1 0 0 0 0 1 0 1;
1001 3001 0.0429500840918572 0.0894433409638538 0.136379232461822 0 0 0 0
0 1;
1003 3001 0.0582715357650746 0.182025284191857 0.123352918316623 0 0 0 0
0 1;
1002 1001 0.0263410520633372 0.0704279622553246 0.00751697994832171 0 0 0
0 0 1;
3001 3002 0.007 0.1 0 0 0 0 1 0 1];
return;
