Loading [MathJax]/jax/element/mml/optable/GeneralPunctuation.js

GS_program_FDM_(v1.1)

1265 days ago by fresl


Metoda gustoća sila

(Force Density Method)


[Ćelije u kojima je prvi redak %auto ne trebate izvoditi ( evaluate ); one se izvode „same” pri otvaranju radnoga lista.  Međutim, ako radni list ponovno pokrenete—ne izlazeći iz njega—s  Restart worksheet  (u izborniku  Action... ), i njih ćete morati izvesti.]


       


Primjer 1.:

 

• čvorovi:

       
       
[[0, 0, 0],
 [2, 0, 0],
 [4, 0, 0],
 [6, 0, 0],
 [8, 0, 0],
 [1, 1, 0],
 [3, 1, 0],
 [5, 1, 0],
 [7, 1, 0],
 [0, 2, 0],
 [2, 2, 0],
 [4, 2, 0],
 [6, 2, 0],
 [8, 2, 0],
 [1, 3, 0],
 [3, 3, 0],
 [5, 3, 0],
 [7, 3, 0],
 [0, 4, 0],
 [2, 4, 0],
 [4, 4, 0],
 [6, 4, 0],
 [8, 4, 0],
 [1, 5, 0],
 [3, 5, 0],
 [5, 5, 0],
 [7, 5, 0],
 [0, 6, 0],
 [2, 6, 0],
 [4, 6, 0],
 [6, 6, 0],
 [8, 6, 0],
 [1, 7, 0],
 [3, 7, 0],
 [5, 7, 0],
 [7, 7, 0],
 [0, 8, 0],
 [2, 8, 0],
 [4, 8, 0],
 [6, 8, 0],
 [8, 8, 0]]
[[0, 0, 0],
 [2, 0, 0],
 [4, 0, 0],
 [6, 0, 0],
 [8, 0, 0],
 [1, 1, 0],
 [3, 1, 0],
 [5, 1, 0],
 [7, 1, 0],
 [0, 2, 0],
 [2, 2, 0],
 [4, 2, 0],
 [6, 2, 0],
 [8, 2, 0],
 [1, 3, 0],
 [3, 3, 0],
 [5, 3, 0],
 [7, 3, 0],
 [0, 4, 0],
 [2, 4, 0],
 [4, 4, 0],
 [6, 4, 0],
 [8, 4, 0],
 [1, 5, 0],
 [3, 5, 0],
 [5, 5, 0],
 [7, 5, 0],
 [0, 6, 0],
 [2, 6, 0],
 [4, 6, 0],
 [6, 6, 0],
 [8, 6, 0],
 [1, 7, 0],
 [3, 7, 0],
 [5, 7, 0],
 [7, 7, 0],
 [0, 8, 0],
 [2, 8, 0],
 [4, 8, 0],
 [6, 8, 0],
 [8, 8, 0]]
       
       
0 :  [0, 0, 0]
1 :  [2, 0, 0]
2 :  [4, 0, 0]
3 :  [6, 0, 0]
4 :  [8, 0, 0]
5 :  [1, 1, 0]
6 :  [3, 1, 0]
7 :  [5, 1, 0]
8 :  [7, 1, 0]
9 :  [0, 2, 0]
10 :  [2, 2, 0]
11 :  [4, 2, 0]
12 :  [6, 2, 0]
13 :  [8, 2, 0]
14 :  [1, 3, 0]
15 :  [3, 3, 0]
16 :  [5, 3, 0]
17 :  [7, 3, 0]
18 :  [0, 4, 0]
19 :  [2, 4, 0]
20 :  [4, 4, 0]
21 :  [6, 4, 0]
22 :  [8, 4, 0]
23 :  [1, 5, 0]
24 :  [3, 5, 0]
25 :  [5, 5, 0]
26 :  [7, 5, 0]
27 :  [0, 6, 0]
28 :  [2, 6, 0]
29 :  [4, 6, 0]
30 :  [6, 6, 0]
31 :  [8, 6, 0]
32 :  [1, 7, 0]
33 :  [3, 7, 0]
34 :  [5, 7, 0]
35 :  [7, 7, 0]
36 :  [0, 8, 0]
37 :  [2, 8, 0]
38 :  [4, 8, 0]
39 :  [6, 8, 0]
40 :  [8, 8, 0]
0 :  [0, 0, 0]
1 :  [2, 0, 0]
2 :  [4, 0, 0]
3 :  [6, 0, 0]
4 :  [8, 0, 0]
5 :  [1, 1, 0]
6 :  [3, 1, 0]
7 :  [5, 1, 0]
8 :  [7, 1, 0]
9 :  [0, 2, 0]
10 :  [2, 2, 0]
11 :  [4, 2, 0]
12 :  [6, 2, 0]
13 :  [8, 2, 0]
14 :  [1, 3, 0]
15 :  [3, 3, 0]
16 :  [5, 3, 0]
17 :  [7, 3, 0]
18 :  [0, 4, 0]
19 :  [2, 4, 0]
20 :  [4, 4, 0]
21 :  [6, 4, 0]
22 :  [8, 4, 0]
23 :  [1, 5, 0]
24 :  [3, 5, 0]
25 :  [5, 5, 0]
26 :  [7, 5, 0]
27 :  [0, 6, 0]
28 :  [2, 6, 0]
29 :  [4, 6, 0]
30 :  [6, 6, 0]
31 :  [8, 6, 0]
32 :  [1, 7, 0]
33 :  [3, 7, 0]
34 :  [5, 7, 0]
35 :  [7, 7, 0]
36 :  [0, 8, 0]
37 :  [2, 8, 0]
38 :  [4, 8, 0]
39 :  [6, 8, 0]
40 :  [8, 8, 0]

• kabeli i štapovi:

       
0 :  [27, 32, 37]
1 :  [18, 23, 28, 33, 38]
2 :  [9, 14, 19, 24, 29, 34, 39]
3 :  [0, 5, 10, 15, 20, 25, 30, 35, 40]
4 :  [1, 6, 11, 16, 21, 26, 31]
5 :  [2, 7, 12, 17, 22]
6 :  [3, 8, 13]
7 :  [1, 5, 9]
8 :  [2, 6, 10, 14, 18]
9 :  [3, 7, 11, 15, 19, 23, 27]
10 :  [4, 8, 12, 16, 20, 24, 28, 32, 36]
11 :  [13, 17, 21, 25, 29, 33, 37]
12 :  [22, 26, 30, 34, 38]
13 :  [31, 35, 39]
0 :  [27, 32, 37]
1 :  [18, 23, 28, 33, 38]
2 :  [9, 14, 19, 24, 29, 34, 39]
3 :  [0, 5, 10, 15, 20, 25, 30, 35, 40]
4 :  [1, 6, 11, 16, 21, 26, 31]
5 :  [2, 7, 12, 17, 22]
6 :  [3, 8, 13]
7 :  [1, 5, 9]
8 :  [2, 6, 10, 14, 18]
9 :  [3, 7, 11, 15, 19, 23, 27]
10 :  [4, 8, 12, 16, 20, 24, 28, 32, 36]
11 :  [13, 17, 21, 25, 29, 33, 37]
12 :  [22, 26, 30, 34, 38]
13 :  [31, 35, 39]
       
0 :  [27, 32]
1 :  [32, 37]
2 :  [18, 23]
3 :  [23, 28]
4 :  [28, 33]
5 :  [33, 38]
6 :  [9, 14]
7 :  [14, 19]
8 :  [19, 24]
9 :  [24, 29]
10 :  [29, 34]
11 :  [34, 39]
12 :  [0, 5]
13 :  [5, 10]
14 :  [10, 15]
15 :  [15, 20]
16 :  [20, 25]
17 :  [25, 30]
18 :  [30, 35]
19 :  [35, 40]
20 :  [1, 6]
21 :  [6, 11]
22 :  [11, 16]
23 :  [16, 21]
24 :  [21, 26]
25 :  [26, 31]
26 :  [2, 7]
27 :  [7, 12]
28 :  [12, 17]
29 :  [17, 22]
30 :  [3, 8]
31 :  [8, 13]
32 :  [1, 5]
33 :  [5, 9]
34 :  [2, 6]
35 :  [6, 10]
36 :  [10, 14]
37 :  [14, 18]
38 :  [3, 7]
39 :  [7, 11]
40 :  [11, 15]
41 :  [15, 19]
42 :  [19, 23]
43 :  [23, 27]
44 :  [4, 8]
45 :  [8, 12]
46 :  [12, 16]
47 :  [16, 20]
48 :  [20, 24]
49 :  [24, 28]
50 :  [28, 32]
51 :  [32, 36]
52 :  [13, 17]
53 :  [17, 21]
54 :  [21, 25]
55 :  [25, 29]
56 :  [29, 33]
57 :  [33, 37]
58 :  [22, 26]
59 :  [26, 30]
60 :  [30, 34]
61 :  [34, 38]
62 :  [31, 35]
63 :  [35, 39]
0 :  [27, 32]
1 :  [32, 37]
2 :  [18, 23]
3 :  [23, 28]
4 :  [28, 33]
5 :  [33, 38]
6 :  [9, 14]
7 :  [14, 19]
8 :  [19, 24]
9 :  [24, 29]
10 :  [29, 34]
11 :  [34, 39]
12 :  [0, 5]
13 :  [5, 10]
14 :  [10, 15]
15 :  [15, 20]
16 :  [20, 25]
17 :  [25, 30]
18 :  [30, 35]
19 :  [35, 40]
20 :  [1, 6]
21 :  [6, 11]
22 :  [11, 16]
23 :  [16, 21]
24 :  [21, 26]
25 :  [26, 31]
26 :  [2, 7]
27 :  [7, 12]
28 :  [12, 17]
29 :  [17, 22]
30 :  [3, 8]
31 :  [8, 13]
32 :  [1, 5]
33 :  [5, 9]
34 :  [2, 6]
35 :  [6, 10]
36 :  [10, 14]
37 :  [14, 18]
38 :  [3, 7]
39 :  [7, 11]
40 :  [11, 15]
41 :  [15, 19]
42 :  [19, 23]
43 :  [23, 27]
44 :  [4, 8]
45 :  [8, 12]
46 :  [12, 16]
47 :  [16, 20]
48 :  [20, 24]
49 :  [24, 28]
50 :  [28, 32]
51 :  [32, 36]
52 :  [13, 17]
53 :  [17, 21]
54 :  [21, 25]
55 :  [25, 29]
56 :  [29, 33]
57 :  [33, 37]
58 :  [22, 26]
59 :  [26, 30]
60 :  [30, 34]
61 :  [34, 38]
62 :  [31, 35]
63 :  [35, 39]
       
       

• ležajevi:

       
       
       
       


• drugi način zadavanja čvorova i ležajeva:

       
0 :  []
1 :  []
2 :  []
3 :  []
4 :  []
5 :  []
6 :  []
7 :  []
8 :  []
9 :  []
10 :  []
11 :  []
12 :  []
13 :  []
14 :  []
15 :  []
16 :  []
17 :  []
18 :  []
19 :  []
20 :  []
21 :  []
22 :  []
23 :  []
24 :  []
25 :  []
26 :  []
27 :  []
28 :  []
29 :  []
30 :  []
31 :  []
32 :  []
33 :  []
34 :  []
35 :  []
36 :  []
37 :  []
38 :  []
39 :  []
40 :  []
0 :  []
1 :  []
2 :  []
3 :  []
4 :  []
5 :  []
6 :  []
7 :  []
8 :  []
9 :  []
10 :  []
11 :  []
12 :  []
13 :  []
14 :  []
15 :  []
16 :  []
17 :  []
18 :  []
19 :  []
20 :  []
21 :  []
22 :  []
23 :  []
24 :  []
25 :  []
26 :  []
27 :  []
28 :  []
29 :  []
30 :  []
31 :  []
32 :  []
33 :  []
34 :  []
35 :  []
36 :  []
37 :  []
38 :  []
39 :  []
40 :  []
       
