My MATLAB program successfully run for IEEE14,30and57 system and give result but when the same program is applied for 75 bus system which is a Indian system then it shows warning error(matrix is close to singular or badly scaled,result may be inaccurate) and does not stop . Please give me any suggestions regarding this matter.
% Returns Initial Bus datas of the system...
function busdt = busdatas(num)
% Type....
% 1 - Slack Bus..
% 2 - PV Bus..
% 3 - PQ Bus..
% |Bus | Type | Vsp | theta | PGi | QGi | PLi | QLi | Qmin | Qmax |
busdat14 = [1 1 1.060 0 0 0 0 0 0 0;
2 2 1.045 0 40 42.4 21.7 12.7 -40 50;
3 2 1.010 0 0 23.4 94.2 19.0 0 40;
4 3 1.0 0 0 0 47.8 -3.9 0 0;
5 3 1.0 0 0 0 7.6 1.6 0 0;
6 2 1.070 0 0 12.2 11.2 7.5 -6 24;
7 3 1.0 0 0 0 0.0 0.0 0 0;
8 2 1.090 0 0 17.4 0.0 0.0 -6 24;
9 3 1.0 0 0 0 29.5 16.6 0 0;
10 3 1.0 0 0 0 9.0 5.8 0 0;
11 3 1.0 0 0 0 3.5 1.8 0 0;
12 3 1.0 0 0 0 6.1 1.6 0 0;
13 3 1.0 0 0 0 13.5 5.8 0 0;
14 3 1.0 0 0 0 14.9 5.0 0 0;];
% |Bus | Type | Vsp | theta | PGi | QGi | PLi | QLi | Qmin | Qmax |
busdat30 = [1 1 1.06 0 0 0 0 0 0 0;
2 2 1.043 0 40 50.0 21.7 12.7 -40 50;
3 3 1.0 0 0 0 2.4 1.2 0 0;
4 3 1.06 0 0 0 7.6 1.6 0 0;
5 2 1.01 0 0 37.0 94.2 19.0 -40 40;
6 3 1.0 0 0 0 0.0 0.0 0 0;
7 3 1.0 0 0 0 22.8 10.9 0 0;
8 2 1.01 0 0 37.3 30.0 30.0 -10 40;
9 3 1.0 0 0 0 0.0 0.0 0 0;
10 3 1.0 0 0 19.0 5.8 2.0 0 0;
11 2 1.082 0 0 16.2 0.0 0.0 -6 24;
12 3 1.0 0 0 0 11.2 7.5 0 0;
13 2 1.071 0 0 10.6 0.0 0.0 -6 24;
14 3 1.0 0 0 0 6.2 1.6 0 0;
15 3 1.0 0 0 0 8.2 2.5 0 0;
16 3 1.0 0 0 0 3.5 1.8 0 0;
17 3 1.0 0 0 0 9.0 5.8 0 0;
18 3 1.0 0 0 0 3.2 0.9 0 0;
19 3 1.0 0 0 0 9.5 3.4 0 0;
20 3 1.0 0 0 0 2.2 0.7 0 0;
21 3 1.0 0 0 0 17.5 11.2 0 0;
22 3 1.0 0 0 0 0.0 0.0 0 0;
23 3 1.0 0 0 0 3.2 1.6 0 0;
24 3 1.0 0 0 4.3 8.7 6.7 0 0;
25 3 1.0 0 0 0 0.0 0.0 0 0;
26 3 1.0 0 0 0 3.5 2.3 0 0;
27 3 1.0 0 0 0 0.0 0.0 0 0;
28 3 1.0 0 0 0 0.0 0.0 0 0;
29 3 1.0 0 0 0 2.4 0.9 0 0;
30 3 1.0 0 0 0 10.6 1.9 0 0 ];
% |Bus | Type | Vsp | theta | PGi | QGi | PLi | QLi | Qmin | Qmax |
busdat57 = [1 1 1.040 0 0.0 0.0 0.0 0.0 0.0 0.0;
2 2 1.010 0 3.0 88.0 0.0 -0.8 50.0 -17.0;
3 2 0.985 0 41.0 21.0 40.0 -1.0 60.0 -10.0;
4 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
5 3 1.000 0 13.0 4.0 0.0 0.0 0.0 0.0;
6 2 0.980 0 75.0 2.0 0.0 0.8 25.0 -8.0;
7 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
8 2 1.005 0 150.0 22.0 450.0 62.1 200.0 -140.0;
9 2 0.980 0 121.0 26.0 0.0 2.2 9.0 -3.0;
10 3 1.000 0 5.0 2.0 0.0 0.0 0.0 0.0;
11 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
12 2 1.015 0 377.0 24.0 310.0 128.5 155.0 -150.0;
13 3 1.000 0 18.0 2.3 0.0 0.0 0.0 0.0;
14 3 1.000 0 10.5 5.3 0.0 0.0 0.0 0.0;
15 3 1.000 0 22.0 5.0 0.0 0.0 0.0 0.0;
16 3 1.000 0 43.0 3.0 0.0 0.0 0.0 0.0;
17 3 1.000 0 42.0 8.0 0.0 0.0 0.0 0.0;
18 3 1.000 0 27.2 9.8 0.0 0.0 0.0 0.0;
19 3 1.000 0 3.3 0.6 0.0 0.0 0.0 0.0;
20 3 1.000 0 2.3 1.0 0.0 0.0 0.0 0.0;
21 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
22 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
23 3 1.000 0 6.3 2.1 0.0 0.0 0.0 0.0;
24 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
25 3 1.000 0 6.3 3.2 0.0 0.0 0.0 0.0;
26 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
27 3 1.000 0 9.3 0.5 0.0 0.0 0.0 0.0;
28 3 1.000 0 4.6 2.3 0.0 0.0 0.0 0.0;
29 3 1.000 0 17.0 2.6 0.0 0.0 0.0 0.0;
30 3 1.000 0 3.6 1.8 0.0 0.0 0.0 0.0;
31 3 1.000 0 5.8 2.9 0.0 0.0 0.0 0.0;
32 3 1.000 0 1.6 0.8 0.0 0.0 0.0 0.0;
33 3 1.000 0 3.8 1.9 0.0 0.0 0.0 0.0;
34 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
35 3 1.000 0 6.0 3.0 0.0 0.0 0.0 0.0;
36 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
37 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
38 3 1.000 0 14.0 7.0 0.0 0.0 0.0 0.0;
39 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
40 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
41 3 1.000 0 6.3 3.0 0.0 0.0 0.0 0.0;
42 3 1.000 0 7.1 4.4 0.0 0.0 0.0 0.0;
43 3 1.000 0 2.0 1.0 0.0 0.0 0.0 0.0;
44 3 1.000 0 12.0 1.8 0.0 0.0 0.0 0.0;
45 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
46 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
47 3 1.000 0 29.7 11.6 0.0 0.0 0.0 0.0;
48 3 1.000 0 0.0 0.0 0.0 0.0 0.0 0.0;
49 3 1.000 0 18.0 8.5 0.0 0.0 0.0 0.0;
50 3 1.000 0 21.0 10.5 0.0 0.0 0.0 0.0;
51 3 1.000 0 18.0 5.3 0.0 0.0 0.0 0.0;
52 3 1.000 0 4.9 2.2 0.0 0.0 0.0 0.0;
53 3 1.000 0 20.0 10.0 0.0 0.0 0.0 0.0;
54 3 1.000 0 4.1 1.4 0.0 0.0 0.0 0.0;
55 3 1.000 0 6.8 3.4 0.0 0.0 0.0 0.0;
56 3 1.000 0 7.6 2.2 0.0 0.0 0.0 0.0;
57 3 1.000 0 6.7 2.0 0.0 0.0 0.0
0.0];
busdat75=[ 1 1 1.03 0 7.25 0 0 0 0
4
2 2 1.05 0 2.6 0 0 0 0 0.96
3 2 1.03 0 1.8 0 0 0 0 0.83
4 2 1.05 0 1 0 0 0 0 0.6
5 2 1.05 0 1.8 0 0 0 0 0.31
6 2 1.05 0 1.2 0 0 0 0 0.2
7 2 1.05 0 0.6 0 0 0 0 0.19
8 2 1.05 0 0.8 0 0 0 0 0.68
9 2 1.05 0 5.5 0 0 0 0 2.5
10 2 1.02 0 0.8 0 0 0 0 0.56
11 2 1.02 0 1.09 0 0 0 0 1.05
12 2 1.05 0 18 0 0 0 0 3.44
13 2 1.05 0 9 0 0 0 0 2.8
14 2 1.03 0 1.5 0 0 0 0 0.84
15 2 1.01 0 4.54 0 0 0 -0.3 0.35
16 3 1 0 0 0 -0.5869 0.2757 0 0
17 3 1 0 0 0 0 0 0 0
18 3 1 0 0 0 0 0 0 0
19 3 1 0 0 0 0 0 0 0
20 3 1 0 0 0 1.5637 0.3393 0 0
21 3 1 0 0 0 0 0 0 0
22 3 1 0 0 0 0 0 0 0
23 3 1 0 0 0 0 0 0 0
24 3 1 0 0 0 2.2795 0.4053 0 0
25 3 1 0 0 0 2.1048 0.4343 0 0
26 3 1 0 0 0 0 0 0 0
27 3 1 0 0 0 3.06 0.407 0 0
28 3 1 0 0 0 1.2775 0.2835 0 0
29 3 1 0 0 0 0 0 0 0
30 3 1 0 0 0 2.2646 0.4424 0 0
31 3 1 0 0 0 0 0 0 0
32 3 1 0 0 0 0.7811 0.1159 0 0
33 3 1 0 0 0 0 0 0 0
34 3 1 0 0 0 0.817 0.8384 0 0
35 3 1 0 0 0 0 0 0 0
36 3 1 0 0 0 0 0 0 0
37 3 1 0 0 0 1.4428 0.4093 0 0
38 3 1 0 0 0 0 0 0 0
39 3 1 0 0 0 0.8515 0.2946 0 0
40 3 1 0 0 0 0 0 0 0
41 3 1 0 0 0 0 0 0 0
42 3 1 0 0 0 10 0 0 0
43 3 1 0 0 0 0 0 0 0
44 3 1 0 0 0 0 0 0 0
45 3 1 0 0 0 0 0 0 0
46 3 1 0 0 0 1.5632 0.8336 0 0
47 3 1 0 0 0 0.7455 0.0877 0 0
48 3 1 0 0 0 0.5083 0.1397 0 0
49 3 1 0 0 0 0.8572 0.5158 0 0
50 3 1 0 0 0 2.021 0.7494 0 0
51 3 1 0 0 0 0.5776 0.0062 0 0
52 3 1 0 0 0 0.9827 -0.1815 0 0
53 3 1 0 0 0 0.8263 -0.33 0 0
54 3 1 0 0 0 1.6195 0.29 0 0
55 3 1 0 0 0 2.7423 0.3108 0 0
56 3 1 0 0 0 1.4128 0.3432 0 0
57 3 1 0 0 0 1.5278 0.1853 0 0
58 3 1 0 0 0 0.9419 0.1129 0 0
59 3 1 0 0 0 1.2189 0.1101 0 0
60 3 1 0 0 0 0.742 0.0742 0 0
61 3 1 0 0 0 1.065 0.0658 0 0
62 3 1 0 0 0 1.0018 0.1813 0 0
63 3 1 0 0 0 0.5801 0.0531 0 0
64 3 1 0 0 0 0.5679 0.1333 0 0
65 3 1 0 0 0 1.4784 0.1281 0 0
66 3 1 0 0 0 0.3174 0.1518 0 0
67 3 1 0 0 0 0.9643 0.1055 0 0
68 3 1 0 0 0 0.4287 0.336 0 0
69 3 1 0 0 0 0.5594 0.3253 0 0
70 3 1 0 0 0 0.2334 0.023 0 0
71 3 1 0 0 0 0.9173 0.2136 0 0
72 3 1 0 0 0 0.5252 0.1176 0 0
73 3 1 0 0 0 4.47 0 0 0
74 3 1 0 0 0 2.88 0 0 0
75 3 1 0 0 0 -4.44 0 0 0];
switch num
case 14
busdt = busdat14;
case 30
busdt = busdat30;
case 57
busdt = busdat57;
case 75
busdt=busdat75;
end% Returns Line datas of the system...
function linedt = linedatas(num)
% | From | To | R | X | B/2 | X'mer |
% | Bus | Bus | pu | pu | pu | TAP (a) |
linedat14 = [1 2 0.01938 0.05917 0.0264 1
1 5 0.05403 0.22304 0.0246 1
2 3 0.04699 0.19797 0.0219 1
2 4 0.05811 0.17632 0.0170 1
2 5 0.05695 0.17388 0.0173 1
3 4 0.06701 0.17103 0.0064 1
4 5 0.01335 0.04211 0.0 1
4 7 0.0 0.20912 0.0 0.978
4 9 0.0 0.55618 0.0 0.969
5 6 0.0 0.25202 0.0 0.932
6 11 0.09498 0.19890 0.0 1
6 12 0.12291 0.25581 0.0 1
6 13 0.06615 0.13027 0.0 1
7 8 0.0 0.17615 0.0 1
7 9 0.0 0.11001 0.0 1
9 10 0.03181 0.08450 0.0 1
9 14 0.12711 0.27038 0.0 1
10 11 0.08205 0.19207 0.0 1
12 13 0.22092 0.19988 0.0 1
13 14 0.17093 0.34802 0.0 1 ];
% | From | To | R | X | B/2 | X'mer |
% | Bus | Bus | pu | pu | pu | TAP (a) |
linedat30 = [1 2 0.0192 0.0575 0.0264 1
1 3 0.0452 0.1652 0.0204 1
2 4 0.0570 0.1737 0.0184 1
3 4 0.0132 0.0379 0.0042 1
2 5 0.0472 0.1983 0.0209 1
2 6 0.0581 0.1763 0.0187 1
4 6 0.0119 0.0414 0.0045 1
5 7 0.0460 0.1160 0.0102 1
6 7 0.0267 0.0820 0.0085 1
6 8 0.0120 0.0420 0.0045 1
6 9 0.0 0.2080 0.0 0.978
6 10 0.0 0.5560 0.0 0.969
9 11 0.0 0.2080 0.0 1
9 10 0.0 0.1100 0.0 1
4 12 0.0 0.2560 0.0 0.932
12 13 0.0 0.1400 0.0 1
12 14 0.1231 0.2559 0.0 1
12 15 0.0662 0.1304 0.0 1
12 16 0.0945 0.1987 0.0 1
14 15 0.2210 0.1997 0.0 1
16 17 0.0824 0.1923 0.0 1
15 18 0.1073 0.2185 0.0 1
18 19 0.0639 0.1292 0.0 1
19 20 0.0340 0.0680 0.0 1
10 20 0.0936 0.2090 0.0 1
10 17 0.0324 0.0845 0.0 1
10 21 0.0348 0.0749 0.0 1
10 22 0.0727 0.1499 0.0 1
21 23 0.0116 0.0236 0.0 1
15 23 0.1000 0.2020 0.0 1
22 24 0.1150 0.1790 0.0 1
23 24 0.1320 0.2700 0.0 1
24 25 0.1885 0.3292 0.0 1
25 26 0.2544 0.3800 0.0 1
25 27 0.1093 0.2087 0.0 1
28 27 0.0 0.3960 0.0 0.968
27 29 0.2198 0.4153 0.0 1
27 30 0.3202 0.6027 0.0 1
29 30 0.2399 0.4533 0.0 1
8 28 0.0636 0.2000 0.0214 1
6 28 0.0169 0.0599 0.065 1 ];
% | From | To | R | X | B/2 | X'mer |
% | Bus | Bus | pu | pu | pu | TAP (a) |
linedat57 = [ 1 2 0.0083 0.0280 0.0645 1
2 3 0.0298 0.0850 0.0409 1
3 4 0.0112 0.0366 0.0190 1
4 5 0.0625 0.1320 0.0129 1
4 6 0.0430 0.1480 0.0174 1
6 7 0.0200 0.1020 0.0138 1
6 8 0.0339 0.1730 0.0235 1
8 9 0.0099 0.0505 0.0274 1
9 10 0.0369 0.1679 0.0220 1
9 11 0.0258 0.0848 0.0109 1
9 12 0.0648 0.2950 0.0386 1
9 13 0.0481 0.1580 0.0203 1
13 14 0.0132 0.0434 0.0055 1
13 15 0.0269 0.0869 0.0115 1
1 15 0.0178 0.0910 0.0494 1
1 16 0.0454 0.2060 0.0273 1
1 17 0.0238 0.1080 0.0143 1
3 15 0.0162 0.0530 0.0272 1
4 18 0.0 0.5550 0.0 0.970
4 18 0.0 0.4300 0.0 0.978
5 6 0.0302 0.0641 0.0062 1
7 8 0.0139 0.0712 0.0097 1
10 12 0.0277 0.1262 0.0164 1
11 13 0.0223 0.0732 0.0094 1
12 13 0.0178 0.0580 0.0302 1
12 16 0.0180 0.0813 0.0108 1
12 17 0.0397 0.1790 0.0238 1
14 15 0.0171 0.0547 0.0074 1
18 19 0.4610 0.6850 0.0 1
19 20 0.2830 0.4340 0.0 1
21 20 0.0 0.7767 0.0 1.043
21 22 0.0736 0.1170 0.0 1
22 23 0.0099 0.0152 0.0 1
23 24 0.1660 0.2560 0.0042 1
24 25 0.0 1.1820 0.0 1
24 25 0.0 1.2300 0.0 1
24 26 0.0 0.0473 0.0 1.043
26 27 0.1650 0.2540 0.0 1
27 28 0.0618 0.0954 0.0 1
28 29 0.0418 0.0587 0.0 1
7 29 0.0 0.0648 0.0 0.967
25 30 0.1350 0.2020 0.0 1
30 31 0.3260 0.4970 0.0 1
31 32 0.5070 0.7550 0.0 1
32 33 0.0392 0.0360 0.0 1
34 32 0.0 0.9530 0.0 0.975
34 35 0.0520 0.0780 0.0016 1
35 36 0.0430 0.0537 0.0008 1
36 37 0.0290 0.0366 0.0 1
37 38 0.0651 0.1009 0.0010 1
37 39 0.0239 0.0379 0.0 1
36 40 0.0300 0.0466 0.0 1
22 38 0.0192 0.0295 0.0 1
11 41 0.0 0.7490 0.0 0.955
41 42 0.2070 0.3520 0.0 1
41 43 0.0 0.4120 0.0 1
38 44 0.0289 0.0585 0.0010 1
15 45 0.0 0.1042 0.0 0.955
14 46 0.0 0.0735 0.0 0.900
46 47 0.0230 0.0680 0.0016 1
47 48 0.0182 0.0233 0.0 1
48 49 0.0834 0.1290 0.0024 1
49 50 0.0801 0.1280 0.0 1
50 51 0.1386 0.2200 0.0 1
10 51 0.0 0.0712 0.0 0.930
13 49 0.0 0.1910 0.0 0.895
29 52 0.1442 0.1870 0.0 1
52 53 0.0762 0.0984 0.0 1
53 54 0.1878 0.2320 0.0 1
54 55 0.1732 0.2265 0.0 1
11 43 0.0 0.1530 0.0 0.958
44 45 0.0624 0.1242 0.0020 1
40 56 0.0 1.1950 0.0 0.958
56 41 0.5530 0.5490 0.0 1
56 42 0.2125 0.3540 0.0 1
39 57 0.0 1.3550 0.0 0.980
57 56 0.1740 0.2600 0.0 1
38 49 0.1150 0.1770 0.0015 1
38 48 0.0312 0.0482 0.0 1
9 55 0.0 0.1205 0.0 0.940];
linedat75=[ 19 19 19 19 19 19
16 2 0.00123 0.02469 0 1
18 3 0 0.02917 0 1
17 16 0.00065 0.02604 0 1
17 1 0.00073 0.0146 0 1
22 25 0.00065 0.02604 0 1
23 24 0.00065 0.026 0 1
26 27 0.00065 0.02604 0 1
28 4 0.00306 0.06135 0 1
29 30 0.00043 0.01736 0 1
31 5 0.00235 0.0471 0 1
32 6 0.00514 0.10285 0 1
33 7 0.00549 0.10978 0 1
34 8 0 0.0486 0 1
35 9 0.00049 0.01943 0 1
36 37 0.00065 0.02604 0 1
38 39 0.0013 0.05208 0 1
24 10 0.00243 0.0486 0 1
40 11 0.0077 0.0272 0 1
41 12 0.00016 0.00591 0 1
42 13 0.0003 0.01199 0 1
43 14 0 0.02841 0 1
44 15 0 0.02273 0 1
45 44 0.00056 0.02222 0 1
16 46 0.0081 0.0388 0.2806 1
16 50 0.00993 0.04746 0.77286 1
17 19 0.00468 0.0477 1.259 1
17 23 0.00785 0.0799 2.09476 1
23 29 0.00806 0.08169 2.13616 1
20 64 0.0183 0.0927 0.1478 1
19 26 0.00294 0.02997 0.78412 1
47 50 0.01093 0.05221 0.37784 1
35 41 0.00031 0.0031 0.08112 1
47 67 0.00662 0.03164 0.22902 1
24 27 0.00505 0.02416 0.1746 1
24 54 0.01291 0.06171 0.44654 1
27 51 0.016 0.081 0.1288 1
51 52 0.0155 0.0794 0.126 1
25 60 0.0166 0.0843 0.1344 1
25 43 0.0127 0.0641 0.1044 1
34 54 0.0177 0.0851 0.6052 1
54 28 0.0106 0.0506 0.3664 1
28 43 0.0058 0.029 0.0474 1
28 56 0.0037 0.0178 0.0644 1
56 30 0.0049 0.0237 0.1718 1
30 57 0.0075 0.0384 0.0622 1
53 30 0.00679 0.03412 0.05564 1
53 61 0.00666 0.0339 0.05344 1
30 61 0.0144 0.0731 0.117 1
57 58 0.0067 0.0339 0.0534 1
57 59 0.00583 0.02956 0.04692 1
59 39 0.0141 0.0718 0.114 1
39 31 0.0144 0.0725 0.118 1
39 33 0.00705 0.0359 0.228 1
54 63 0.0099 0.0509 0.0802 1
55 63 0.0078 0.0398 0.0628 1
61 62 0.0116 0.0583 0.095 1
62 32 0.0069 0.035 0.2252 1
31 32 0.0005 0.00253 0.0161 1
35 36 0.00479 0.0488 1.27228 1
46 37 0.01732 0.08784 0.13964 1
19 36 0.00254 0.02584 0.67596 1
17 35 0.00051 0.00517 0.1352 1
20 40 0.0058 0.0294 0.1884 1
40 48 0.0083 0.0424 0.0668 1
74 41 0.00927 0.09429 0.46586 1
74 41 0.00833 0.08478 2.2171 1
74 73 0.00559 0.05686 1.48708 1
26 22 0.0065 0.06617 1.73042 1
29 22 0.0026 0.02646 0.6922 1
26 41 0.00823 0.08375 0.19106 1
48 49 0.0093 0.0475 0.0748 1
49 40 0.0133 0.0668 0.1084 1
38 29 0.0037 0.03762 0.9774 1
38 22 0.00325 0.03307 0.86528 1
18 47 0.00437 0.02552 0.18798 1
30 65 0.00248 0.01186 0.08588 1
41 42 0.00031 0.0031 0.08112 1
42 74 0.00918 0.09306 2.4336 1
20 66 0.01325 0.06667 0.10832 1
23 74 0.00015 0.00155 0.05408 1
24 67 0.00124 0.00593 0.04294 1
18 68 0.00336 0.01963 0.03716 1
18 71 0.01344 0.07852 0.1446 1
27 68 0.01344 0.07852 0.1446 1
27 71 0.00336 0.01963 0.03616 1
25 72 0.01598 0.08108 0.12872 1
43 58 0.01315 0.06696 0.11556 1
39 31 0.0144 0.0725 0.118 1
43 56 0.00499 0.02397 0.17046 1
55 44 0.00998 0.04794 0.34092 1
73 45 0.00121 0.01109 1.4563 1
29 75 0.00051 0.00517 0.1352 1
37 69 0.01212 0.061 0.09912 1
70 72 0.00878 0.0443 0.0716 1
21 65 0.00083 0.00396 0.02862 1
21 30 0.00695 0.035 0.05686 1];
switch num
case 14
linedt = linedat14;
case 30
linedt = linedat30;
case 57
linedt = linedat57;
case 75
linedt=linedat75;
end
