I think there is some confusion about bus numbering. There are external bus numbers and there are row indices into the bus matrix. In the original case file, these need not be the same. We call the bus numbers the external bus numbers. They need not be consecutive, i.e. numbering from 101 to 509 is fine. However, these are not row indices into the bus matrix (which has only 59 rows). So assigning a value to bus(101, PD) simply creates a bunch of new blank rows in the bus matrix and assigns a value to the PD column of row 101. But now you have an invalid bus matrix. Specifically, lots of rows with the BUS_I column equal to zero.
Internally, MATPOWER calls ext2int() to convert from external bus numbers to consecutive integers (1–59 in your case). After this conversion, bus numbers and row indices into the bus matrix are identical. So the bottom line is, when changing the load, you need to find the row index that corresponds to the bus number you are trying to change and use that to index the bus matrix. I often do the following to create and external to internal bus mapping. nb = size(mpc.bus, 1); max_bus_num = max(mpc.bus(:, BUS_I)); e2i = sparse(mpc.bus(:, BUS_I), 1, 1:nb, max_bus_num, 1); % ext to int bus number map Then you could use something like the following to update your loads … mpc.bus(e2i(bus_nbr), PD) = MW Hope this helps, Ray > On Apr 10, 2019, at 7:49 AM, Tayeb Meridji <[email protected]> > wrote: > > Hi Ray, > > One last note. So below i took a snapshot of the bug. > > - when i run the OPF on the original base case, it works fine; > - when I attempt to change the load at bus, it allows me to do that (at least > i think so since no error is shown). But when right after i attempt to run > the OPF again, it gives me the following error message (see below). > > Thanks. > > > <pastedImage.png> > > > From: Tayeb Meridji > Sent: Wednesday, April 10, 2019 8:43 AM > To: MATPOWER discussion forum > Subject: Re: Problems creating a new case in Matpower - 14 Generator Model > > Hi Ray, > > One other note. When I simply run an OPF on my networks fine, and the buses > correspond to what I have defined. So it really seems to be a problem with > the load change function, which for some reason is not working. > > Thanks. > > > From: Tayeb Meridji > Sent: Wednesday, April 10, 2019 8:36 AM > To: MATPOWER discussion forum > Subject: Re: Problems creating a new case in Matpower - 14 Generator Model > > Hi Ray, > > No, bus numbers in the base case i have created go from 101 to 509. And that > is the strange thing because the error message mentions some bus numbers that > are not even in my base case. I don't know if you can open files sent by > users, but in case you can I have attached my base case. > > Thank you. > > > > From: [email protected] > <mailto:[email protected]> > <[email protected] > <mailto:[email protected]>> on behalf of Ray > Zimmerman <[email protected] <mailto:[email protected]>> > Sent: Tuesday, April 9, 2019 5:32 PM > To: MATPOWER discussion forum > Subject: Re: Problems creating a new case in Matpower - 14 Generator Model > > And can you confirm that bus_nbr is always a number between 1 and 59? > > Ray > > >> On Apr 9, 2019, at 3:36 PM, Tayeb Meridji <[email protected] >> <mailto:[email protected]>> wrote: >> >> Hi, >> Thanks for the reply. To change the load i used: mpc.bus(bus_nbr, PD) = MW. >> Thanks for the tips. >> Regards. >> >> Get Outlook for Android <https://aka.ms/ghei36> >> From: [email protected] >> <mailto:[email protected]> >> <[email protected] >> <mailto:[email protected]>> on behalf of Ray >> Zimmerman <[email protected] <mailto:[email protected]>> >> Sent: Tuesday, April 9, 2019 11:27:35 AM >> To: MATPOWER discussion forum >> Subject: Re: Problems creating a new case in Matpower - 14 Generator Model >> >> You haven’t shown how you “attempt to change the load at a given bus” but I >> suspect there is a problem there. Your original case has 59 rows in the bus >> matrix. Make sure that mpc.bus still has only 59 rows after updating the >> load. >> >> Two other tips … >> If you aren’t already, you might consider using the scale_load() >> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/scale_load.html> >> or apply_changes() >> <http://www.pserc.cornell.edu/matpower/docs/ref/matpower6.0/apply_changes.html> >> functions to update your loads. >> I notice that your line ratings are all set to 9999. If you actually mean >> “unconstrained” you should simply use 0, rather than a large number. This >> will eliminate line constraints completely and speed up the optimization. >> Ray >> >> >>> On Apr 2, 2019, at 10:35 AM, Tayeb Meridji <[email protected] >>> <mailto:[email protected]>> wrote: >>> >>> >>> Hi, >>> >>> I am trying to create a new MAtpower case for the 14 generator system >>> representing an equivalent of the Australian network. >>> >>> When I run an OPF of the created network, it seems to work fine. However, >>> when I attempt to change the load at a given bus, I get an error message >>> saying the following (see below). The problem is that the buses mentioned >>> in the error message are not even part of the newly created network. So I >>> am really puzzled. Please note that at the end of this message I have >>> included the input of my new network. >>> >>> Your guidance would be greatly appreciated. >>> >>> Many thanks. >>> >>> Tayeb. >>> >>> >>> >>> Error messages: >>> >>> error: ext2int: bus 60 has an invalid BUS_TYPEext2int: bus 61 has an >>> invalid BUS_TYPEext2int: bus 62 has an invalid BUS_TYPEext2int: bus 63 has >>> an invalid >>> BUS_TYPEext2int: bus 64 has an invalid BUS_TYPEext2int: bus 65 has an >>> invalid BUS_TYPEext2int: bus 66 has an invalid BUS_TYPEext2int: bus 67 has >>> an invalid >>> BUS_TYPEext2int: bus 68 has an invalid BUS_TYPEext2int: bus 69 has an >>> invalid BUS_TYPEext2int: bus 70 has an invalid BUS_TYPEext2int: bus 71 has >>> an invali >>> d BUS_TYPEext2int: bus 72 has an invalid BUS_TYPEext2int: bus 73 has an >>> invalid BUS_TYPEext2int: bus 74 has an invalid BUS_TYPEext2int: bus 75 has >>> an inval >>> id BUS_TYPEext2int: bus 76 has an invalid BUS_TYPEext2int: bus 77 has an >>> invalid BUS_TYPEext2int: bus 78 has an invalid BUS_TYPEext2int: bus 79 has >>> an inva >>> lid BUS_TYPEext2int: bus 80 has an invalid BUS_TYPEext2int: bus 81 has an >>> invalid BUS_TYPEext2int: bus 82 has an invalid BUS_TYPEext2int: bus 83 has >>> an inv >>> alid BUS_TYPEext2int: bus 84 has an invalid BUS_TYPEext2int: bus 85 has an >>> invalid BUS_TYPEext2int: bus 86 has an invalid BUS_TYPEext2int: bus 87 has >>> an in >>> valid BUS_TYPEext2int: bus 88 has an invalid BUS_TYPEext2int: bus 89 has an >>> invalid BUS_TYPEext2int: bus 90 has an invalid BUS_TYPEext2int: bus 91 has >>> an i >>> nvalid BUS_TYPEext2int: bus 92 has an invalid BUS_TYPEext2int: bus 93 has >>> an invalid BUS_TYPEext2int: bus 94 has an invalid BUS_TYPEext2int: bus 95 >>> has an >>> invalid BUS_TYPEext2int: bus 96 has an invalid BUS_TYPEext2int: bus 97 has >>> an invalid BUS_TYPEext2int: bus 98 has an invalid BUS_TYPEext2int: bus 99 >>> has an >>> invalid BUS_TYPEext2int: bus 100 has an invalid BUS_TYPEext2int: bus 101 >>> has an invalid BUS_TYPEext2int: bus 102 has an invalid BUS_TYPEext2int: bus >>> 103 h >>> as an invalid BUS_TYPEext2int: bus 104 has an invalid BUS_TYPEext2int: bus >>> 105 has an invalid BUS_TYPEext2int: bus 106 has an invalid BUS_TYPEext2int: >>> bus >>> 107 has an invalid BUS_TYPEext2int: bus 108 has an invalid BUS_TYPEext2int: >>> bus 109 has an invalid BUS_TYPEext2int: bus 110 has an invalid >>> BUS_TYPEext2int: >>> bus 111 has an invalid BUS_TYPEext2int: bus 112 has an invalid >>> BUS_TYPEext2int: bus 113 has an invalid BUS_TYPEext2int: bus 114 has an >>> invalid BUS_TYPEext >>> 2int: bus 115 has an invalid BUS_TYPEext2int: bus 116 has an invalid >>> BUS_TYPEext2int: bus 117 has an invalid BUS_TYPEext2int: bus 118 has an >>> invalid BUS_TY >>> PEext2int: bus 119 has an invalid BUS_TYPEext2int: bus 120 has an invalid >>> BUS_TYPEext2int: bus 121 has an invalid BUS_TYPEext2int: bus 122 has an >>> invalid B >>> US_TYPEext2int: bus 123 has an invalid BUS_TYPEext2int: bus 124 has an >>> invalid BUS_TYPEext2int: bus 125 has an invalid BUS_TYPEext2int: bus 126 >>> has an inva >>> lid BUS_TYPEext2int: bus 127 has an invalid BUS_TYPEext2int: bus 128 has an >>> invalid BUS_TYPEext2int: bus 129 has an invalid BUS_TYPEext2int: bus 130 >>> has an >>> invalid BUS_TYPEext2int: bus 131 has an invalid BUS_TYPEext2int: bus 132 >>> has an invalid BUS_TYPEext2int: bus 133 has an invalid BUS_TYPEext2int: bus >>> 134 h >>> as an invalid BUS_TYPEext2int: bus 135 has an invalid BUS_TYPEext2int: bus >>> 136 has an invalid BUS_TYPEext2int: bus 137 has an invalid BUS_TYPEext2int: >>> bus >>> 138 has an invalid BUS_TYPEext2int: bus 139 has an invalid BUS_TYPEext2int: >>> bus 140 has an invalid BUS_TYPEext2int: bus 141 has an invalid >>> BUS_TYPEext2int: >>> bus 142 has an invalid BUS_TYPEext2int: bus 143 has an invalid >>> BUS_TYPEext2int: bus 144 has an invalid BUS_TYPEext2int: bus 145 has an >>> invalid BUS_TYPEext >>> 2int: bus 146 has an invalid BUS_TYPEext2int: bus 147 has an invalid >>> BUS_TYPEext2int: bus 148 has an invalid BUS_TYPEext2int: bus 149 has an >>> invalid BUS_TY >>> PEext2int: bus 150 has an invalid BUS_TYPEext2int: bus 151 has an invalid >>> BUS_TYPEext2int: bus 152 has an invalid BUS_TYPEext2int: bus 153 has an >>> invalid B >>> US_TYPEext2int: bus 154 has an invalid BUS_TYPEext2int: bus 155 has an >>> invalid BUS_TYPEext2int: bus 156 has an invalid BUS_TYPEext2int: bus 157 >>> has an inva >>> lid BUS_TYPEext2int: bus 158 has an invalid BUS_TYPEext2int: bus 159 has an >>> invalid BUS_TYPEext2int: bus 160 has an invalid BUS_TYPEext2int: bus 161 >>> has an >>> invalid BUS_TYPEext2int: bus 162 has an invalid BUS_TYPEext2int: bus 163 >>> has an invalid BUS_TYPEext2int: bus 164 has an invalid BUS_TYPEext2int: bus >>> 165 h >>> as an invalid BUS_TYPEext2int: bus 166 has an invalid BUS_TYPEext2int: bus >>> 167 has an invalid BUS_TYPEext2int: bus 168 has an invalid BUS_TYPEext2int: >>> bus >>> 169 has an invalid BUS_TYPEext2int: bus 170 has an invalid BUS_TYPEext2int: >>> bus 171 has an invalid BUS_TYPEext2int: bus 172 has an invalid >>> BUS_TYPEext2int: >>> bus 173 has an invalid BUS_TYPEext2int: bus 174 has an invalid >>> BUS_TYPEext2int: bus 175 has an invalid BUS_TYPEext2int: bus 176 has an >>> invalid BUS_TYPEext >>> 2int: bus 177 has an invalid BUS_TYPEext2int: bus 178 has an invalid >>> BUS_TYPEext2int: bus 179 has an invalid BUS_TYPEext2int: bus 180 has an >>> invalid BUS_TY >>> PEext2int: bus 181 has an invalid BUS_TYPEext2int: bus 182 has an invalid >>> BUS_TYPEext2int: bus 183 has an invalid BUS_TYPEext2int: bus 184 has an >>> invalid B >>> US_TYPEext2int: bus 185 has an invalid BUS_TYPEext2int: bus 186 has an >>> invalid BUS_TYPEext2int: bus 187 has an invalid BUS_TYPEext2int: bus 188 >>> has an inva >>> lid BUS_TYPEext2int: bus 189 has an invalid BUS_TYPEext2int: bus 190 has an >>> invalid BUS_TYPEext2int: bus 191 has an invalid BUS_TYPEext2int: bus 192 >>> has an >>> invalid BUS_TYPEext2int: bus 193 has an invalid BUS_TYPEext2int: bus 194 >>> has an invalid BUS_TYPEext2int: bus 195 has an invalid BUS_TYPEext2int: bus >>> 196 h >>> as an invalid BUS_TYPEext2int: bus 197 has an invalid BUS_TYPEext2int: bus >>> 198 has an invalid BUS_TYPEext2int: bus 199 has an invalid BUS_TYPEext2int: >>> bus >>> 200 has an invalid BUS_TYPEext2int: bus 201 has an invalid BUS_TYPEext2int: >>> bus 202 has an invalid BUS_TYPEext2int: bus 203 has an invalid >>> BUS_TYPEext2int: >>> bus 204 has an invalid BUS_TYPEext2int: bus 205 has an invalid >>> BUS_TYPEext2int: bus 206 has an invalid BUS_TYPEext2int: bus 207 has an >>> invalid BUS_TYPEext >>> 2int: bus 208 has an invalid BUS_TYPE >>> error: called from >>> ext2int at line 135 column 17 >>> >>> >>> >>> >>> >>> 14 Gen Matpower Model >>> >>> >>> function mpc = case14gen >>> %CASE39 Power flow data for 14 Gen Australian equivalent system. >>> % Please see CASEFORMAT for details on the case file format. >>> >>> >>> % MATPOWER >>> >>> %% MATPOWER Case Format : Version 2 >>> mpc.version = '2'; >>> >>> %%----- Power Flow Data -----%% >>> %% system MVA base >>> mpc.baseMVA = 100; >>> >>> %% bus data >>> % bus_i >>> type >>> Pd Qd >>> Gs Bs >>> area Vm >>> Va baseKV >>> zone Vmax >>> Vmin >>> mpc.bus = [ >>> 101 3 >>> 0 >>> 0 0 >>> 0 1 >>> 1 >>> 0 15 >>> 1 1.06 >>> 0.94; >>> 102 1 >>> 450 >>> 45 0 >>> 0 1 >>> 1.039 -1.4 >>> 330 1 >>> 1.06 >>> 0.94; >>> 201 2 >>> 0 >>> 0 0 >>> 0 2 >>> 1 >>> 48.9 20 >>> 1 1.06 >>> 0.94; >>> 202 2 >>> 0 >>> 0 0 >>> 0 2 >>> 1 >>> 37.9 20 >>> 1 1.06 >>> 0.94; >>> 203 2 >>> 0 >>> 0 0 >>> 0 2 >>> 1 >>> 32.24 20 >>> 1 1.06 >>> 0.94; >>> 204 2 >>> 0 >>> 0 0 >>> 0 2 >>> 1 >>> 39.21 20 >>> 1 1.06 >>> 0.94; >>> 205 2 >>> 390 >>> 39 0 >>> 0 2 >>> 1.055 44.02 >>> 330 1 >>> 1.06 >>> 0.94; >>> 206 1 >>> 130 >>> 13 0 >>> 0 2 >>> 1.0472 41.45 >>> 330 1 >>> 1.06 >>> 0.94; >>> 207 1 >>> 1880 >>> 188 0 >>> 0 2 >>> 1.0221 27.69 >>> 330 1 >>> 1.06 >>> 0.94; >>> 208 1 >>> 210 >>> 21 0 >>> 0 2 >>> 1.0324 26.25 >>> 330 1 >>> 1.06 >>> 0.94; >>> 209 1 >>> 0 >>> 0 0 >>> 0 2 >>> 1.0379 >>> 30.41 330 >>> 1 1.06 >>> 0.94; >>> 210 1 >>> 0 >>> 0 0 >>> 0 2 >>> 1.0594 >>> 26.68 500 >>> 1 1.06 >>> 0.94; >>> 211 1 >>> 1700 >>> 170 0 >>> 0 2 >>> 1.0081 18.87 >>> 330 1 >>> 1.06 >>> 0.94; >>> 212 1 >>> 1660 >>> 166 0 >>> 0 2 >>> 1.0084 19.13 >>> 330 1 >>> 1.06 >>> 0.94; >>> 213 1 >>> 0 >>> 0 0 >>> 0 2 >>> 1.0426 >>> 22.67 500 >>> 1 1.06 >>> 0.94; >>> 214 1 >>> 0 >>> 0 0 >>> 0 2 >>> 1.0257 >>> 18.75 330 >>> 1 1.06 >>> 0.94; >>> 215 1 >>> 480 >>> 48 0 >>> 0 2 >>> 1.0433 33.11 >>> 330 1 >>> 1.06 >>> 0.94; >>> 216 1 >>> 1840 >>> 184 0 >>> 0 2 >>> 1.013 16.1 >>> 330 1 >>> 1.06 >>> 0.94; >>> 217 1 >>> 1260 >>> 126 0 >>> 0 2 >>> 1.0021 9.39 >>> 330 1 >>> 1.06 >>> 0.94; >>> 301 2 >>> 0 >>> 0 0 >>> 0 3 >>> 1 >>> -3.39 20 >>> 1 1.06 >>> 0.94; >>> 302 2 >>> 0 >>> 0 0 >>> 0 3 >>> 1 >>> -18.02 20 >>> 1 1.06 >>> 0.94; >>> 303 1 >>> 0 >>> 0 0 >>> 0 3 >>> 1.0417 >>> -10.85 500 >>> 1 1.06 >>> 0.94; >>> 304 1 >>> 0 >>> 0 0 >>> 0 3 >>> 1.0093 >>> -21.75 500 >>> 1 1.06 >>> 0.94; >>> 305 1 >>> 0 >>> 0 0 >>> 0 3 >>> 1.0164 >>> -22.76 500 >>> 1 1.06 >>> 0.94; >>> 306 1 >>> 1230 >>> 123 0 >>> 0 3 >>> 1.016 -24.92 >>> 500 1 >>> 1.06 >>> 0.94; >>> 307 1 >>> 650 >>> 65 0 >>> 0 3 >>> 1.0182 -24.97 >>> 500 1 >>> 1.06 >>> 0.94; >>> 308 1 >>> 655 >>> 66 0 >>> 0 3 >>> 1.0359 -35.41 >>> 500 1 >>> 1.06 >>> 0.94; >>> 309 1 >>> 195 >>> 20 0 >>> 0 3 >>> 1.0291 -9.34 >>> 330 1 >>> 1.06 >>> 0.94; >>> 310 1 >>> 0 >>> 0 0 >>> 0 3 >>> 1.0198 >>> -24.81 330 >>> 1 1.06 >>> 0.94; >>> 311 1 >>> 0 >>> 0 0 >>> 0 3 >>> 1.0136 >>> -17.66 330 >>> 1 1.06 >>> 0.94; >>> 312 1 >>> 115 >>> 12 0 >>> 0 3 >>> 1.0383 -23.59 >>> 220 1 >>> 1.06 >>> 0.94; >>> 313 2 >>> 2405 >>> 240 0 >>> 0 3 >>> 1.015 -30.23 >>> 220 1 >>> 1.06 >>> 0.94; >>> 314 1 >>> 250 >>> 25 0 >>> 0 3 >>> 1.0158 -28.52 >>> 220 1 >>> 1.06 >>> 0.94; >>> 315 1 >>> 0 >>> 0 0 >>> 0 3 >>> 1.0376 >>> -39.01 275 >>> 1 1.06 >>> 0.94; >>> 401 2 >>> 0 >>> 0 0 >>> 0 4 >>> 1 >>> 74.34 20 >>> 1 1.06 >>> 0.94; >>> 402 2 >>> 0 >>> 0 0 >>> 0 4 >>> 1 >>> 107.94 20 >>> 1 1.06 >>> 0.94; >>> 403 2 >>> 0 >>> 0 0 >>> 0 4 >>> 1 >>> 113.34 20 >>> 1 1.06 >>> 0.94; >>> 404 2 >>> 0 >>> 0 0 >>> 0 4 >>> 1 >>> 106.53 20 >>> 1 1.06 >>> 0.94; >>> 405 1 >>> 989.9999 >>> 99 0 >>> 0 4 >>> 1.0312 99.54 >>> 275 1 >>> 1.06 >>> 0.94; >>> 406 1 >>> 740 >>> 74 0 >>> 0 4 >>> 1.0342 104.79 >>> 275 1 >>> 1.06 >>> 0.94; >>> 407 1 >>> 0 >>> 0 0 >>> 0 4 >>> 1.0393 >>> 107.11 275 >>> 1 1.06 >>> 0.94; >>> 408 1 >>> 150 >>> 15 0 >>> 0 4 >>> 1.0301 100.38 >>> 275 1 >>> 1.06 >>> 0.94; >>> 409 1 >>> 260 >>> 26 0 >>> 0 4 >>> 0.9855 76.61 >>> 275 1 >>> 1.06 >>> 0.94; >>> 410 1 >>> 530 >>> 53 0 >>> 0 4 >>> 1.0301 68.14 >>> 275 1 >>> 1.06 >>> 0.94; >>> 411 1 >>> 575 >>> 58 0 >>> 0 4 >>> 0.993 62.25 >>> 275 1 >>> 1.06 >>> 0.94; >>> 412 2 >>> 1255 >>> 126 0 >>> 0 4 >>> 1 60.97 >>> 275 1 >>> 1.06 >>> 0.94; >>> 413 1 >>> 0 >>> 0 0 >>> 0 4 >>> 1.0398 >>> 60.95 275 >>> 1 1.06 >>> 0.94; >>> 414 1 >>> 0 >>> 0 0 >>> 0 4 >>> 1.0418 >>> 60.2 330 >>> 1 1.06 >>> 0.94; >>> 415 1 >>> 0 >>> 0 0 >>> 0 4 >>> 1.0474 >>> 56.66 330 >>> 1 1.06 >>> 0.94; >>> 416 1 >>> 0 >>> 0 0 >>> 0 4 >>> 1.0552 >>> 50.27 330 >>> 1 1.06 >>> 0.94; >>> 501 2 >>> 0 >>> 0 0 >>> 0 5 >>> 1 >>> -55.32 20 >>> 1 1.06 >>> 0.94; >>> 502 2 >>> 0 >>> 0 0 >>> 0 5 >>> 1 >>> -53.95 15 >>> 1 1.06 >>> 0.94; >>> 503 2 >>> 0 >>> 0 0 >>> 0 5 >>> 1 >>> -56.28 15 >>> 1 1.06 >>> 0.94; >>> 504 1 >>> 300 >>> 60 0 >>> 0 5 >>> 1.0478 -63.31 >>> 275 1 >>> 1.06 >>> 0.94; >>> 505 1 >>> 0 >>> 0 0 >>> 0 5 >>> 1.0221 >>> -60.87 275 >>> 1 1.06 >>> 0.94; >>> 506 1 >>> 0 >>> 0 0 >>> 0 5 >>> 1.0204 >>> -62.18 275 >>> 1 1.06 >>> 0.94; >>> 507 2 >>> 1000 >>> 200 0 >>> 0 5 >>> 1.015 -63.45 >>> 275 1 >>> 1.06 >>> 0.94; >>> 508 1 >>> 800 >>> 160 0 >>> 0 5 >>> 1.0104 -64.66 >>> 275 1 >>> 1.06 >>> 0.94; >>> 509 2 >>> 200 >>> 40 0 >>> 0 5 >>> 1.03 -45.74 >>> 275 1 >>> 1.06 >>> 0.94; >>> >>> >>> ]; >>> >>> %% generator data >>> % bus >>> Pg >>> Qg Qmax >>> Qmin Vg >>> mBase status >>> Pmax Pmin >>> Pc1 Pc2 >>> Qc1min Qc1max >>> Qc2min Qc2max >>> ramp_agc ramp_10 >>> ramp_30 ramp_q >>> apf >>> mpc.gen = [ >>> 101 300.8091 >>> 311.4755 581.128 >>> -581.128 1 >>> 100 1 >>> 1200 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 201 3600 >>> 573.6404 1743.648 >>> -1743.648 1 >>> 100 1 >>> 3600 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 202 2500 >>> 663.3815 1210.9 >>> -1210.9 1 >>> 100 1 >>> 2500 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 203 1500 >>> 531.1833 968.72 >>> -968.72 1 >>> 100 1 >>> 2000 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 204 2950.2 >>> 734.1567 1743.648 >>> -1743.648 1 >>> 100 1 >>> 3600 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 301 4200 >>> 996.5973 2034.256 >>> -2034.256 1 >>> 100 1 >>> 4200 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 302 939.9 >>> 154.5849 581.127 >>> -581.127 1 >>> 100 1 >>> 1200 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 401 1400 >>> 514.8802 774.836 >>> -774.836 1 >>> 100 1 >>> 1600 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 402 837 >>> 177.8523 435.846 >>> -435.846 1 >>> 100 1 >>> 900 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0; >>> 403 1400 >>> 209.1907 774.836 >>> -774.836 1 >>> 100 1 >>> 1600 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 404 1549.8 >>> 326.8041 871.692 >>> -871.692 1 >>> 100 1 >>> 1800 >>> 0 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0; >>> 501 600 >>> 50.694 290.564 >>> -290.564 1 >>> 100 1 >>> 600 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0; >>> 502 800 >>> 160.2667 600 >>> -600 1 >>> 100 >>> 1 800 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0; >>> 503 436 >>> 100.8554 290.652 >>> -290.652 1 >>> 100 1 >>> 600 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0; >>> 205 0 >>> -68.2624 >>> 430 -220 >>> 1.055 100 >>> 1 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0; >>> 313 0 >>> 71.5048 >>> 600 -200 >>> 1.015 100 >>> 1 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0; >>> 412 0 >>> 58.1613 >>> 1100 -330 >>> 1 100 >>> 1 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0; >>> 507 0 >>> 22.6476 >>> 320 -180 >>> 1.015 100 >>> 1 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0; >>> 509 0 >>> 10.569 >>> 400 -150 >>> 1.03 100 >>> 1 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0 >>> 0 >>> 0 0 >>> 0 0; >>> >>> >>> ]; >>> >>> %% branch data >>> % fbus >>> tbus >>> r x >>> b rateA >>> rateB rateC >>> ratio angle >>> status angmin >>> angmax >>> mpc.branch = [ >>> 102 217 >>> 0.0084 0.0667 >>> 0.817 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 102 217 >>> 0.0084 0.0667 >>> 0.817 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 102 217 >>> 0.0078 0.062 >>> 0.76 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 102 217 >>> 0.0078 0.062 >>> 0.76 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 102 309 >>> 0.0045 0.0356 >>> 0.437 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 102 309 >>> 0.0045 0.0356 >>> 0.437 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 102 309 >>> 0.0109 0.0868 >>> 0.76 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 205 206 >>> 0.0096 0.076 >>> 0.931 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 205 206 >>> 0.0096 0.076 >>> 0.931 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 205 416 >>> 0.0037 0.046 >>> 0.73 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 205 416 >>> 0.0037 0.046 >>> 0.73 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 206 207 >>> 0.0045 0.0356 >>> 0.437 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 206 207 >>> 0.0045 0.0356 >>> 0.437 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 206 212 >>> 0.0066 0.0527 >>> 0.646 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 206 212 >>> 0.0066 0.0527 >>> 0.646 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 206 215 >>> 0.0066 0.0527 >>> 0.646 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 206 215 >>> 0.0066 0.0527 >>> 0.646 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 207 208 >>> 0.0018 0.014 >>> 0.171 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 207 208 >>> 0.0018 0.014 >>> 0.171 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 207 209 >>> 0.0008 0.0062 >>> 0.076 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 208 211 >>> 0.0031 0.0248 >>> 0.304 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 208 211 >>> 0.0031 0.0248 >>> 0.304 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 208 211 >>> 0.0031 0.0248 >>> 0.304 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 209 212 >>> 0.0045 0.0356 >>> 0.437 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 210 213 >>> 0.001 0.0145 >>> 1.54 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 210 213 >>> 0.001 0.0145 >>> 1.54 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 211 212 >>> 0.0014 0.0108 >>> 0.133 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 211 212 >>> 0.0014 0.0108 >>> 0.133 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 211 214 >>> 0.0019 0.0155 >>> 0.19 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 212 217 >>> 0.007 0.0558 >>> 0.684 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 214 216 >>> 0.001 0.0077 >>> 0.095 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 214 217 >>> 0.0049 0.0388 >>> 0.475 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 215 216 >>> 0.0051 0.0403 >>> 0.494 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 215 216 >>> 0.0051 0.0403 >>> 0.494 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 215 217 >>> 0.0072 0.0574 >>> 0.703 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 215 217 >>> 0.0072 0.0574 >>> 0.703 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 216 217 >>> 0.0051 0.0403 >>> 0.494 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 303 304 >>> 0.002 0.028 >>> 0.74 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 303 304 >>> 0.002 0.028 >>> 0.74 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 303 305 >>> 0.0011 0.016 >>> 1.7 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 303 305 >>> 0.0011 0.016 >>> 1.7 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 304 305 >>> 0.0003 0.004 >>> 0.424 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 305 306 >>> 0.0002 0.003 >>> 0.32 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 305 307 >>> 0.0003 0.0045 >>> 0.447 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 305 307 >>> 0.0003 0.0045 >>> 0.447 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 306 307 >>> 0.0001 0.0012 >>> 0.127 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 307 308 >>> 0.0023 0.0325 >>> 3.445 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 307 308 >>> 0.0023 0.0325 >>> 3.445 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 309 310 >>> 0.0135 0.10695 >>> 0.58267 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 309 310 >>> 0.0135 0.10695 >>> 0.58267 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 309 310 >>> 0.0135 0.10695 >>> 0.58267 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 310 311 >>> 0 -0.0337 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 310 311 >>> 0 -0.0337 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 312 313 >>> 0.006 0.045 >>> 0.3 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 312 313 >>> 0.006 0.045 >>> 0.3 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 312 313 >>> 0.006 0.045 >>> 0.3 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 313 314 >>> 0.001 0.01 >>> 0.26 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 313 314 >>> 0.001 0.01 >>> 0.26 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 315 509 >>> 0.007 0.05 >>> 0.19 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 315 509 >>> 0.007 0.05 >>> 0.19 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 405 406 >>> 0.0039 0.0475 >>> 0.381 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 405 406 >>> 0.0039 0.0475 >>> 0.381 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 405 408 >>> 0.0054 0.05 >>> 0.189 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 405 409 >>> 0.018 0.122 >>> 0.79 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 405 409 >>> 0.018 0.122 >>> 0.79 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 405 409 >>> 0.018 0.122 >>> 0.79 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 406 407 >>> 0.0006 0.0076 >>> 0.062 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 406 407 >>> 0.0006 0.0076 >>> 0.062 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 407 408 >>> 0.0042 0.0513 >>> 0.412 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 408 410 >>> 0.0165 0.192 >>> 0.67333 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 408 410 >>> 0.0165 0.192 >>> 0.67333 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 408 410 >>> 0.0165 0.192 >>> 0.67333 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 409 411 >>> 0.0103 0.0709 >>> 0.46 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 409 411 >>> 0.0103 0.0709 >>> 0.46 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 410 411 >>> 0.0043 0.0532 >>> 0.427 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 410 412 >>> 0.0043 0.0532 >>> 0.427 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 410 412 >>> 0.0043 0.0532 >>> 0.427 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 410 412 >>> 0.0043 0.0532 >>> 0.427 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 410 412 >>> 0.0043 0.0532 >>> 0.427 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 410 413 >>> 0.004 0.0494 >>> 0.4 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 410 413 >>> 0.004 0.0494 >>> 0.4 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 411 412 >>> 0.0012 0.0152 >>> 0.122 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 411 412 >>> 0.0012 0.0152 >>> 0.122 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 414 415 >>> 0.002 0.025 >>> 0.39 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 414 415 >>> 0.002 0.025 >>> 0.39 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 415 416 >>> 0.0037 0.046 >>> 0.73 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 415 416 >>> 0.0037 0.046 >>> 0.73 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 504 507 >>> 0.023 0.15 >>> 0.56 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 504 507 >>> 0.023 0.15 >>> 0.56 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 504 508 >>> 0.026 0.019 >>> 0.87 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 504 508 >>> 0.026 0.019 >>> 0.87 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 505 507 >>> 0.0016 0.017 >>> 0.03 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 505 507 >>> 0.0016 0.017 >>> 0.03 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 505 508 >>> 0.0025 0.028 >>> 0.17 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 506 507 >>> 0.0016 0.017 >>> 0.03 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 506 507 >>> 0.0016 0.017 >>> 0.03 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 506 508 >>> 0.003 0.028 >>> 0.14 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 507 508 >>> 0.002 0.019 >>> 0.09 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 507 509 >>> 0.09 0.66 >>> 0.3 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 507 509 >>> 0.09 0.66 >>> 0.3 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 507 509 >>> 0.09 0.66 >>> 0.3 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 507 509 >>> 0.09 0.66 >>> 0.3 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 507 509 >>> 0.09 0.66 >>> 0.3 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 507 509 >>> 0.09 0.66 >>> 0.3 9999 >>> 9999 9999 >>> 0 0 >>> 1 >>> -360 360; >>> 101 102 >>> 0 0.009 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 201 206 >>> 0 0.004 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 202 209 >>> 0 0.00576 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 203 208 >>> 0 0.00765 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 204 215 >>> 0 0.004 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 209 210 >>> 0 0.00000503 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 213 214 >>> 0 0.00000503 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 301 303 >>> 0 0.00343 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 302 312 >>> 0 0.01127 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 304 313 >>> 0 0.016 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 305 311 >>> 0 0.00048600 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 305 314 >>> 0 0.012 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 308 315 >>> 0 0.0135 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 401 410 >>> 0 0.00845 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 402 408 >>> 0 0.017 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 403 407 >>> 0 0.00845 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 404 405 >>> 0 0.0085 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 413 414 >>> 0 0.00002133 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 501 504 >>> 0 0.0255 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 502 505 >>> 0 0.016 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> 503 506 >>> 0 0.025 >>> 0 9999 >>> 9999 >>> 9999 0 >>> 0 1 >>> -360 360; >>> >>> >>> ]; >>> >>> %%----- OPF Data -----%% >>> %% generator cost data >>> % 1 >>> startup >>> shutdown n >>> x1 y1 >>> ... xn >>> yn >>> % 2 >>> startup >>> shutdown n >>> c(n-1) ... >>> c0 >>> mpc.gencost = [ >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> 2 0 >>> 0 >>> 3 0.01 >>> 0.3 0.2; >>> >>> ];