0 :  [0, 0, 0]
1 :  [2, 0, 0]
2 :  [4, 0, 0]
3 :  [6, 0, 0]
4 :  [8, 0, 0]
5 :  []
6 :  []
7 :  []
8 :  []
9 :  [0, 2, 0]
10 :  []
11 :  []
12 :  []
13 :  [8, 2, 1]
14 :  []
15 :  []
16 :  []
17 :  []
18 :  [0, 4, 0]
19 :  []
20 :  []
21 :  []
22 :  [8, 4, 2]
23 :  []
24 :  []
25 :  []
26 :  []
27 :  [0, 6, 0]
28 :  []
29 :  []
30 :  []
31 :  [8, 6, 3]
32 :  []
33 :  []
34 :  []
35 :  []
36 :  [0, 8, 0]
37 :  [2, 8, 1]
38 :  [4, 8, 2]
39 :  [6, 8, 3]
40 :  [8, 8, 4]
0 :  [0, 0, 0]
1 :  [2, 0, 0]
2 :  [4, 0, 0]
3 :  [6, 0, 0]
4 :  [8, 0, 0]
5 :  []
6 :  []
7 :  []
8 :  []
9 :  [0, 2, 0]
10 :  []
11 :  []
12 :  []
13 :  [8, 2, 1]
14 :  []
15 :  []
16 :  []
17 :  []
18 :  [0, 4, 0]
19 :  []
20 :  []
21 :  []
22 :  [8, 4, 2]
23 :  []
24 :  []
25 :  []
26 :  []
27 :  [0, 6, 0]
28 :  []
29 :  []
30 :  []
31 :  [8, 6, 3]
32 :  []
33 :  []
34 :  []
35 :  []
36 :  [0, 8, 0]
37 :  [2, 8, 1]
38 :  [4, 8, 2]
39 :  [6, 8, 3]
40 :  [8, 8, 4]
       
[0, 1, 2, 3, 4, 9, 13, 18, 22, 27, 31, 36, 37, 38, 39, 40]
[0, 1, 2, 3, 4, 9, 13, 18, 22, 27, 31, 36, 37, 38, 39, 40]
       
