p102.res

I[1].M.Smith所著的《有限元方法编程》第三版Fortran程序

Global coordinates
Node    1      0.0000E+00  0.0000E+00
Node    2      0.0000E+00 -0.5000E+00
Node    3      0.0000E+00 -0.1000E+01
Node    4      0.4000E+00  0.0000E+00
Node    5      0.4000E+00 -0.1000E+01
Node    6      0.8000E+00  0.0000E+00
Node    7      0.8000E+00 -0.5000E+00
Node    8      0.8000E+00 -0.1000E+01
Node    9      0.1200E+01  0.0000E+00
Node   10      0.1200E+01 -0.1000E+01
Node   11      0.1600E+01  0.0000E+00
Node   12      0.1600E+01 -0.5000E+00
Node   13      0.1600E+01 -0.1000E+01
Node   14      0.2000E+01  0.0000E+00
Node   15      0.2000E+01 -0.1000E+01
Node   16      0.2400E+01  0.0000E+00
Node   17      0.2400E+01 -0.5000E+00
Node   18      0.2400E+01 -0.1000E+01
Node   19      0.2800E+01  0.0000E+00
Node   20      0.2800E+01 -0.1000E+01
Node   21      0.3200E+01  0.0000E+00
Node   22      0.3200E+01 -0.5000E+00
Node   23      0.3200E+01 -0.1000E+01
Node   24      0.3600E+01  0.0000E+00
Node   25      0.3600E+01 -0.1000E+01
Node   26      0.4000E+01  0.0000E+00
Node   27      0.4000E+01 -0.5000E+00
Node   28      0.4000E+01 -0.1000E+01
Global node numbers
Element    1          3    2    1    4    6    7    8    5
Element    2          8    7    6    9   11   12   13   10
Element    3         13   12   11   14   16   17   18   15
Element    4         18   17   16   19   21   22   23   20
Element    5         23   22   21   24   26   27   28   25
There are    50  equations and the half-bandwidth is   15
The global mass diagonal is:
  0.1600E+00  0.1600E+00  0.1600E+00  0.1600E+00  0.8000E-01  0.8000E-01
  0.3200E+00  0.3200E+00  0.8000E-01  0.8000E-01  0.1600E+00  0.1600E+00
  0.1600E+00  0.1600E+00  0.8000E-01  0.8000E-01  0.3200E+00  0.3200E+00
  0.8000E-01  0.8000E-01  0.1600E+00  0.1600E+00  0.1600E+00  0.1600E+00
  0.8000E-01  0.8000E-01  0.3200E+00  0.3200E+00  0.8000E-01  0.8000E-01
  0.1600E+00  0.1600E+00  0.1600E+00  0.1600E+00  0.8000E-01  0.8000E-01
  0.3200E+00  0.3200E+00  0.8000E-01  0.8000E-01  0.1600E+00  0.1600E+00
  0.1600E+00  0.1600E+00  0.4000E-01  0.4000E-01  0.1600E+00  0.1600E+00
  0.4000E-01  0.4000E-01
The eigenvalues are:
  0.4083E-02  0.9840E-01  0.1706E+00  0.5002E+00  0.1240E+01  0.1337E+01
  0.1606E+01  0.2232E+01  0.3035E+01  0.3071E+01  0.3617E+01  0.4170E+01
  0.4351E+01  0.5301E+01  0.5505E+01  0.6241E+01  0.6451E+01  0.6959E+01
  0.7002E+01  0.8034E+01  0.8105E+01  0.8305E+01  0.8730E+01  0.8908E+01
  0.9296E+01  0.9747E+01  0.9865E+01  0.1021E+02  0.1043E+02  0.1063E+02
  0.1114E+02  0.1193E+02  0.1223E+02  0.1313E+02  0.1452E+02  0.1611E+02
  0.2046E+02  0.2224E+02  0.2460E+02  0.2610E+02  0.2764E+02  0.3194E+02
  0.3224E+02  0.4417E+02  0.4521E+02  0.5879E+02  0.6235E+02  0.6759E+02
  0.8007E+02  0.9102E+02
Eigenvector number   1 is: 
  0.3450E-01 -0.2445E-01 -0.3450E-01 -0.2445E-01  0.6826E-01 -0.7617E-01
  0.4841E-13 -0.6472E-01 -0.6826E-01 -0.7617E-01  0.9816E-01 -0.1475E+00
 -0.9816E-01 -0.1475E+00  0.1212E+00 -0.2398E+00  0.9336E-13 -0.2352E+00
 -0.1212E+00 -0.2398E+00  0.1384E+00 -0.3493E+00 -0.1384E+00 -0.3493E+00
  0.1511E+00 -0.4704E+00  0.7832E-13 -0.4674E+00 -0.1511E+00 -0.4704E+00
  0.1595E+00 -0.5985E+00 -0.1595E+00 -0.5985E+00  0.1643E+00 -0.7314E+00
  0.8503E-13 -0.7305E+00 -0.1643E+00 -0.7314E+00  0.1664E+00 -0.8656E+00
 -0.1664E+00 -0.8656E+00  0.1671E+00 -0.1000E+01  0.8154E-13 -0.9999E+00
 -0.1671E+00 -0.1000E+01
