| 
 | 
 
代码: 
    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)的值,不知道是哪的问题 
 
 |   
 
 
 
 |