Traceback (click to the left of this block for traceback)
...
IndexError: list index out of range
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "_sage_input_20.py", line 10, in <module>
    exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8 -*-\\n" + _support_.preparse_worksheet_cell(base64.b64decode("cGxvdDNkX21lc2ggKG5kcywgZWxzKQ=="),globals())+"\\n"); execfile(os.path.abspath("___code___.py"))
  File "", line 1, in <module>
    
  File "/tmp/tmpvRJ4hq/___code___.py", line 2, in <module>
    exec compile(u'plot3d_mesh (nds, els)
  File "", line 1, in <module>
    
  File "<string>", line 599, in plot3d_mesh
  File "<string>", line 574, in plot3d_nodes
  File "/opt/SageMath/local/lib/python2.7/site-packages/sage/misc/decorators.py", line 653, in wrapper
    return func(*args, **kwds)
  File "/opt/SageMath/local/lib/python2.7/site-packages/sage/plot/plot3d/shapes2.py", line 1175, in point3d
    A = sum([Point(z, size, **kwds) for z in v])
  File "/opt/SageMath/local/lib/python2.7/site-packages/sage/plot/plot3d/shapes2.py", line 742, in __init__
    self.loc = (float(center[0]), float(center[1]), float(center[2]))
IndexError: list index out of range


pripadnost elemenata kabelima:

       
0 :  [0, 1]
1 :  [2, 3, 4, 5]
2 :  [6, 7, 8, 9, 10, 11]
3 :  [12, 13, 14, 15, 16, 17, 18, 19]
4 :  [20, 21, 22, 23, 24, 25]
5 :  [26, 27, 28, 29]
6 :  [30, 31]
7 :  [32, 33]
8 :  [34, 35, 36, 37]
9 :  [38, 39, 40, 41, 42, 43]
10 :  [44, 45, 46, 47, 48, 49, 50, 51]
11 :  [52, 53, 54, 55, 56, 57]
12 :  [58, 59, 60, 61]
13 :  [62, 63]
0 :  [0, 1]
1 :  [2, 3, 4, 5]
2 :  [6, 7, 8, 9, 10, 11]
3 :  [12, 13, 14, 15, 16, 17, 18, 19]
4 :  [20, 21, 22, 23, 24, 25]
5 :  [26, 27, 28, 29]
6 :  [30, 31]
7 :  [32, 33]
8 :  [34, 35, 36, 37]
9 :  [38, 39, 40, 41, 42, 43]
10 :  [44, 45, 46, 47, 48, 49, 50, 51]
11 :  [52, 53, 54, 55, 56, 57]
12 :  [58, 59, 60, 61]
13 :  [62, 63]

 

• gustoće sila:

       
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
       
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  1
15 :  1
16 :  1
17 :  1
18 :  1
19 :  1
20 :  1
21 :  1
22 :  1
23 :  1
24 :  1
25 :  1
26 :  1
27 :  1
28 :  1
29 :  1
30 :  1
31 :  1
32 :  1
33 :  1
34 :  1
35 :  1
36 :  1
37 :  1
38 :  1
39 :  1
40 :  1
41 :  1
42 :  1
43 :  1
44 :  1
45 :  1
46 :  1
47 :  1
48 :  1
49 :  1
50 :  1
51 :  1
52 :  1
53 :  1
54 :  1
55 :  1
56 :  1
57 :  1
58 :  1
59 :  1
60 :  1
61 :  1
62 :  1
63 :  1
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  1
15 :  1
16 :  1
17 :  1
18 :  1
19 :  1
20 :  1
21 :  1
22 :  1
23 :  1
24 :  1
25 :  1
26 :  1
27 :  1
28 :  1
29 :  1
30 :  1
31 :  1
32 :  1
33 :  1
34 :  1
35 :  1
36 :  1
37 :  1
38 :  1
39 :  1
40 :  1
41 :  1
42 :  1
43 :  1
44 :  1
45 :  1
46 :  1
47 :  1
48 :  1
49 :  1
50 :  1
51 :  1
52 :  1
53 :  1
54 :  1
55 :  1
56 :  1
57 :  1
58 :  1
59 :  1
60 :  1
61 :  1
62 :  1
63 :  1

 

• proračun metodom gustoća sila:

ncoords = FDM (nodes, elems, supports, qs)

parametri:

nodes — lista koordinata čvorova; zadane moraju biti koordinate ležajnih čvorova

elems — lista štapova određenih indeksima početnoga i krajnjeg čvora

supports — lista indeksa ležajnih čvorova

qs — lista gustoća sila u štapovima

rezultat:

ncoords — matrica tipa nˇc×3 koja sadrži koordinate čvorova upisane po redcima

       
       
0 :  (0.0, 0.0, 0.0)
1 :  (2.0, 0.0, 0.0)
2 :  (4.0, 0.0, 0.0)
3 :  (6.0, 0.0, 0.0)
4 :  (8.0, 0.0, 0.0)
5 :  (1.0, 0.9999999999999999, 0.062499999999999986)
6 :  (3.0, 0.9999999999999998, 0.18749999999999997)
7 :  (5.0, 0.9999999999999998, 0.3125)
8 :  (7.0, 1.0, 0.4375)
9 :  (0.0, 2.0, 0.0)
10 :  (2.0, 1.9999999999999996, 0.24999999999999994)
11 :  (3.9999999999999996, 1.9999999999999996, 0.49999999999999994)
12 :  (5.999999999999999, 1.9999999999999998, 0.7500000000000001)
13 :  (8.0, 2.0, 1.0)
14 :  (0.9999999999999999, 2.9999999999999996, 0.18749999999999997)
15 :  (2.9999999999999996, 2.999999999999999, 0.5624999999999998)
16 :  (5.0, 2.999999999999999, 0.9375)
17 :  (6.999999999999999, 3.0000000000000004, 1.3125)
18 :  (0.0, 4.0, 0.0)
19 :  (1.9999999999999998, 3.9999999999999996, 0.49999999999999994)
20 :  (3.9999999999999996, 4.0, 0.9999999999999998)
21 :  (6.0, 4.0, 1.4999999999999998)
22 :  (8.0, 4.0, 2.0)
23 :  (1.0, 4.999999999999999, 0.31250000000000006)
24 :  (3.0, 5.000000000000001, 0.9374999999999999)
25 :  (5.000000000000001, 5.0, 1.5624999999999998)
26 :  (7.000000000000001, 5.0, 2.1875)
27 :  (0.0, 6.0, 0.0)
28 :  (2.0000000000000004, 6.0, 0.7500000000000001)
29 :  (4.0, 6.0, 1.4999999999999998)
30 :  (6.000000000000002, 5.999999999999999, 2.25)
31 :  (8.0, 6.0, 3.0)
32 :  (1.0000000000000002, 6.999999999999999, 0.43750000000000006)
33 :  (3.0000000000000004, 7.0, 1.3125000000000002)
34 :  (5.0, 7.0, 2.1875)
35 :  (7.000000000000001, 7.000000000000001, 3.0625)
36 :  (0.0, 8.0, 0.0)
37 :  (2.0, 8.0, 1.0)
38 :  (4.0, 8.0, 2.0)
39 :  (6.0, 8.0, 3.0)
40 :  (8.0, 8.0, 4.0)
0 :  (0.0, 0.0, 0.0)
1 :  (2.0, 0.0, 0.0)
2 :  (4.0, 0.0, 0.0)
3 :  (6.0, 0.0, 0.0)
4 :  (8.0, 0.0, 0.0)
5 :  (1.0, 0.9999999999999999, 0.062499999999999986)
6 :  (3.0, 0.9999999999999998, 0.18749999999999997)
7 :  (5.0, 0.9999999999999998, 0.3125)
8 :  (7.0, 1.0, 0.4375)
9 :  (0.0, 2.0, 0.0)
10 :  (2.0, 1.9999999999999996, 0.24999999999999994)
11 :  (3.9999999999999996, 1.9999999999999996, 0.49999999999999994)
12 :  (5.999999999999999, 1.9999999999999998, 0.7500000000000001)
13 :  (8.0, 2.0, 1.0)
14 :  (0.9999999999999999, 2.9999999999999996, 0.18749999999999997)
15 :  (2.9999999999999996, 2.999999999999999, 0.5624999999999998)
16 :  (5.0, 2.999999999999999, 0.9375)
17 :  (6.999999999999999, 3.0000000000000004, 1.3125)
18 :  (0.0, 4.0, 0.0)
19 :  (1.9999999999999998, 3.9999999999999996, 0.49999999999999994)
20 :  (3.9999999999999996, 4.0, 0.9999999999999998)
21 :  (6.0, 4.0, 1.4999999999999998)
22 :  (8.0, 4.0, 2.0)
23 :  (1.0, 4.999999999999999, 0.31250000000000006)
24 :  (3.0, 5.000000000000001, 0.9374999999999999)
25 :  (5.000000000000001, 5.0, 1.5624999999999998)
26 :  (7.000000000000001, 5.0, 2.1875)
27 :  (0.0, 6.0, 0.0)
28 :  (2.0000000000000004, 6.0, 0.7500000000000001)
29 :  (4.0, 6.0, 1.4999999999999998)
30 :  (6.000000000000002, 5.999999999999999, 2.25)
31 :  (8.0, 6.0, 3.0)
32 :  (1.0000000000000002, 6.999999999999999, 0.43750000000000006)
33 :  (3.0000000000000004, 7.0, 1.3125000000000002)
34 :  (5.0, 7.0, 2.1875)
35 :  (7.000000000000001, 7.000000000000001, 3.0625)
36 :  (0.0, 8.0, 0.0)
37 :  (2.0, 8.0, 1.0)
38 :  (4.0, 8.0, 2.0)
39 :  (6.0, 8.0, 3.0)
40 :  (8.0, 8.0, 4.0)
       
       
0 :  1.4803399102908759
1 :  1.5219744577357406
2 :  1.4483287782820577
3 :  1.4803399102908774
4 :  1.5219744577357401
5 :  1.5724682031761401
6 :  1.4265890263141654
7 :  1.4483287782820584
8 :  1.4803399102908774
9 :  1.5219744577357395
10 :  1.5724682031761408
11 :  1.630998543837486
12 :  1.4155939566132654
13 :  1.4265890263141656
14 :  1.4483287782820577
15 :  1.480339910290877
16 :  1.521974457735741
17 :  1.5724682031761408
18 :  1.6309985438374865
19 :  1.6967339950622775
20 :  1.4265890263141658
21 :  1.448328778282058
22 :  1.4803399102908763
23 :  1.5219744577357406
24 :  1.5724682031761412
25 :  1.6309985438374854
26 :  1.4483287782820582
27 :  1.4803399102908759
28 :  1.5219744577357406
29 :  1.5724682031761408
30 :  1.4803399102908763
31 :  1.5219744577357401
32 :  1.4155939566132654
33 :  1.4155939566132656
34 :  1.4265890263141658
35 :  1.4155939566132654
36 :  1.4155939566132656
37 :  1.4265890263141663
38 :  1.4483287782820582
39 :  1.426589026314166
40 :  1.4155939566132651
41 :  1.4155939566132658
42 :  1.4265890263141654
43 :  1.448328778282059
44 :  1.4803399102908763
45 :  1.4483287782820589
46 :  1.4265890263141647
47 :  1.4155939566132665
48 :  1.4155939566132658
49 :  1.426589026314165
50 :  1.448328778282058
51 :  1.4803399102908772
52 :  1.4483287782820593
53 :  1.426589026314165
54 :  1.415593956613265
55 :  1.4155939566132663
56 :  1.4265890263141656
57 :  1.4483287782820586
58 :  1.4265890263141654
59 :  1.4155939566132643
60 :  1.4155939566132674
61 :  1.4265890263141658
62 :  1.4155939566132656
63 :  1.4155939566132656
0 :  1.4803399102908759
1 :  1.5219744577357406
2 :  1.4483287782820577
3 :  1.4803399102908774
4 :  1.5219744577357401
5 :  1.5724682031761401
6 :  1.4265890263141654
7 :  1.4483287782820584
8 :  1.4803399102908774
9 :  1.5219744577357395
10 :  1.5724682031761408
11 :  1.630998543837486
12 :  1.4155939566132654
13 :  1.4265890263141656
14 :  1.4483287782820577
15 :  1.480339910290877
16 :  1.521974457735741
17 :  1.5724682031761408
18 :  1.6309985438374865
19 :  1.6967339950622775
20 :  1.4265890263141658
21 :  1.448328778282058
22 :  1.4803399102908763
23 :  1.5219744577357406
24 :  1.5724682031761412
25 :  1.6309985438374854
26 :  1.4483287782820582
27 :  1.4803399102908759
28 :  1.5219744577357406
29 :  1.5724682031761408
30 :  1.4803399102908763
31 :  1.5219744577357401
32 :  1.4155939566132654
33 :  1.4155939566132656
34 :  1.4265890263141658
35 :  1.4155939566132654
36 :  1.4155939566132656
37 :  1.4265890263141663
38 :  1.4483287782820582
39 :  1.426589026314166
40 :  1.4155939566132651
41 :  1.4155939566132658
42 :  1.4265890263141654
43 :  1.448328778282059
44 :  1.4803399102908763
45 :  1.4483287782820589
46 :  1.4265890263141647
47 :  1.4155939566132665
48 :  1.4155939566132658
49 :  1.426589026314165
50 :  1.448328778282058
51 :  1.4803399102908772
52 :  1.4483287782820593
53 :  1.426589026314165
54 :  1.415593956613265
55 :  1.4155939566132663
56 :  1.4265890263141656
57 :  1.4483287782820586
58 :  1.4265890263141654
59 :  1.4155939566132643
60 :  1.4155939566132674
61 :  1.4265890263141658
62 :  1.4155939566132656
63 :  1.4155939566132656
       
0 :  1.4803399102908759
1 :  1.5219744577357406
2 :  1.4483287782820577
3 :  1.4803399102908774
4 :  1.5219744577357401
5 :  1.5724682031761401
6 :  1.4265890263141654
7 :  1.4483287782820584
8 :  1.4803399102908774
9 :  1.5219744577357395
10 :  1.5724682031761408
11 :  1.630998543837486
12 :  1.4155939566132654
13 :  1.4265890263141656
14 :  1.4483287782820577
15 :  1.480339910290877
16 :  1.521974457735741
17 :  1.5724682031761408
18 :  1.6309985438374865
19 :  1.6967339950622775
20 :  1.4265890263141658
21 :  1.448328778282058
22 :  1.4803399102908763
23 :  1.5219744577357406
24 :  1.5724682031761412
25 :  1.6309985438374854
26 :  1.4483287782820582
27 :  1.4803399102908759
28 :  1.5219744577357406
29 :  1.5724682031761408
30 :  1.4803399102908763
31 :  1.5219744577357401
32 :  1.4155939566132654
33 :  1.4155939566132656
34 :  1.4265890263141658
35 :  1.4155939566132654
36 :  1.4155939566132656
37 :  1.4265890263141663
38 :  1.4483287782820582
39 :  1.426589026314166
40 :  1.4155939566132651
41 :  1.4155939566132658
42 :  1.4265890263141654
43 :  1.448328778282059
44 :  1.4803399102908763
45 :  1.4483287782820589
46 :  1.4265890263141647
47 :  1.4155939566132665
48 :  1.4155939566132658
49 :  1.426589026314165
50 :  1.448328778282058
51 :  1.4803399102908772
52 :  1.4483287782820593
53 :  1.426589026314165
54 :  1.415593956613265
55 :  1.4155939566132663
56 :  1.4265890263141656
57 :  1.4483287782820586
58 :  1.4265890263141654
59 :  1.4155939566132643
60 :  1.4155939566132674
61 :  1.4265890263141658
62 :  1.4155939566132656
63 :  1.4155939566132656
0 :  1.4803399102908759
1 :  1.5219744577357406
2 :  1.4483287782820577
3 :  1.4803399102908774
4 :  1.5219744577357401
5 :  1.5724682031761401
6 :  1.4265890263141654
7 :  1.4483287782820584
8 :  1.4803399102908774
9 :  1.5219744577357395
10 :  1.5724682031761408
11 :  1.630998543837486
12 :  1.4155939566132654
13 :  1.4265890263141656
14 :  1.4483287782820577
15 :  1.480339910290877
16 :  1.521974457735741
17 :  1.5724682031761408
18 :  1.6309985438374865
19 :  1.6967339950622775
20 :  1.4265890263141658
21 :  1.448328778282058
22 :  1.4803399102908763
23 :  1.5219744577357406
24 :  1.5724682031761412
25 :  1.6309985438374854
26 :  1.4483287782820582
27 :  1.4803399102908759
28 :  1.5219744577357406
29 :  1.5724682031761408
30 :  1.4803399102908763
31 :  1.5219744577357401
32 :  1.4155939566132654
33 :  1.4155939566132656
34 :  1.4265890263141658
35 :  1.4155939566132654
36 :  1.4155939566132656
37 :  1.4265890263141663
38 :  1.4483287782820582
39 :  1.426589026314166
40 :  1.4155939566132651
41 :  1.4155939566132658
42 :  1.4265890263141654
43 :  1.448328778282059
44 :  1.4803399102908763
45 :  1.4483287782820589
46 :  1.4265890263141647
47 :  1.4155939566132665
48 :  1.4155939566132658
49 :  1.426589026314165
50 :  1.448328778282058
51 :  1.4803399102908772
52 :  1.4483287782820593
53 :  1.426589026314165
54 :  1.415593956613265
55 :  1.4155939566132663
56 :  1.4265890263141656
57 :  1.4483287782820586
58 :  1.4265890263141654
59 :  1.4155939566132643
60 :  1.4155939566132674
61 :  1.4265890263141658
62 :  1.4155939566132656
63 :  1.4155939566132656

 

• zadavanje vrijednostî sila u kabelima:

       
0 :  1.4
1 :  1.4
2 :  1.4
3 :  1.4
4 :  1.4
5 :  1.4
6 :  1.4
7 :  1.4
8 :  1.4
9 :  1.4
10 :  1.4
11 :  1.4
12 :  1.4
13 :  1.4
0 :  1.4
1 :  1.4
2 :  1.4
3 :  1.4
4 :  1.4
5 :  1.4
6 :  1.4
7 :  1.4
8 :  1.4
9 :  1.4
10 :  1.4
11 :  1.4
12 :  1.4
13 :  1.4
       
0 :  1.4
1 :  1.4
2 :  1.4
3 :  1.4
4 :  1.4
5 :  1.4
6 :  1.4
7 :  1.4
8 :  1.4
9 :  1.4
10 :  1.4
11 :  1.4
12 :  1.4
13 :  1.4
14 :  1.4
15 :  1.4
16 :  1.4
17 :  1.4
18 :  1.4
19 :  1.4
20 :  1.4
21 :  1.4
22 :  1.4
23 :  1.4
24 :  1.4
25 :  1.4
26 :  1.4
27 :  1.4
28 :  1.4
29 :  1.4
30 :  1.4
31 :  1.4
32 :  1.4
33 :  1.4
34 :  1.4
35 :  1.4
36 :  1.4
37 :  1.4
38 :  1.4
39 :  1.4
40 :  1.4
41 :  1.4
42 :  1.4
43 :  1.4
44 :  1.4
45 :  1.4
46 :  1.4
47 :  1.4
48 :  1.4
49 :  1.4
50 :  1.4
51 :  1.4
52 :  1.4
53 :  1.4
54 :  1.4
55 :  1.4
56 :  1.4
57 :  1.4
58 :  1.4
59 :  1.4
60 :  1.4
61 :  1.4
62 :  1.4
63 :  1.4
0 :  1.4
1 :  1.4
2 :  1.4
3 :  1.4
4 :  1.4
5 :  1.4
6 :  1.4
7 :  1.4
8 :  1.4
9 :  1.4
10 :  1.4
11 :  1.4
12 :  1.4
13 :  1.4
14 :  1.4
15 :  1.4
16 :  1.4
17 :  1.4
18 :  1.4
19 :  1.4
20 :  1.4
21 :  1.4
22 :  1.4
23 :  1.4
24 :  1.4
25 :  1.4
26 :  1.4
27 :  1.4
28 :  1.4
29 :  1.4
30 :  1.4
31 :  1.4
32 :  1.4
33 :  1.4
34 :  1.4
35 :  1.4
36 :  1.4
37 :  1.4
38 :  1.4
39 :  1.4
40 :  1.4
41 :  1.4
42 :  1.4
43 :  1.4
44 :  1.4
45 :  1.4
46 :  1.4
47 :  1.4
48 :  1.4
49 :  1.4
50 :  1.4
51 :  1.4
52 :  1.4
53 :  1.4
54 :  1.4
55 :  1.4
56 :  1.4
57 :  1.4
58 :  1.4
59 :  1.4
60 :  1.4
61 :  1.4
62 :  1.4
63 :  1.4
       
[(0, 1.4),
 (1, 1.4),
 (2, 1.4),
 (3, 1.4),
 (4, 1.4),
 (5, 1.4),
 (6, 1.4),
 (7, 1.4),
 (8, 1.4),
 (9, 1.4),
 (10, 1.4),
 (11, 1.4),
 (12, 1.4),
 (13, 1.4),
 (14, 1.4),
 (15, 1.4),
 (16, 1.4),
 (17, 1.4),
 (18, 1.4),
 (19, 1.4),
 (20, 1.4),
 (21, 1.4),
 (22, 1.4),
 (23, 1.4),
 (24, 1.4),
 (25, 1.4),
 (26, 1.4),
 (27, 1.4),
 (28, 1.4),
 (29, 1.4),
 (30, 1.4),
 (31, 1.4),
 (32, 1.4),
 (33, 1.4),
 (34, 1.4),
 (35, 1.4),
 (36, 1.4),
 (37, 1.4),
 (38, 1.4),
 (39, 1.4),
 (40, 1.4),
 (41, 1.4),
 (42, 1.4),
 (43, 1.4),
 (44, 1.4),
 (45, 1.4),
 (46, 1.4),
 (47, 1.4),
 (48, 1.4),
 (49, 1.4),
 (50, 1.4),
 (51, 1.4),
 (52, 1.4),
 (53, 1.4),
 (54, 1.4),
 (55, 1.4),
 (56, 1.4),
 (57, 1.4),
 (58, 1.4),
 (59, 1.4),
 (60, 1.4),
 (61, 1.4),
 (62, 1.4),
 (63, 1.4)]
[(0, 1.4),
 (1, 1.4),
 (2, 1.4),
 (3, 1.4),
 (4, 1.4),
 (5, 1.4),
 (6, 1.4),
 (7, 1.4),
 (8, 1.4),
 (9, 1.4),
 (10, 1.4),
 (11, 1.4),
 (12, 1.4),
 (13, 1.4),
 (14, 1.4),
 (15, 1.4),
 (16, 1.4),
 (17, 1.4),
 (18, 1.4),
 (19, 1.4),
 (20, 1.4),
 (21, 1.4),
 (22, 1.4),
 (23, 1.4),
 (24, 1.4),
 (25, 1.4),
 (26, 1.4),
 (27, 1.4),
 (28, 1.4),
 (29, 1.4),
 (30, 1.4),
 (31, 1.4),
 (32, 1.4),
 (33, 1.4),
 (34, 1.4),
 (35, 1.4),
 (36, 1.4),
 (37, 1.4),
 (38, 1.4),
 (39, 1.4),
 (40, 1.4),
 (41, 1.4),
 (42, 1.4),
 (43, 1.4),
 (44, 1.4),
 (45, 1.4),
 (46, 1.4),
 (47, 1.4),
 (48, 1.4),
 (49, 1.4),
 (50, 1.4),
 (51, 1.4),
 (52, 1.4),
 (53, 1.4),
 (54, 1.4),
 (55, 1.4),
 (56, 1.4),
 (57, 1.4),
 (58, 1.4),
 (59, 1.4),
 (60, 1.4),
 (61, 1.4),
 (62, 1.4),
 (63, 1.4)]

 

• iteracijska primjena metode gustoća sila odabrani broj puta:

ncoords, forces, f_dens = multistepFDM (nodes, elems, supports, qs, fcs = [], lcs = [], l0cs = [], steps = 250)

parametri:

nodes — lista koordinata čvorova; zadane moraju biti koordinate ležajnih čvorova

elems — lista štapova određenih indeksima početnoga i krajnjeg čvora

supports — lista indeksa ležajnih čvorova

qs — lista gustoća sila u štapovima za prvi korak

fcs — lista zadanih vrijednosti sila u štapovima

komponente liste: parovi (idx, f); idx – indeks štapa, f – vrijednost sile

lcs — lista zadanih duljina štapova

komponente liste: parovi (idx, l); idx – indeks štapa, l – duljina

l0cs — lista zadanih nerastegnutih duljina štapova

komponente liste: parovi (idx, (l0, AE))idx – indeks štapa, l0 – duljina, AE – uzdužna krutost 

steps — broj koraka

rezultati:

ncoords — matrica tipa nˇc×3 koja sadrži po redcima svrstane koordinate čvorova

forces — lista konačnih vrijednosti sila u štapovima

f_dens — lista gustoća sila u štapovima izračunanih u posljednjem koraku

       
       
0 :  1.3990568080406078
1 :  1.4007728237337864
2 :  1.4003698477611368
3 :  1.3979199182680202
4 :  1.4015769404285785
5 :  1.39986120140528
6 :  1.4001792288151318
7 :  1.4003804551383674
8 :  1.397270037700222
9 :  1.402042598351432
10 :  1.3999661890159862
11 :  1.399846330250519
12 :  1.3999587154366229
13 :  1.4002938643650253
14 :  1.4003587208645938
15 :  1.397068436131475
16 :  1.4021925989565944
17 :  1.4000171501685441
18 :  1.3997594204720143
19 :  1.4000239580997809
20 :  1.400179228815131
21 :  1.4003804551383665
22 :  1.3972700377002214
23 :  1.4020425983514322
24 :  1.3999661890159867
25 :  1.3998463302505202
26 :  1.4003698477611342
27 :  1.3979199182680206
28 :  1.4015769404285794
29 :  1.3998612014052814
30 :  1.3990568080406054
31 :  1.400772823733788
32 :  1.4000000430746729
33 :  1.4000000430746729
34 :  1.4001313709819598
35 :  1.3998628744226316
36 :  1.399862874422632
37 :  1.4001313709819596
38 :  1.3999522197444154
39 :  1.4002700413299265
40 :  1.399761015465588
41 :  1.399761015465588
42 :  1.400270041329928
43 :  1.3999522197444134
44 :  1.3999873044162154
45 :  1.3999253052138612
46 :  1.4003380824269656
47 :  1.3997254955735527
48 :  1.3997254955735527
49 :  1.4003380824269667
50 :  1.3999253052138623
51 :  1.399987304416213
52 :  1.3999427699560911
53 :  1.4003021873623076
54 :  1.3997541576677954
55 :  1.3997541576677957
56 :  1.4003021873623056
57 :  1.3999427699560916
58 :  1.4001394115848915
59 :  1.3998567908108037
60 :  1.3998567908108042
61 :  1.4001394115848893
62 :  1.4000000888273185
63 :  1.4000000888273136
0 :  1.3990568080406078
1 :  1.4007728237337864
2 :  1.4003698477611368
3 :  1.3979199182680202
4 :  1.4015769404285785
5 :  1.39986120140528
6 :  1.4001792288151318
7 :  1.4003804551383674
8 :  1.397270037700222
9 :  1.402042598351432
10 :  1.3999661890159862
11 :  1.399846330250519
12 :  1.3999587154366229
13 :  1.4002938643650253
14 :  1.4003587208645938
15 :  1.397068436131475
16 :  1.4021925989565944
17 :  1.4000171501685441
18 :  1.3997594204720143
19 :  1.4000239580997809
20 :  1.400179228815131
21 :  1.4003804551383665
22 :  1.3972700377002214
23 :  1.4020425983514322
24 :  1.3999661890159867
25 :  1.3998463302505202
26 :  1.4003698477611342
27 :  1.3979199182680206
28 :  1.4015769404285794
29 :  1.3998612014052814
30 :  1.3990568080406054
31 :  1.400772823733788
32 :  1.4000000430746729
33 :  1.4000000430746729
34 :  1.4001313709819598
35 :  1.3998628744226316
36 :  1.399862874422632
37 :  1.4001313709819596
38 :  1.3999522197444154
39 :  1.4002700413299265
40 :  1.399761015465588
41 :  1.399761015465588
42 :  1.400270041329928
43 :  1.3999522197444134
44 :  1.3999873044162154
45 :  1.3999253052138612
46 :  1.4003380824269656
47 :  1.3997254955735527
48 :  1.3997254955735527
49 :  1.4003380824269667
50 :  1.3999253052138623
51 :  1.399987304416213
52 :  1.3999427699560911
53 :  1.4003021873623076
54 :  1.3997541576677954
55 :  1.3997541576677957
56 :  1.4003021873623056
57 :  1.3999427699560916
58 :  1.4001394115848915
59 :  1.3998567908108037
60 :  1.3998567908108042
61 :  1.4001394115848893
62 :  1.4000000888273185
63 :  1.4000000888273136
       
0 :  1.0352037288121443
1 :  0.8484556861574744
2 :  1.0264501792827603
3 :  1.0381770336697995
4 :  0.8698610263674347
5 :  0.8223220883026721
6 :  1.0028352613508433
7 :  1.0415811966779887
8 :  1.0389575433676579
9 :  0.8833001896158216
10 :  0.8271369223140367
11 :  0.8144211836231726
12 :  0.9923946429328483
13 :  1.0066277939354111
14 :  1.046507166960226
15 :  1.0390678548265233
16 :  0.8877754948393836
17 :  0.8286766729966295
18 :  0.8160007402172065
19 :  0.8094385871235791
20 :  1.0028352613508427
21 :  1.0415811966779878
22 :  1.0389575433676579
23 :  0.8833001896158224
24 :  0.8271369223140375
25 :  0.8144211836231723
26 :  1.026450179282759
27 :  1.0381770336698
28 :  0.8698610263674346
29 :  0.8223220883026731
30 :  1.0352037288121436
31 :  0.8484556861574747
32 :  0.9890875727047087
33 :  0.9890875727047087
34 :  0.9668956510821726
35 :  1.005231937365145
36 :  1.0052319373651457
37 :  0.9668956510821717
38 :  0.950088076242526
39 :  0.9778893096321719
40 :  1.0147479573698444
41 :  1.0147479573698452
42 :  0.9778893096321715
43 :  0.9500880762425249
44 :  0.9414071083607118
45 :  0.9533088939307233
46 :  0.9820866321793824
47 :  1.0182546027288517
48 :  1.0182546027288513
49 :  0.9820866321793819
50 :  0.9533088939307227
51 :  0.9414071083607115
52 :  0.9505504105462814
53 :  0.9790767904313928
54 :  1.0154778322245286
55 :  1.0154778322245281
56 :  0.9790767904313921
57 :  0.950550410546281
58 :  0.9679752599610358
59 :  1.0056386832227215
60 :  1.0056386832227209
61 :  0.9679752599610352
62 :  0.9893795548555568
63 :  0.9893795548555534
0 :  1.0352037288121443
1 :  0.8484556861574744
2 :  1.0264501792827603
3 :  1.0381770336697995
4 :  0.8698610263674347
5 :  0.8223220883026721
6 :  1.0028352613508433
7 :  1.0415811966779887
8 :  1.0389575433676579
9 :  0.8833001896158216
10 :  0.8271369223140367
11 :  0.8144211836231726
12 :  0.9923946429328483
13 :  1.0066277939354111
14 :  1.046507166960226
15 :  1.0390678548265233
16 :  0.8877754948393836
17 :  0.8286766729966295
18 :  0.8160007402172065
19 :  0.8094385871235791
20 :  1.0028352613508427
21 :  1.0415811966779878
22 :  1.0389575433676579
23 :  0.8833001896158224
24 :  0.8271369223140375
25 :  0.8144211836231723
26 :  1.026450179282759
27 :  1.0381770336698
28 :  0.8698610263674346
29 :  0.8223220883026731
30 :  1.0352037288121436
31 :  0.8484556861574747
32 :  0.9890875727047087
33 :  0.9890875727047087
34 :  0.9668956510821726
35 :  1.005231937365145
36 :  1.0052319373651457
37 :  0.9668956510821717
38 :  0.950088076242526
39 :  0.9778893096321719
40 :  1.0147479573698444
41 :  1.0147479573698452
42 :  0.9778893096321715
43 :  0.9500880762425249
44 :  0.9414071083607118
45 :  0.9533088939307233
46 :  0.9820866321793824
47 :  1.0182546027288517
48 :  1.0182546027288513
49 :  0.9820866321793819
50 :  0.9533088939307227
51 :  0.9414071083607115
52 :  0.9505504105462814
53 :  0.9790767904313928
54 :  1.0154778322245286
55 :  1.0154778322245281
56 :  0.9790767904313921
57 :  0.950550410546281
58 :  0.9679752599610358
59 :  1.0056386832227215
60 :  1.0056386832227209
61 :  0.9679752599610352
62 :  0.9893795548555568
63 :  0.9893795548555534
       

 

• iteracijska primjena metode gustoća sila do postizanja zadane točnosti:

ncoords, forces = multistepFDM_wtol (nodes, elems, supports, qs, fcs = [], lcs = [], l0cs = [], tol_f = 1e-3, tol_l = 1e-3, abs_error = True, steps = 250)

parametri:

nodes — lista koordinata čvorova; zadane moraju biti koordinate ležajnih čvorova

elems — lista štapova određenih indeksima početnoga i krajnjeg čvora

supports — lista indeksa ležajnih čvorova

qs — lista gustoća sila u štapovima za prvi korak

fcs — lista zadanih vrijednosti sila u štapovima

komponente liste: parovi (idx, f); idx – indeks štapa, f – vrijednost sile

lcs — lista zadanih duljina štapova

komponente liste: parovi (idx, l); idx – indeks štapa, l – duljina

l0cs — lista zadanih nerastegnutih duljina štapova

komponente liste: parovi (idx, (l0, AE))idx – indeks štapa, l0 – duljina, AE – uzdužna krutost 

tol_f — točnost za vrijednosti sila

tol_l — točnost za duljine

abs_error — ako je True, izračunava se apsolutna, a ako je False, relativna pogreška

primjerice, za vrijednosti sila:  eabs=  i  e_{\mathrm{rel}} = \dfrac{\Big\|S_{\{i,j\}}^{(k)} - \bar{S}_{\{i,j\}} \Big\|_{\mathrm{max}}}{\bar{S}_{\{i,j\}} }

steps — broj koraka

rezultati:

ncoords — matrica tipa n_{\check{c}}\!\times 3 koja sadrži po redcima svrstane koordinate čvorova

forces — lista konačnih vrijednosti sila u štapovima

       
steps: 82
maximal force error:  9.98477589584379e-05
steps: 82
maximal force error:  9.98477589584379e-05
       
steps: 62
maximal force error:  9.98477589584379e-05
steps: 62
maximal force error:  9.98477589584379e-05
       
0 :  1.3999756024971828
1 :  1.4000196211174494
2 :  1.4000161539507356
3 :  1.3999365563626371
4 :  1.4000466733858947
5 :  1.3999910136841984
6 :  1.3999969092350149
7 :  1.4000306868753198
8 :  1.3999095161441921
9 :  1.4000651986580923
10 :  1.3999834631052133
11 :  1.4000013379188834
12 :  1.4000001277249048
13 :  1.3999956574265136
14 :  1.400035797742343
15 :  1.3999001522410415
16 :  1.4000715887729343
17 :  1.3999808355882508
18 :  1.4000018616629346
19 :  1.399999951807579
20 :  1.3999969092350157
21 :  1.4000306868753207
22 :  1.3999095161441926
23 :  1.4000651986580925
24 :  1.3999834631052126
25 :  1.4000013379188823
26 :  1.4000161539507376
27 :  1.3999365563626387
28 :  1.4000466733858936
29 :  1.3999910136841967
30 :  1.3999756024971832
31 :  1.4000196211174494
32 :  1.3999999985829343
33 :  1.3999999985829343
34 :  1.4000001405542881
35 :  1.3999998559200084
36 :  1.399999855920009
37 :  1.4000001405542888
38 :  1.400000106770757
39 :  1.400000328217458
40 :  1.399999357245428
41 :  1.399999357245429
42 :  1.4000003282174582
43 :  1.4000001067707577
44 :  1.3999999573385042
45 :  1.399999837582753
46 :  1.4000005621542306
47 :  1.3999996074864378
48 :  1.3999996074864387
49 :  1.4000005621542306
50 :  1.399999837582754
51 :  1.3999999573385047
52 :  1.3999997921400842
53 :  1.4000005451886088
54 :  1.3999998167627585
55 :  1.399999816762758
56 :  1.4000005451886088
57 :  1.399999792140086
58 :  1.4000002334751445
59 :  1.3999997692900439
60 :  1.3999997692900437
61 :  1.400000233475145
62 :  1.4000000015978915
63 :  1.4000000015978904
0 :  1.3999756024971828
1 :  1.4000196211174494
2 :  1.4000161539507356
3 :  1.3999365563626371
4 :  1.4000466733858947
5 :  1.3999910136841984
6 :  1.3999969092350149
7 :  1.4000306868753198
8 :  1.3999095161441921
9 :  1.4000651986580923
10 :  1.3999834631052133
11 :  1.4000013379188834
12 :  1.4000001277249048
13 :  1.3999956574265136
14 :  1.400035797742343
15 :  1.3999001522410415
16 :  1.4000715887729343
17 :  1.3999808355882508
18 :  1.4000018616629346
19 :  1.399999951807579
20 :  1.3999969092350157
21 :  1.4000306868753207
22 :  1.3999095161441926
23 :  1.4000651986580925
24 :  1.3999834631052126
25 :  1.4000013379188823
26 :  1.4000161539507376
27 :  1.3999365563626387
28 :  1.4000466733858936
29 :  1.3999910136841967
30 :  1.3999756024971832
31 :  1.4000196211174494
32 :  1.3999999985829343
33 :  1.3999999985829343
34 :  1.4000001405542881
35 :  1.3999998559200084
36 :  1.399999855920009
37 :  1.4000001405542888
38 :  1.400000106770757
39 :  1.400000328217458
40 :  1.399999357245428
41 :  1.399999357245429
42 :  1.4000003282174582
43 :  1.4000001067707577
44 :  1.3999999573385042
45 :  1.399999837582753
46 :  1.4000005621542306
47 :  1.3999996074864378
48 :  1.3999996074864387
49 :  1.4000005621542306
50 :  1.399999837582754
51 :  1.3999999573385047
52 :  1.3999997921400842
53 :  1.4000005451886088
54 :  1.3999998167627585
55 :  1.399999816762758
56 :  1.4000005451886088
57 :  1.399999792140086
58 :  1.4000002334751445
59 :  1.3999997692900439
60 :  1.3999997692900437
61 :  1.400000233475145
62 :  1.4000000015978915
63 :  1.4000000015978904
       
       
       
       
0 :  1.0459076986085656
1 :  0.8414002417647986
2 :  1.02070465473782
3 :  1.064443273650525
4 :  0.853544543230639
5 :  0.8248210384297869
6 :  1.002821350397471
7 :  1.0319501812838598
8 :  1.075470794259526
9 :  0.8607582658712909
10 :  0.830991886298519
11 :  0.8146459629428328
12 :  0.9925443368936756
13 :  1.0064279300534424
14 :  1.0356025892151586
15 :  1.0790292643988164
16 :  0.8630843135303687
17 :  0.8329568753595153
18 :  0.8164159626409812
19 :  0.8093633328181367
20 :  1.0028213503974712
21 :  1.0319501812838605
22 :  1.0754707942595265
23 :  0.8607582658712918
24 :  0.8309918862985177
25 :  0.8146459629428324
26 :  1.020704654737822
27 :  1.0644432736505252
28 :  0.8535445432306397
29 :  0.8248210384297848
30 :  1.045907698608567
31 :  0.8414002417647977
32 :  0.989087612864574
33 :  0.9890876128645741
34 :  0.9660372457914239
35 :  1.0061626917734057
36 :  1.006162691773407
37 :  0.9660372457914235
38 :  0.9504464039394668
39 :  0.9759818248294584
40 :  1.0165010064030755
41 :  1.016501006403077
42 :  0.9759818248294562
43 :  0.9504464039394684
44 :  0.9414539281145508
45 :  0.9539532745639763
46 :  0.9796404134623456
47 :  1.0202130931463154
48 :  1.0202130931463158
49 :  0.979640413462344
50 :  0.9539532745639763
51 :  0.9414539281145519
52 :  0.9510393481684906
53 :  0.9769244497225359
54 :  1.0171785202582262
55 :  1.0171785202582275
56 :  0.9769244497225342
57 :  0.9510393481684919
58 :  0.9670377047026227
59 :  1.0066416335753037
60 :  1.006641633575305
61 :  0.9670377047026216
62 :  0.9893786677408988
63 :  0.989378667740898
0 :  1.0459076986085656
1 :  0.8414002417647986
2 :  1.02070465473782
3 :  1.064443273650525
4 :  0.853544543230639
5 :  0.8248210384297869
6 :  1.002821350397471
7 :  1.0319501812838598
8 :  1.075470794259526
9 :  0.8607582658712909
10 :  0.830991886298519
11 :  0.8146459629428328
12 :  0.9925443368936756
13 :  1.0064279300534424
14 :  1.0356025892151586
15 :  1.0790292643988164
16 :  0.8630843135303687
17 :  0.8329568753595153
18 :  0.8164159626409812
19 :  0.8093633328181367
20 :  1.0028213503974712
21 :  1.0319501812838605
22 :  1.0754707942595265
23 :  0.8607582658712918
24 :  0.8309918862985177
25 :  0.8146459629428324
26 :  1.020704654737822
27 :  1.0644432736505252
28 :  0.8535445432306397
29 :  0.8248210384297848
30 :  1.045907698608567
31 :  0.8414002417647977
32 :  0.989087612864574
33 :  0.9890876128645741
34 :  0.9660372457914239
35 :  1.0061626917734057
36 :  1.006162691773407
37 :  0.9660372457914235
38 :  0.9504464039394668
39 :  0.9759818248294584
40 :  1.0165010064030755
41 :  1.016501006403077
42 :  0.9759818248294562
43 :  0.9504464039394684
44 :  0.9414539281145508
45 :  0.9539532745639763
46 :  0.9796404134623456
47 :  1.0202130931463154
48 :  1.0202130931463158
49 :  0.979640413462344
50 :  0.9539532745639763
51 :  0.9414539281145519
52 :  0.9510393481684906
53 :  0.9769244497225359
54 :  1.0171785202582262
55 :  1.0171785202582275
56 :  0.9769244497225342
57 :  0.9510393481684919
58 :  0.9670377047026227
59 :  1.0066416335753037
60 :  1.006641633575305
61 :  0.9670377047026216
62 :  0.9893786677408988
63 :  0.989378667740898


• povećanje vrijednosti sile u jednom kabelu:

       
0 :  [27, 32, 37]
1 :  [18, 23, 28, 33, 38]
2 :  [9, 14, 19, 24, 29, 34, 39]
3 :  [0, 5, 10, 15, 20, 25, 30, 35, 40]
4 :  [1, 6, 11, 16, 21, 26, 31]
5 :  [2, 7, 12, 17, 22]
6 :  [3, 8, 13]
7 :  [1, 5, 9]
8 :  [2, 6, 10, 14, 18]
9 :  [3, 7, 11, 15, 19, 23, 27]
10 :  [4, 8, 12, 16, 20, 24, 28, 32, 36]
11 :  [13, 17, 21, 25, 29, 33, 37]
12 :  [22, 26, 30, 34, 38]
13 :  [31, 35, 39]
0 :  [27, 32, 37]
1 :  [18, 23, 28, 33, 38]
2 :  [9, 14, 19, 24, 29, 34, 39]
3 :  [0, 5, 10, 15, 20, 25, 30, 35, 40]
4 :  [1, 6, 11, 16, 21, 26, 31]
5 :  [2, 7, 12, 17, 22]
6 :  [3, 8, 13]
7 :  [1, 5, 9]
8 :  [2, 6, 10, 14, 18]
9 :  [3, 7, 11, 15, 19, 23, 27]
10 :  [4, 8, 12, 16, 20, 24, 28, 32, 36]
11 :  [13, 17, 21, 25, 29, 33, 37]
12 :  [22, 26, 30, 34, 38]
13 :  [31, 35, 39]
       
       
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  1
15 :  1
16 :  1
17 :  1
18 :  1
19 :  1
20 :  1
21 :  1
22 :  1
23 :  1
24 :  1
25 :  1
26 :  1
27 :  1
28 :  1
29 :  1
30 :  1
31 :  1
32 :  1
33 :  1
34 :  1
35 :  1
36 :  1
37 :  1
38 :  1
39 :  1
40 :  1
41 :  1
42 :  1
43 :  1
44 :  10
45 :  10
46 :  10
47 :  10
48 :  10
49 :  10
50 :  10
51 :  10
52 :  1
53 :  1
54 :  1
55 :  1
56 :  1
57 :  1
58 :  1
59 :  1
60 :  1
61 :  1
62 :  1
63 :  1
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  1
15 :  1
16 :  1
17 :  1
18 :  1
19 :  1
20 :  1
21 :  1
22 :  1
23 :  1
24 :  1
25 :  1
26 :  1
27 :  1
28 :  1
29 :  1
30 :  1
31 :  1
32 :  1
33 :  1
34 :  1
35 :  1
36 :  1
37 :  1
38 :  1
39 :  1
40 :  1
41 :  1
42 :  1
43 :  1
44 :  10
45 :  10
46 :  10
47 :  10
48 :  10
49 :  10
50 :  10
51 :  10
52 :  1
53 :  1
54 :  1
55 :  1
56 :  1
57 :  1
58 :  1
59 :  1
60 :  1
61 :  1
62 :  1
63 :  1
       
       
       
0 :  1.4291212486183933
1 :  1.621919491642616
2 :  1.421839054737862
3 :  1.4291212486183926
4 :  1.6219194916426158
5 :  1.6515170754548774
6 :  1.4169634671112834
7 :  1.4218390547378612
8 :  1.4291212486183928
9 :  1.6219194916426156
10 :  1.6515170754548774
11 :  1.682651145927977
12 :  1.4145193713424764
13 :  1.4169634671112827
14 :  1.4218390547378605
15 :  1.4291212486183935
16 :  1.621919491642617
17 :  1.651517075454878
18 :  1.6826511459279758
19 :  1.7152380366850977
20 :  1.4169634671112827
21 :  1.4218390547378605
22 :  1.4291212486183922
23 :  1.621919491642618
24 :  1.6515170754548767
25 :  1.682651145927977
26 :  1.421839054737861
27 :  1.4291212486183926
28 :  1.6219194916426167
29 :  1.6515170754548774
30 :  1.429121248618393
31 :  1.621919491642616
32 :  1.4145193713424764
33 :  1.4145193713424764
34 :  1.4169634671112827
35 :  1.4145193713424764
36 :  1.4145193713424769
37 :  1.4169634671112823
38 :  1.421839054737861
39 :  1.4169634671112827
40 :  1.4145193713424764
41 :  1.4145193713424769
42 :  1.4169634671112827
43 :  1.4218390547378608
44 :  14.29121248618393
45 :  14.218390547378618
46 :  14.169634671112831
47 :  14.145193713424757
48 :  14.145193713424757
49 :  14.169634671112828
50 :  14.218390547378615
51 :  14.291212486183934
52 :  1.4218390547378612
53 :  1.4169634671112834
54 :  1.4145193713424764
55 :  1.4145193713424757
56 :  1.4169634671112827
57 :  1.4218390547378612
58 :  1.4169634671112827
59 :  1.4145193713424764
60 :  1.4145193713424775
61 :  1.4169634671112816
62 :  1.4145193713424764
63 :  1.4145193713424764
0 :  1.4291212486183933
1 :  1.621919491642616
2 :  1.421839054737862
3 :  1.4291212486183926
4 :  1.6219194916426158
5 :  1.6515170754548774
6 :  1.4169634671112834
7 :  1.4218390547378612
8 :  1.4291212486183928
9 :  1.6219194916426156
10 :  1.6515170754548774
11 :  1.682651145927977
12 :  1.4145193713424764
13 :  1.4169634671112827
14 :  1.4218390547378605
15 :  1.4291212486183935
16 :  1.621919491642617
17 :  1.651517075454878
18 :  1.6826511459279758
19 :  1.7152380366850977
20 :  1.4169634671112827
21 :  1.4218390547378605
22 :  1.4291212486183922
23 :  1.621919491642618
24 :  1.6515170754548767
25 :  1.682651145927977
26 :  1.421839054737861
27 :  1.4291212486183926
28 :  1.6219194916426167
29 :  1.6515170754548774
30 :  1.429121248618393
31 :  1.621919491642616
32 :  1.4145193713424764
33 :  1.4145193713424764
34 :  1.4169634671112827
35 :  1.4145193713424764
36 :  1.4145193713424769
37 :  1.4169634671112823
38 :  1.421839054737861
39 :  1.4169634671112827
40 :  1.4145193713424764
41 :  1.4145193713424769
42 :  1.4169634671112827
43 :  1.4218390547378608
44 :  14.29121248618393
45 :  14.218390547378618
46 :  14.169634671112831
47 :  14.145193713424757
48 :  14.145193713424757
49 :  14.169634671112828
50 :  14.218390547378615
51 :  14.291212486183934
52 :  1.4218390547378612
53 :  1.4169634671112834
54 :  1.4145193713424764
55 :  1.4145193713424757
56 :  1.4169634671112827
57 :  1.4218390547378612
58 :  1.4169634671112827
59 :  1.4145193713424764
60 :  1.4145193713424775
61 :  1.4169634671112816
62 :  1.4145193713424764
63 :  1.4145193713424764
       
       
[(0, 1.4),
 (1, 1.4),
 (2, 1.4),
 (3, 1.4),
 (4, 1.4),
 (5, 1.4),
 (6, 1.4),
 (7, 1.4),
 (8, 1.4),
 (9, 1.4),
 (10, 1.4),
 (11, 1.4),
 (12, 1.4),
 (13, 1.4),
 (14, 1.4),
 (15, 1.4),
 (16, 1.4),
 (17, 1.4),
 (18, 1.4),
 (19, 1.4),
 (20, 1.4),
 (21, 1.4),
 (22, 1.4),
 (23, 1.4),
 (24, 1.4),
 (25, 1.4),
 (26, 1.4),
 (27, 1.4),
 (28, 1.4),
 (29, 1.4),
 (30, 1.4),
 (31, 1.4),
 (32, 1.4),
 (33, 1.4),
 (34, 1.4),
 (35, 1.4),
 (36, 1.4),
 (37, 1.4),
 (38, 1.4),
 (39, 1.4),
 (40, 1.4),
 (41, 1.4),
 (42, 1.4),
 (43, 1.4),
 (44, 14),
 (45, 14),
 (46, 14),
 (47, 14),
 (48, 14),
 (49, 14),
 (50, 14),
 (51, 14),
 (52, 1.4),
 (53, 1.4),
 (54, 1.4),
 (55, 1.4),
 (56, 1.4),
 (57, 1.4),
 (58, 1.4),
 (59, 1.4),
 (60, 1.4),
 (61, 1.4),
 (62, 1.4),
 (63, 1.4)]
[(0, 1.4),
 (1, 1.4),
 (2, 1.4),
 (3, 1.4),
 (4, 1.4),
 (5, 1.4),
 (6, 1.4),
 (7, 1.4),
 (8, 1.4),
 (9, 1.4),
 (10, 1.4),
 (11, 1.4),
 (12, 1.4),
 (13, 1.4),
 (14, 1.4),
 (15, 1.4),
 (16, 1.4),
 (17, 1.4),
 (18, 1.4),
 (19, 1.4),
 (20, 1.4),
 (21, 1.4),
 (22, 1.4),
 (23, 1.4),
 (24, 1.4),
 (25, 1.4),
 (26, 1.4),
 (27, 1.4),
 (28, 1.4),
 (29, 1.4),
 (30, 1.4),
 (31, 1.4),
 (32, 1.4),
 (33, 1.4),
 (34, 1.4),
 (35, 1.4),
 (36, 1.4),
 (37, 1.4),
 (38, 1.4),
 (39, 1.4),
 (40, 1.4),
 (41, 1.4),
 (42, 1.4),
 (43, 1.4),
 (44, 14),
 (45, 14),
 (46, 14),
 (47, 14),
 (48, 14),
 (49, 14),
 (50, 14),
 (51, 14),
 (52, 1.4),
 (53, 1.4),
 (54, 1.4),
 (55, 1.4),
 (56, 1.4),
 (57, 1.4),
 (58, 1.4),
 (59, 1.4),
 (60, 1.4),
 (61, 1.4),
 (62, 1.4),
 (63, 1.4)]
       
steps: 131
maximal force error:  9.975756378999279e-05
steps: 131
maximal force error:  9.975756378999279e-05
       
0 :  1.3999999995209256
1 :  1.3999999991377836
2 :  1.3999999774702634
3 :  1.4000000280032108
4 :  1.4000000312059218
5 :  1.3999999822479074
6 :  1.4000000220857656
7 :  1.4000000289190395
8 :  1.3999999447106728
9 :  1.399999960839987
10 :  1.4000000119919067
11 :  1.4000000161246315
12 :  1.3999999992687584
13 :  1.3999999953055775
14 :  1.4000000124257228
15 :  1.3999999558330503
16 :  1.3999999737569855
17 :  1.3999999998638324
18 :  1.399999998810049
19 :  1.3999999994949202
20 :  1.400000022085765
21 :  1.4000000289190402
22 :  1.3999999447106737
23 :  1.3999999608399858
24 :  1.4000000119919078
25 :  1.4000000161246309
26 :  1.399999977470263
27 :  1.4000000280032119
28 :  1.4000000312059198
29 :  1.3999999822479088
30 :  1.3999999995209254
31 :  1.3999999991377836
32 :  1.4000000003013142
33 :  1.400000000301314
34 :  1.4000024949801944
35 :  1.399997451442824
36 :  1.3999974514428246
37 :  1.4000024949801937
38 :  1.39999841412883
39 :  1.400006519468276
40 :  1.3999949909543579
41 :  1.3999949909543579
42 :  1.4000065194682771
43 :  1.3999984141288282
44 :  14.000002147839108
45 :  13.999968225618266
46 :  14.000099757563783
47 :  13.999928912972736
48 :  13.999928912972738
49 :  14.00009975756379
50 :  13.999968225618257
51 :  14.0000021478391
52 :  1.3999984009557924
53 :  1.4000065666868529
54 :  1.3999949569349976
55 :  1.399994956934999
56 :  1.4000065666868537
57 :  1.3999984009557902
58 :  1.40000252448358
59 :  1.3999974216654103
60 :  1.3999974216654134
61 :  1.4000025244835808
62 :  1.4000000002210904
63 :  1.4000000002210917
0 :  1.3999999995209256
1 :  1.3999999991377836
2 :  1.3999999774702634
3 :  1.4000000280032108
4 :  1.4000000312059218
5 :  1.3999999822479074
6 :  1.4000000220857656
7 :  1.4000000289190395
8 :  1.3999999447106728
9 :  1.399999960839987
10 :  1.4000000119919067
11 :  1.4000000161246315
12 :  1.3999999992687584
13 :  1.3999999953055775
14 :  1.4000000124257228
15 :  1.3999999558330503
16 :  1.3999999737569855
17 :  1.3999999998638324
18 :  1.399999998810049
19 :  1.3999999994949202
20 :  1.400000022085765
21 :  1.4000000289190402
22 :  1.3999999447106737
23 :  1.3999999608399858
24 :  1.4000000119919078
25 :  1.4000000161246309
26 :  1.399999977470263
27 :  1.4000000280032119
28 :  1.4000000312059198
29 :  1.3999999822479088
30 :  1.3999999995209254
31 :  1.3999999991377836
32 :  1.4000000003013142
33 :  1.400000000301314
34 :  1.4000024949801944
35 :  1.399997451442824
36 :  1.3999974514428246
37 :  1.4000024949801937
38 :  1.39999841412883
39 :  1.400006519468276
40 :  1.3999949909543579
41 :  1.3999949909543579
42 :  1.4000065194682771
43 :  1.3999984141288282
44 :  14.000002147839108
45 :  13.999968225618266
46 :  14.000099757563783
47 :  13.999928912972736
48 :  13.999928912972738
49 :  14.00009975756379
50 :  13.999968225618257
51 :  14.0000021478391
52 :  1.3999984009557924
53 :  1.4000065666868529
54 :  1.3999949569349976
55 :  1.399994956934999
56 :  1.4000065666868537
57 :  1.3999984009557902
58 :  1.40000252448358
59 :  1.3999974216654103
60 :  1.3999974216654134
61 :  1.4000025244835808
62 :  1.4000000002210904
63 :  1.4000000002210917
       
       
       

 

Primjer 2.:

 

       
0 :  [0, 0, 0]
1 :  []
2 :  []
3 :  []
4 :  [8, 0, 0]
5 :  []
6 :  []
7 :  []
8 :  []
9 :  []
10 :  []
11 :  []
12 :  []
13 :  []
14 :  []
15 :  []
16 :  []
17 :  []
18 :  []
19 :  []
20 :  []
21 :  []
22 :  []
23 :  []
24 :  []
25 :  []
26 :  []
27 :  []
28 :  []
29 :  []
30 :  []
31 :  []
32 :  []
33 :  []
34 :  []
35 :  []
36 :  [0, 8, 0]
37 :  []
38 :  []
39 :  []
40 :  [8, 8, 4]
0 :  [0, 0, 0]
1 :  []
2 :  []
3 :  []
4 :  [8, 0, 0]
5 :  []
6 :  []
7 :  []
8 :  []
9 :  []
10 :  []
11 :  []
12 :  []
13 :  []
14 :  []
15 :  []
16 :  []
17 :  []
18 :  []
19 :  []
20 :  []
21 :  []
22 :  []
23 :  []
24 :  []
25 :  []
26 :  []
27 :  []
28 :  []
29 :  []
30 :  []
31 :  []
32 :  []
33 :  []
34 :  []
35 :  []
36 :  [0, 8, 0]
37 :  []
38 :  []
39 :  []
40 :  [8, 8, 4]
       
[0, 4, 36, 40]
[0, 4, 36, 40]
       
0 :  [27, 32, 37]
1 :  [18, 23, 28, 33, 38]
2 :  [9, 14, 19, 24, 29, 34, 39]
3 :  [0, 5, 10, 15, 20, 25, 30, 35, 40]
4 :  [1, 6, 11, 16, 21, 26, 31]
5 :  [2, 7, 12, 17, 22]
6 :  [3, 8, 13]
7 :  [1, 5, 9]
8 :  [2, 6, 10, 14, 18]
9 :  [3, 7, 11, 15, 19, 23, 27]
10 :  [4, 8, 12, 16, 20, 24, 28, 32, 36]
11 :  [13, 17, 21, 25, 29, 33, 37]
12 :  [22, 26, 30, 34, 38]
13 :  [31, 35, 39]
14 :  [0, 1, 2, 3, 4]
15 :  [36, 37, 38, 39, 40]
16 :  [0, 9, 18, 27, 36]
17 :  [4, 13, 22, 31, 40]
0 :  [27, 32, 37]
1 :  [18, 23, 28, 33, 38]
2 :  [9, 14, 19, 24, 29, 34, 39]
3 :  [0, 5, 10, 15, 20, 25, 30, 35, 40]
4 :  [1, 6, 11, 16, 21, 26, 31]
5 :  [2, 7, 12, 17, 22]
6 :  [3, 8, 13]
7 :  [1, 5, 9]
8 :  [2, 6, 10, 14, 18]
9 :  [3, 7, 11, 15, 19, 23, 27]
10 :  [4, 8, 12, 16, 20, 24, 28, 32, 36]
11 :  [13, 17, 21, 25, 29, 33, 37]
12 :  [22, 26, 30, 34, 38]
13 :  [31, 35, 39]
14 :  [0, 1, 2, 3, 4]
15 :  [36, 37, 38, 39, 40]
16 :  [0, 9, 18, 27, 36]
17 :  [4, 13, 22, 31, 40]
       
       
0 :  [0, 1]
1 :  [2, 3, 4, 5]
2 :  [6, 7, 8, 9, 10, 11]
3 :  [12, 13, 14, 15, 16, 17, 18, 19]
4 :  [20, 21, 22, 23, 24, 25]
5 :  [26, 27, 28, 29]
6 :  [30, 31]
7 :  [32, 33]
8 :  [34, 35, 36, 37]
9 :  [38, 39, 40, 41, 42, 43]
10 :  [44, 45, 46, 47, 48, 49, 50, 51]
11 :  [52, 53, 54, 55, 56, 57]
12 :  [58, 59, 60, 61]
13 :  [62, 63]
14 :  [64, 65, 66, 67]
15 :  [68, 69, 70, 71]
16 :  [72, 73, 74, 75]
17 :  [76, 77, 78, 79]
0 :  [0, 1]
1 :  [2, 3, 4, 5]
2 :  [6, 7, 8, 9, 10, 11]
3 :  [12, 13, 14, 15, 16, 17, 18, 19]
4 :  [20, 21, 22, 23, 24, 25]
5 :  [26, 27, 28, 29]
6 :  [30, 31]
7 :  [32, 33]
8 :  [34, 35, 36, 37]
9 :  [38, 39, 40, 41, 42, 43]
10 :  [44, 45, 46, 47, 48, 49, 50, 51]
11 :  [52, 53, 54, 55, 56, 57]
12 :  [58, 59, 60, 61]
13 :  [62, 63]
14 :  [64, 65, 66, 67]
15 :  [68, 69, 70, 71]
16 :  [72, 73, 74, 75]
17 :  [76, 77, 78, 79]
       
       
       
       
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]
[14, 15, 16, 17]
[0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]
[14, 15, 16, 17]
       