Eigenvector number   2 is: 
  0.1321E+00 -0.1834E+00 -0.1321E+00 -0.1834E+00  0.1872E+00 -0.4372E+00
 -0.5371E-14 -0.4192E+00 -0.1872E+00 -0.4372E+00  0.1606E+00 -0.6530E+00
 -0.1606E+00 -0.6530E+00  0.5414E-01 -0.7995E+00 -0.1042E-13 -0.8395E+00
 -0.5414E-01 -0.7995E+00 -0.9898E-01 -0.8097E+00  0.9898E-01 -0.8097E+00
 -0.2631E+00 -0.6731E+00 -0.1093E-13 -0.7150E+00  0.2631E+00 -0.6731E+00
 -0.4104E+00 -0.3761E+00  0.4104E+00 -0.3761E+00 -0.5182E+00  0.3220E-01
 -0.1236E-13  0.8696E-02  0.5182E+00  0.3220E-01 -0.5800E+00  0.5061E+00
  0.5800E+00  0.5061E+00 -0.5995E+00  0.1000E+01 -0.1191E-13  0.9983E+00
  0.5995E+00  0.1000E+01
Eigenvector number   3 is: 
  0.1428E+00 -0.7847E-01  0.1428E+00  0.7847E-01  0.2959E+00 -0.6998E-01
  0.2830E+00 -0.1783E-14  0.2959E+00  0.6998E-01  0.4308E+00 -0.7509E-01
  0.4308E+00  0.7509E-01  0.5598E+00 -0.6559E-01  0.5560E+00 -0.7784E-14
  0.5598E+00  0.6559E-01  0.6767E+00 -0.5690E-01  0.6767E+00  0.5690E-01
  0.7825E+00 -0.4875E-01  0.7659E+00 -0.1094E-13  0.7825E+00  0.4875E-01
  0.8539E+00 -0.4031E-01  0.8539E+00  0.4031E-01  0.9111E+00 -0.2711E-01
  0.9113E+00 -0.6111E-14  0.9111E+00  0.2711E-01  0.9600E+00 -0.5175E-02
  0.9600E+00  0.5175E-02  0.1000E+01 -0.4300E-01  0.9408E+00  0.3155E-14
  0.1000E+01  0.4300E-01
Eigenvector number   4 is: 
  0.1908E+00 -0.4993E+00 -0.1908E+00 -0.4993E+00  0.9003E-01 -0.9723E+00
  0.5804E-15 -0.1000E+01 -0.9003E-01 -0.9723E+00 -0.1767E+00 -0.9975E+00
  0.1767E+00 -0.9975E+00 -0.4435E+00 -0.6848E+00  0.7939E-15 -0.7485E+00
  0.4435E+00 -0.6848E+00 -0.5330E+00 -0.6995E-01  0.5330E+00 -0.6995E-01
 -0.3336E+00  0.5097E+00  0.1126E-13  0.6005E+00  0.3336E+00  0.5097E+00
  0.6264E-01  0.7222E+00 -0.6264E-01  0.7222E+00  0.5046E+00  0.5093E+00
  0.1123E-13  0.6246E+00 -0.5046E+00  0.5093E+00  0.8327E+00 -0.1429E+00
 -0.8327E+00 -0.1429E+00  0.9446E+00 -0.9532E+00  0.1572E-13 -0.9474E+00
 -0.9446E+00 -0.9532E+00
Eigenvector number   5 is: 
  0.4024E-01 -0.6070E+00 -0.4024E-01 -0.6070E+00 -0.3088E+00 -0.9026E+00
  0.4105E-12 -0.1000E+01  0.3088E+00 -0.9026E+00 -0.6171E+00 -0.3276E+00
  0.6171E+00 -0.3276E+00 -0.5013E+00  0.4258E+00  0.5099E-12  0.5117E+00
  0.5013E+00  0.4258E+00 -0.9460E-01  0.6869E+00  0.9460E-01  0.6869E+00
  0.2750E+00  0.4249E+00  0.2069E-12  0.4992E+00 -0.2750E+00  0.4249E+00
  0.2845E+00 -0.2511E+00 -0.2845E+00 -0.2511E+00 -0.1927E+00 -0.6623E+00
 -0.4585E-12 -0.8507E+00  0.1927E+00 -0.6623E+00 -0.7458E+00 -0.2447E+00
  0.7458E+00 -0.2447E+00 -0.9514E+00  0.5520E+00 -0.2074E-12  0.5574E+00
  0.9514E+00  0.5520E+00