《有限元方法编程》程序4.3结果问题

piteqiu

我按照程序调试，结果如下：
There are    12  equations and the half-bandwidth is     3

   Node   Displacement Moment
    1     0.9031E+00 -0.5722E-03
    2     0.6617E+00 -0.1794E+01
    3     0.4206E+00 -0.2595E+01
    4     0.1800E+00 -0.2114E+01
    5    -0.6012E-01 -0.8324E+00
    6    -0.3001E+00 -0.1088E-03
======================================================
书中结果如下：
There are    12  equations and the half-bandwidth is     3

   Node   Displacement Moment
    1     0.9031E+00 -0.2961E-03
    2     0.6617E+00 -0.1794E+01
    3     0.4206E+00 -0.2595E+01
    4     0.1800E+00 -0.2114E+01
    5    -0.6012E-01 -0.8324E+00
    6    -0.3001E+00  0.2977E-03

通过比较发现，在第一个点和最后一个点的弯矩不一样。但是数据文件和源程序是一样的，却得出不同结果。

哪位同学算过这个算例的出来说说是怎么回事吧。

孤叶

软件版本的原因