[64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79]
[64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79]
       
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  6.25
15 :  6.25
16 :  6.25
17 :  6.25
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  6.25
15 :  6.25
16 :  6.25
17 :  6.25
       
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  1
15 :  1
16 :  1
17 :  1
18 :  1
19 :  1
20 :  1
21 :  1
22 :  1
23 :  1
24 :  1
25 :  1
26 :  1
27 :  1
28 :  1
29 :  1
30 :  1
31 :  1
32 :  1
33 :  1
34 :  1
35 :  1
36 :  1
37 :  1
38 :  1
39 :  1
40 :  1
41 :  1
42 :  1
43 :  1
44 :  1
45 :  1
46 :  1
47 :  1
48 :  1
49 :  1
50 :  1
51 :  1
52 :  1
53 :  1
54 :  1
55 :  1
56 :  1
57 :  1
58 :  1
59 :  1
60 :  1
61 :  1
62 :  1
63 :  1
64 :  6.25
65 :  6.25
66 :  6.25
67 :  6.25
68 :  6.25
69 :  6.25
70 :  6.25
71 :  6.25
72 :  6.25
73 :  6.25
74 :  6.25
75 :  6.25
76 :  6.25
77 :  6.25
78 :  6.25
79 :  6.25
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  1
15 :  1
16 :  1
17 :  1
18 :  1
19 :  1
20 :  1
21 :  1
22 :  1
23 :  1
24 :  1
25 :  1
26 :  1
27 :  1
28 :  1
29 :  1
30 :  1
31 :  1
32 :  1
33 :  1
34 :  1
35 :  1
36 :  1
37 :  1
38 :  1
39 :  1
40 :  1
41 :  1
42 :  1
43 :  1
44 :  1
45 :  1
46 :  1
47 :  1
48 :  1
49 :  1
50 :  1
51 :  1
52 :  1
53 :  1
54 :  1
55 :  1
56 :  1
57 :  1
58 :  1
59 :  1
60 :  1
61 :  1
62 :  1
63 :  1
64 :  6.25
65 :  6.25
66 :  6.25
67 :  6.25
68 :  6.25
69 :  6.25
70 :  6.25
71 :  6.25
72 :  6.25
73 :  6.25
74 :  6.25
75 :  6.25
76 :  6.25
77 :  6.25
78 :  6.25
79 :  6.25
       
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  1
15 :  1
16 :  1
17 :  1
18 :  1
19 :  1
20 :  1
21 :  1
22 :  1
23 :  1
24 :  1
25 :  1
26 :  1
27 :  1
28 :  1
29 :  1
30 :  1
31 :  1
32 :  1
33 :  1
34 :  1
35 :  1
36 :  1
37 :  1
38 :  1
39 :  1
40 :  1
41 :  1
42 :  1
43 :  1
44 :  1
45 :  1
46 :  1
47 :  1
48 :  1
49 :  1
50 :  1
51 :  1
52 :  1
53 :  1
54 :  1
55 :  1
56 :  1
57 :  1
58 :  1
59 :  1
60 :  1
61 :  1
62 :  1
63 :  1
64 :  6.25
65 :  6.25
66 :  6.25
67 :  6.25
68 :  6.25
69 :  6.25
70 :  6.25
71 :  6.25
72 :  6.25
73 :  6.25
74 :  6.25
75 :  6.25
76 :  6.25
77 :  6.25
78 :  6.25
79 :  6.25
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  1
15 :  1
16 :  1
17 :  1
18 :  1
19 :  1
20 :  1
21 :  1
22 :  1
23 :  1
24 :  1
25 :  1
26 :  1
27 :  1
28 :  1
29 :  1
30 :  1
31 :  1
32 :  1
33 :  1
34 :  1
35 :  1
36 :  1
37 :  1
38 :  1
39 :  1
40 :  1
41 :  1
42 :  1
43 :  1
44 :  1
45 :  1
46 :  1
47 :  1
48 :  1
49 :  1
50 :  1
51 :  1
52 :  1
53 :  1
54 :  1
55 :  1
56 :  1
57 :  1
58 :  1
59 :  1
60 :  1
61 :  1
62 :  1
63 :  1
64 :  6.25
65 :  6.25
66 :  6.25
67 :  6.25
68 :  6.25
69 :  6.25
70 :  6.25
71 :  6.25
72 :  6.25
73 :  6.25
74 :  6.25
75 :  6.25
76 :  6.25
77 :  6.25
78 :  6.25
79 :  6.25
       
       
       
