|
|
代码：
call expectek(b1, b2, b, helenev, expectk)
do i=1,nnode
call ghmatrix(nelem, nnode, ek)
enddo
do i=1,nnode
expect0(i)= expect0(i)+bcp(i)*expectk
write(13,*) i, expect0(i)
enddo
有关expect0(i)的输出文件
1 (16.7373554191698,63.5228987726871)
2 (16.7373554191700,63.5228987726868)
3 (16.7373554191701,63.5228987726865)
4 (16.7373554191696,63.5228987726868)
5 (-69.1054022723600,187.713203864395)
6 (-69.1054022723600,187.713203864395)
7 (16.7373554191698,63.5228987726866)
8 (-69.1054022723595,187.713203864395)
9 (16.7373554191697,63.5228987726868)
10 (-69.1054022723597,187.713203864395)
11 (16.7373554191696,63.5228987726872)
12 (16.7373554191696,63.5228987726870)
13 (16.7373554191695,63.5228987726871)
14 (16.7373554191698,63.5228987726868)
15 (16.7373554191695,63.5228987726873)
16 (-69.1054022723598,187.713203864395)
17 (16.7373554191697,63.5228987726868)
18 (16.7373554191697,63.5228987726868)
19 (16.7373554191701,63.5228987726865)
20 (16.7373554191697,63.5228987726873)
21 (16.7373554191700,63.5228987726869)
22 (16.7373554191697,63.5228987726874)
23 (16.7373554191700,63.5228987726868)
24 (-69.1054022723596,187.713203864395)
25 (-69.1054022723596,187.713203864395)
26 (16.7373554191701,63.5228987726867)
27 (16.7373554191699,63.5228987726868)
28 (16.7373554191705,63.5228987726864)
29 (16.7373554191695,63.5228987726872)
30 (-69.1054022723602,187.713203864396)
31 (16.7373554191699,63.5228987726868)
32 (16.7373554191698,63.5228987726868)
33 (16.7373554191690,63.5228987726874)
34 (16.7373554191704,63.5228987726863)
35 (16.7373554191695,63.5228987726869)
36 (16.7373554191692,63.5228987726873)
37 (16.7373554191704,63.5228987726861)
39 (16.7373554191692,63.5228987726873)
40 (16.7373554191696,63.5228987726867)
41 (16.7373554191699,63.5228987726867)
42 (16.7373554191703,63.5228987726862)
43 (-69.1054022723594,187.713203864395)
44 (16.7373554191695,63.5228987726870)
45 (-69.1054022723596,187.713203864395)
46 (16.7373554191698,63.5228987726871)
47 (16.7373554191701,63.5228987726870)
48 (16.7373554191708,63.5228987726860)
49 (16.7373554191698,63.5228987726871)
50 (16.7373554191701,63.5228987726870)
51 (16.7373554191702,63.5228987726867)
52 (-69.1054022723594,187.713203864395)
53 (16.7373554191698,63.5228987726869)
54 (16.7373554191689,63.5228987726879)
55 (-69.1054022723600,187.713203864396)
56 (16.7373554191699,63.5228987726870)
57 (16.7373554191701,63.5228987726866)
58 (16.7373554191700,63.5228987726867)
59 (16.7373554191691,63.5228987726875)
60 (16.7373554191696,63.5228987726871)
61 (-69.1054022723600,187.713203864395)
62 (-69.1054022723598,187.713203864395)
63 (-69.1054022723598,187.713203864395)
64 (-69.1054022723599,187.713203864395)
65 (16.7373554191696,63.5228987726867)
66 (16.7373554191684,63.5228987726879)
67 (-69.1054022723597,187.713203864395)
68 (16.7373554191709,63.5228987726856)
69 (16.7373554191697,63.5228987726867)
70 (-69.1054022723601,187.713203864395)
71 (16.7373554191696,63.5228987726871)
72 (16.7373554191682,63.5228987726883)
73 (16.7373554191700,63.5228987726868)
74 (-69.1054022723596,187.713203864395)
75 (16.7373554191701,63.5228987726868)
76 (16.7373554191698,63.5228987726869)
77 (-69.1054022723597,187.713203864395)
78 (16.7373554191699,63.5228987726870)
79 (16.7373554191698,63.5228987726873)
80 (-69.1054022723600,187.713203864395)
81 (16.7373554191710,63.5228987726859)
82 (16.7373554191685,63.5228987726877)
83 (-69.1054022723595,187.713203864395)
84 (-69.1054022723598,187.713203864395)
85 (16.7373554191696,63.5228987726867)
86 (-69.1054022723600,187.713203864395)
87 (16.7373554191698,63.5228987726870)
88 (16.7373554191684,63.5228987726879)
89 (16.7373554191698,63.5228987726868)
90 (-69.1054022723598,187.713203864395)
91 (16.7373554191697,63.5228987726869)
92 (-69.1054022723597,187.713203864395)
93 (16.7373554191700,63.5228987726866)
94 (-69.1054022723596,187.713203864395)
95 (16.7373554191697,63.5228987726870)
96 (16.7373554191710,63.5228987726862)
97 (16.7373554191695,63.5228987726874)
98 (16.7373554191686,63.5228987726881)
99 (16.7373554191701,63.5228987726871)
100 (16.7373554193436,63.5228987725356)
1 (16.7373554191698,63.5228987726871)
2 (16.7373554191700,63.5228987726868)
3 (16.7373554191701,63.5228987726865)
4 (16.7373554191696,63.5228987726868)
5 (-69.1054022723600,187.713203864395)
6 (-69.1054022723600,187.713203864395)
7 (16.7373554191698,63.5228987726866)
8 (-69.1054022723595,187.713203864395)
9 (16.7373554191697,63.5228987726868)
10 (-69.1054022723597,187.713203864395)
11 (16.7373554191696,63.5228987726872)
12 (16.7373554191696,63.5228987726870)
13 (16.7373554191695,63.5228987726871)
14 (16.7373554191698,63.5228987726868)
15 (16.7373554191695,63.5228987726873)
16 (-69.1054022723598,187.713203864395)
17 (16.7373554191697,63.5228987726868)
18 (16.7373554191697,63.5228987726868)
19 (16.7373554191701,63.5228987726865)
20 (16.7373554191697,63.5228987726873)
21 (16.7373554191700,63.5228987726869)
22 (16.7373554191697,63.5228987726874)
23 (16.7373554191700,63.5228987726868)
24 (-69.1054022723596,187.713203864395)
25 (-69.1054022723596,187.713203864395)
26 (16.7373554191701,63.5228987726867)
27 (16.7373554191699,63.5228987726868)
28 (16.7373554191705,63.5228987726864)
29 (16.7373554191695,63.5228987726872)
30 (-69.1054022723602,187.713203864396)
31 (16.7373554191699,63.5228987726868)
32 (16.7373554191698,63.5228987726868)
33 (16.7373554191690,63.5228987726874)
34 (16.7373554191704,63.5228987726863)
35 (16.7373554191695,63.5228987726869)
36 (16.7373554191692,63.5228987726873)
37 (16.7373554191704,63.5228987726861)
38 (16.7373554191698,63.5228987726866)
39 (16.7373554191692,63.5228987726873)
40 (16.7373554191696,63.5228987726867)
41 (16.7373554191699,63.5228987726867)
42 (16.7373554191703,63.5228987726862)
43 (-69.1054022723594,187.713203864395)
44 (16.7373554191695,63.5228987726870)
45 (-69.1054022723596,187.713203864395)
46 (16.7373554191698,63.5228987726871)
47 (16.7373554191701,63.5228987726870)
48 (16.7373554191708,63.5228987726860)
49 (16.7373554191698,63.5228987726871)
50 (16.7373554191701,63.5228987726870)
51 (16.7373554191702,63.5228987726867)
52 (-69.1054022723594,187.713203864395)
53 (16.7373554191698,63.5228987726869)
54 (16.7373554191689,63.5228987726879)
55 (-69.1054022723600,187.713203864396)
56 (16.7373554191699,63.5228987726870)
57 (16.7373554191701,63.5228987726866)
58 (16.7373554191700,63.5228987726867)
59 (16.7373554191691,63.5228987726875)
60 (16.7373554191696,63.5228987726871)
61 (-69.1054022723600,187.713203864395)
62 (-69.1054022723598,187.713203864395)
63 (-69.1054022723598,187.713203864395)
64 (-69.1054022723599,187.713203864395)
65 (16.7373554191696,63.5228987726867)
66 (16.7373554191684,63.5228987726879)
67 (-69.1054022723597,187.713203864395)
68 (16.7373554191709,63.5228987726856)
69 (16.7373554191697,63.5228987726867)
70 (-69.1054022723601,187.713203864395)
71 (16.7373554191696,63.5228987726871)
72 (16.7373554191682,63.5228987726883)
73 (16.7373554191700,63.5228987726868)
74 (-69.1054022723596,187.713203864395)
75 (16.7373554191701,63.5228987726868)
76 (16.7373554191698,63.5228987726869)
77 (-69.1054022723597,187.713203864395)
78 (16.7373554191699,63.5228987726870)
79 (16.7373554191698,63.5228987726873)
80 (-69.1054022723600,187.713203864395)
81 (16.7373554191710,63.5228987726859)
82 (16.7373554191685,63.5228987726877)
83 (-69.1054022723595,187.713203864395)
84 (-69.1054022723598,187.713203864395)
85 (16.7373554191696,63.5228987726867)
86 (-69.1054022723600,187.713203864395)
87 (16.7373554191698,63.5228987726870)
88 (16.7373554191684,63.5228987726879)
89 (16.7373554191698,63.5228987726868)
90 (-69.1054022723598,187.713203864395)
91 (16.7373554191697,63.5228987726869)
92 (-69.1054022723597,187.713203864395)
93 (16.7373554191700,63.5228987726866)
94 (-69.1054022723596,187.713203864395)
95 (16.7373554191697,63.5228987726870)
96 (16.7373554191710,63.5228987726862)
97 (16.7373554191695,63.5228987726874)
98 (16.7373554191686,63.5228987726881)
99 (16.7373554191701,63.5228987726871)
100 (16.7373554193436,63.5228987725356)
输出文件重复两次出现expect0(i)的值，不知道是哪的问题
|
|
F 币
贡献
权杖
窥视卡
发表于 2020-5-7 15:14:31
微信
收藏
点赞
点踩
提升卡
置顶卡
变色卡
千斤顶
显身卡
发表于 2020-5-7 16:08:25
楼主
发表于 2020-5-8 08:44:34