0 :  1.1960272072907627
1 :  1.242602725914725
2 :  1.2728667778699367
3 :  1.2746643741090038
4 :  1.3162615786414187
5 :  1.3922963794239258
6 :  1.392204802828229
7 :  1.336535903473659
8 :  1.3307047282874076
9 :  1.3713077080859617
10 :  1.456446097481846
11 :  1.5884403267000098
12 :  1.6132130073131805
13 :  1.437314068108649
14 :  1.3627076146202934
15 :  1.3515652348139
16 :  1.392173707652228
17 :  1.4837401860206136
18 :  1.6388743028351174
19 :  1.9032121824517505
20 :  1.392204802828228
21 :  1.3365359034736586
22 :  1.3307047282874076
23 :  1.3713077080859613
24 :  1.4564460974818456
25 :  1.5884403267000107
26 :  1.2728667778699354
27 :  1.2746643741090031
28 :  1.3162615786414182
29 :  1.392296379423927
30 :  1.196027207290762
31 :  1.2426027259147263
32 :  1.1489114065001953
33 :  1.1489114065001946
34 :  1.2543742380154266
35 :  1.221081443384317
36 :  1.2210814433843176
37 :  1.2543742380154252
38 :  1.4065197656346569
39 :  1.3168919510964812
40 :  1.2737477454668196
41 :  1.273747745466819
42 :  1.316891951096481
43 :  1.4065197656346557
44 :  1.6699835980990725
45 :  1.4568582825123682
46 :  1.3444238704694962
47 :  1.2937122189659807
48 :  1.2937122189659807
49 :  1.3444238704694964
50 :  1.4568582825123673
51 :  1.669983598099072
52 :  1.417758975955832
53 :  1.3220540926977795
54 :  1.2744185384060047
55 :  1.2744185384060027
56 :  1.3220540926977813
57 :  1.4177589759558284
58 :  1.2663668432950266
59 :  1.2227456329729054
60 :  1.2227456329729078
61 :  1.2663668432950241
62 :  1.1528938884861937
63 :  1.1528938884861972
64 :  12.824882226452411
65 :  12.452687869221489
66 :  12.449626535755758
67 :  12.834083150198516
68 :  14.053191894999737
69 :  13.783982526133245
70 :  13.95905654885747
71 :  14.601141518874904
72 :  12.824882226452408
73 :  12.452687869221492
74 :  12.449626535755758
75 :  12.834083150198516
76 :  14.053191894999728
77 :  13.783982526133244
78 :  13.959056548857474
79 :  14.601141518874911
0 :  1.1960272072907627
1 :  1.242602725914725
2 :  1.2728667778699367
3 :  1.2746643741090038
4 :  1.3162615786414187
5 :  1.3922963794239258
6 :  1.392204802828229
7 :  1.336535903473659
8 :  1.3307047282874076
9 :  1.3713077080859617
10 :  1.456446097481846
11 :  1.5884403267000098
12 :  1.6132130073131805
13 :  1.437314068108649
14 :  1.3627076146202934
15 :  1.3515652348139
16 :  1.392173707652228
17 :  1.4837401860206136
18 :  1.6388743028351174
19 :  1.9032121824517505
20 :  1.392204802828228
21 :  1.3365359034736586
22 :  1.3307047282874076
23 :  1.3713077080859613
24 :  1.4564460974818456
25 :  1.5884403267000107
26 :  1.2728667778699354
27 :  1.2746643741090031
28 :  1.3162615786414182
29 :  1.392296379423927
30 :  1.196027207290762
31 :  1.2426027259147263
32 :  1.1489114065001953
33 :  1.1489114065001946
34 :  1.2543742380154266
35 :  1.221081443384317
36 :  1.2210814433843176
37 :  1.2543742380154252
38 :  1.4065197656346569
39 :  1.3168919510964812
40 :  1.2737477454668196
41 :  1.273747745466819
42 :  1.316891951096481
43 :  1.4065197656346557
44 :  1.6699835980990725
45 :  1.4568582825123682
46 :  1.3444238704694962
47 :  1.2937122189659807
48 :  1.2937122189659807
49 :  1.3444238704694964
50 :  1.4568582825123673
51 :  1.669983598099072
52 :  1.417758975955832
53 :  1.3220540926977795
54 :  1.2744185384060047
55 :  1.2744185384060027
56 :  1.3220540926977813
57 :  1.4177589759558284
58 :  1.2663668432950266
59 :  1.2227456329729054
60 :  1.2227456329729078
61 :  1.2663668432950241
62 :  1.1528938884861937
63 :  1.1528938884861972
64 :  12.824882226452411
65 :  12.452687869221489
66 :  12.449626535755758
67 :  12.834083150198516
68 :  14.053191894999737
69 :  13.783982526133245
70 :  13.95905654885747
71 :  14.601141518874904
72 :  12.824882226452408
73 :  12.452687869221492
74 :  12.449626535755758
75 :  12.834083150198516
76 :  14.053191894999728
77 :  13.783982526133244
78 :  13.959056548857474
79 :  14.601141518874911
       
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  6.25
15 :  6.25
16 :  6.25
17 :  6.25
0 :  1
1 :  1
2 :  1
3 :  1
4 :  1
5 :  1
6 :  1
7 :  1
8 :  1
9 :  1
10 :  1
11 :  1
12 :  1
13 :  1
14 :  6.25
15 :  6.25
16 :  6.25
17 :  6.25
       
       
[(0, 1),
 (1, 1),
 (2, 1),
 (3, 1),
 (4, 1),
 (5, 1),
 (6, 1),
 (7, 1),
 (8, 1),
 (9, 1),
 (10, 1),
 (11, 1),
 (12, 1),
 (13, 1),
 (14, 1),
 (15, 1),
 (16, 1),
 (17, 1),
 (18, 1),
 (19, 1),
 (20, 1),
 (21, 1),
 (22, 1),
 (23, 1),
 (24, 1),
 (25, 1),
 (26, 1),
 (27, 1),
 (28, 1),
 (29, 1),
 (30, 1),
 (31, 1),
 (32, 1),
 (33, 1),
 (34, 1),
 (35, 1),
 (36, 1),
 (37, 1),
 (38, 1),
 (39, 1),
 (40, 1),
 (41, 1),
 (42, 1),
 (43, 1),
 (44, 1),
 (45, 1),
 (46, 1),
 (47, 1),
 (48, 1),
 (49, 1),
 (50, 1),
 (51, 1),
 (52, 1),
 (53, 1),
 (54, 1),
 (55, 1),
 (56, 1),
 (57, 1),
 (58, 1),
 (59, 1),
 (60, 1),
 (61, 1),
 (62, 1),
 (63, 1),
 (64, 6.25),
 (65, 6.25),
 (66, 6.25),
 (67, 6.25),
 (68, 6.25),
 (69, 6.25),
 (70, 6.25),
 (71, 6.25),
 (72, 6.25),
 (73, 6.25),
 (74, 6.25),
 (75, 6.25),
 (76, 6.25),
 (77, 6.25),
 (78, 6.25),
 (79, 6.25)]
[(0, 1),
 (1, 1),
 (2, 1),
 (3, 1),
 (4, 1),
 (5, 1),
 (6, 1),
 (7, 1),
 (8, 1),
 (9, 1),
 (10, 1),
 (11, 1),
 (12, 1),
 (13, 1),
 (14, 1),
 (15, 1),
 (16, 1),
 (17, 1),
 (18, 1),
 (19, 1),
 (20, 1),
 (21, 1),
 (22, 1),
 (23, 1),
 (24, 1),
 (25, 1),
 (26, 1),
 (27, 1),
 (28, 1),
 (29, 1),
 (30, 1),
 (31, 1),
 (32, 1),
 (33, 1),
 (34, 1),
 (35, 1),
 (36, 1),
 (37, 1),
 (38, 1),
 (39, 1),
 (40, 1),
 (41, 1),
 (42, 1),
 (43, 1),
 (44, 1),
 (45, 1),
 (46, 1),
 (47, 1),
 (48, 1),
 (49, 1),
 (50, 1),
 (51, 1),
 (52, 1),
 (53, 1),
 (54, 1),
 (55, 1),
 (56, 1),
 (57, 1),
 (58, 1),
 (59, 1),
 (60, 1),
 (61, 1),
 (62, 1),
 (63, 1),
 (64, 6.25),
 (65, 6.25),
 (66, 6.25),
 (67, 6.25),
 (68, 6.25),
 (69, 6.25),
 (70, 6.25),
 (71, 6.25),
 (72, 6.25),
 (73, 6.25),
 (74, 6.25),
 (75, 6.25),
 (76, 6.25),
 (77, 6.25),
 (78, 6.25),
 (79, 6.25)]
       
WARNING: maximal number of steps (250) reached !
maximal force error:  0.264571921550413
WARNING: maximal number of steps (250) reached !
maximal force error:  0.264571921550413
       
       
       
       
0 :  0.9897244379261613
1 :  0.9986683237585273
2 :  0.9822151659092557
3 :  0.9660583690117344
4 :  0.970524344834898
5 :  0.9873068890641922
6 :  0.9991337293110453
7 :  0.9651953859671321
8 :  0.9528277854484988
9 :  0.9571603967366801
10 :  0.9739846402627143
11 :  1.0041506080111016
12 :  1.0360082985417987
13 :  0.9942719196830935
14 :  0.9617185826112437
15 :  0.9492403468272287
16 :  0.9535831705400167
17 :  0.9715527929174856
18 :  1.0044917714933101
19 :  1.0533088089851677
20 :  0.9991337293110452
21 :  0.9651953859671308
22 :  0.9528277854484981
23 :  0.957160396736681
24 :  0.9739846402627151
25 :  1.0041506080111022
26 :  0.9822151659092542
27 :  0.9660583690117341
28 :  0.9705243448348986
29 :  0.9873068890641944
30 :  0.9897244379261624
31 :  0.9986683237585267
32 :  0.9950013096882427
33 :  0.9950013096882429
34 :  0.985724474606509
35 :  0.9701947676295839
36 :  0.9701947676295839
37 :  0.9857244746065086
38 :  1.0008578395004628
39 :  0.9725078102192829
40 :  0.9591266639568565
41 :  0.9591266639568569
42 :  0.9725078102192826
43 :  1.0008578395004626
44 :  1.0375143372841695
45 :  1.0014453544590403
46 :  0.972535309572493
47 :  0.9585311695400429
48 :  0.9585311695400429
49 :  0.972535309572494
50 :  1.0014453544590394
51 :  1.037514337284169
52 :  1.0083171892261975
53 :  0.9795720508828246
54 :  0.9660054664754946
55 :  0.9660054664754933
56 :  0.9795720508828268
57 :  1.0083171892261968
58 :  1.0018367133116146
59 :  0.9852485753803142
60 :  0.985248575380313
61 :  1.0018367133116153
62 :  1.0182556780250136
63 :  1.0182556780250158
64 :  6.378486399349052
65 :  6.079964473784658
66 :  6.110489436624264
67 :  6.453316270239007
68 :  6.32512106586914
69 :  6.051128300507018
70 :  6.120625284183555
71 :  6.523884365656237
72 :  6.378486399349056
73 :  6.07996447378466
74 :  6.110489436624266
75 :  6.453316270238999
76 :  6.325121065869137
77 :  6.051128300507018
78 :  6.1206252841835544
79 :  6.523884365656237
0 :  0.9897244379261613
1 :  0.9986683237585273
2 :  0.9822151659092557
3 :  0.9660583690117344
4 :  0.970524344834898
5 :  0.9873068890641922
6 :  0.9991337293110453
7 :  0.9651953859671321
8 :  0.9528277854484988
9 :  0.9571603967366801
10 :  0.9739846402627143
11 :  1.0041506080111016
12 :  1.0360082985417987
13 :  0.9942719196830935
14 :  0.9617185826112437
15 :  0.9492403468272287
16 :  0.9535831705400167
17 :  0.9715527929174856
18 :  1.0044917714933101
19 :  1.0533088089851677
20 :  0.9991337293110452
21 :  0.9651953859671308
22 :  0.9528277854484981
23 :  0.957160396736681
24 :  0.9739846402627151
25 :  1.0041506080111022
26 :  0.9822151659092542
27 :  0.9660583690117341
28 :  0.9705243448348986
29 :  0.9873068890641944
30 :  0.9897244379261624
31 :  0.9986683237585267
32 :  0.9950013096882427
33 :  0.9950013096882429
34 :  0.985724474606509
35 :  0.9701947676295839
36 :  0.9701947676295839
37 :  0.9857244746065086
38 :  1.0008578395004628
39 :  0.9725078102192829
40 :  0.9591266639568565
41 :  0.9591266639568569
42 :  0.9725078102192826
43 :  1.0008578395004626
44 :  1.0375143372841695
45 :  1.0014453544590403
46 :  0.972535309572493
47 :  0.9585311695400429
48 :  0.9585311695400429
49 :  0.972535309572494
50 :  1.0014453544590394
51 :  1.037514337284169
52 :  1.0083171892261975
53 :  0.9795720508828246
54 :  0.9660054664754946
55 :  0.9660054664754933
56 :  0.9795720508828268
57 :  1.0083171892261968
58 :  1.0018367133116146
59 :  0.9852485753803142
60 :  0.985248575380313
61 :  1.0018367133116153
62 :  1.0182556780250136
63 :  1.0182556780250158
64 :  6.378486399349052
65 :  6.079964473784658
66 :  6.110489436624264
67 :  6.453316270239007
68 :  6.32512106586914
69 :  6.051128300507018
70 :  6.120625284183555
71 :  6.523884365656237
72 :  6.378486399349056
73 :  6.07996447378466
74 :  6.110489436624266
75 :  6.453316270238999
76 :  6.325121065869137
77 :  6.051128300507018
78 :  6.1206252841835544
79 :  6.523884365656237
       
[(64, 2.202309165299064),
 (65, 1.8526833278305133),
 (66, 1.8872415274385066),
 (67, 2.312144109667196),
 (68, 2.35255117319463),
 (69, 2.022928134339961),
 (70, 2.1323379574837067),
 (71, 2.7339675988485754),
 (72, 2.2023091652990634),
 (73, 1.8526833278305124),
 (74, 1.8872415274385077),
 (75, 2.312144109667196),
 (76, 2.352551173194629),
 (77, 2.022928134339963),
 (78, 2.1323379574837062),
 (79, 2.7339675988485745)]
[(64, 2.202309165299064),
 (65, 1.8526833278305133),
 (66, 1.8872415274385066),
 (67, 2.312144109667196),
 (68, 2.35255117319463),
 (69, 2.022928134339961),
 (70, 2.1323379574837067),
 (71, 2.7339675988485754),
 (72, 2.2023091652990634),
 (73, 1.8526833278305124),
 (74, 1.8872415274385077),
 (75, 2.312144109667196),
 (76, 2.352551173194629),
 (77, 2.022928134339963),
 (78, 2.1323379574837062),
 (79, 2.7339675988485745)]
       
[(0, 1),
 (1, 1),
 (2, 1),
 (3, 1),
 (4, 1),
 (5, 1),
 (6, 1),
 (7, 1),
 (8, 1),
 (9, 1),
 (10, 1),
 (11, 1),
 (12, 1),
 (13, 1),
 (14, 1),
 (15, 1),
 (16, 1),
 (17, 1),
 (18, 1),
 (19, 1),
 (20, 1),
 (21, 1),
 (22, 1),
 (23, 1),
 (24, 1),
 (25, 1),
 (26, 1),
 (27, 1),
 (28, 1),
 (29, 1),
 (30, 1),
 (31, 1),
 (32, 1),
 (33, 1),
 (34, 1),
 (35, 1),
 (36, 1),
 (37, 1),
 (38, 1),
 (39, 1),
 (40, 1),
 (41, 1),
 (42, 1),
 (43, 1),
 (44, 1),
 (45, 1),
 (46, 1),
 (47, 1),
 (48, 1),
 (49, 1),
 (50, 1),
 (51, 1),
 (52, 1),
 (53, 1),
 (54, 1),
 (55, 1),
 (56, 1),
 (57, 1),
 (58, 1),
 (59, 1),
 (60, 1),
 (61, 1),
 (62, 1),
 (63, 1)]
[(0, 1),
 (1, 1),
 (2, 1),
 (3, 1),
 (4, 1),
 (5, 1),
 (6, 1),
 (7, 1),
 (8, 1),
 (9, 1),
 (10, 1),
 (11, 1),
 (12, 1),
 (13, 1),
 (14, 1),
 (15, 1),
 (16, 1),
 (17, 1),
 (18, 1),
 (19, 1),
 (20, 1),
 (21, 1),
 (22, 1),
 (23, 1),
 (24, 1),
 (25, 1),
 (26, 1),
 (27, 1),
 (28, 1),
 (29, 1),
 (30, 1),
 (31, 1),
 (32, 1),
 (33, 1),
 (34, 1),
 (35, 1),
 (36, 1),
 (37, 1),
 (38, 1),
 (39, 1),
 (40, 1),
 (41, 1),
 (42, 1),
 (43, 1),
 (44, 1),
 (45, 1),
 (46, 1),
 (47, 1),
 (48, 1),
 (49, 1),
 (50, 1),
 (51, 1),
 (52, 1),
 (53, 1),
 (54, 1),
 (55, 1),
 (56, 1),
 (57, 1),
 (58, 1),
 (59, 1),
 (60, 1),
 (61, 1),
 (62, 1),
 (63, 1)]
       
[(0, 1),
 (1, 1),
 (2, 1),
 (3, 1),
 (4, 1),
 (5, 1),
 (6, 1),
 (7, 1),
 (8, 1),
 (9, 1),
 (10, 1),
 (11, 1),
 (12, 1),
 (13, 1),
 (14, 1),
 (15, 1),
 (16, 1),
 (17, 1),
 (18, 1),
 (19, 1),
 (20, 1),
 (21, 1),
 (22, 1),
 (23, 1),
 (24, 1),
 (25, 1),
 (26, 1),
 (27, 1),
 (28, 1),
 (29, 1),
 (30, 1),
 (31, 1),
 (32, 1),
 (33, 1),
 (34, 1),
 (35, 1),
 (36, 1),
 (37, 1),
 (38, 1),
 (39, 1),
 (40, 1),
 (41, 1),
 (42, 1),
 (43, 1),
 (44, 1),
 (45, 1),
 (46, 1),
 (47, 1),
 (48, 1),
 (49, 1),
 (50, 1),
 (51, 1),
 (52, 1),
 (53, 1),
 (54, 1),
 (55, 1),
 (56, 1),
 (57, 1),
 (58, 1),
 (59, 1),
 (60, 1),
 (61, 1),
 (62, 1),
 (63, 1)]
[(0, 1),
 (1, 1),
 (2, 1),
 (3, 1),
 (4, 1),
 (5, 1),
 (6, 1),
 (7, 1),
 (8, 1),
 (9, 1),
 (10, 1),
 (11, 1),
 (12, 1),
 (13, 1),
 (14, 1),
 (15, 1),
 (16, 1),
 (17, 1),
 (18, 1),
 (19, 1),
 (20, 1),
 (21, 1),
 (22, 1),
 (23, 1),
 (24, 1),
 (25, 1),
 (26, 1),
 (27, 1),
 (28, 1),
 (29, 1),
 (30, 1),
 (31, 1),
 (32, 1),
 (33, 1),
 (34, 1),
 (35, 1),
 (36, 1),
 (37, 1),
 (38, 1),
 (39, 1),
 (40, 1),
 (41, 1),
 (42, 1),
 (43, 1),
 (44, 1),
 (45, 1),
 (46, 1),
 (47, 1),
 (48, 1),
 (49, 1),
 (50, 1),
 (51, 1),
 (52, 1),
 (53, 1),
 (54, 1),
 (55, 1),
 (56, 1),
 (57, 1),
 (58, 1),
 (59, 1),
 (60, 1),
 (61, 1),
 (62, 1),
 (63, 1)]
       
steps: 116
maximal force error:  9.546959981943104e-05
maximal length error:  0.00017164170608952034
steps: 116
maximal force error:  9.546959981943104e-05
maximal length error:  0.00017164170608952034