integer lwork, lrwork, liwork, info parameter ( lwork = 5000, lrwork = 5000, liwork = 500 ) integer iwork(liwork) double complex aa(40, 40), work(lwork) double precision w(40), rwork(lrwork) C data aa / C # (-0.435191, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-6.45929, -0.264506), C # (0.000654894, 0.0494788), C # (0.0010041, 0.0458484), C # (0.0150791, -0.372499), C # (-0.012103, 0.30548), C # (-0.000274567, -0.0405584), C # (-0.00056093, -0.0375814), C # (1.01677d-08, -1.05704d-07), C # (-4.79086d-07, 1.41263d-05), C # (-7.9517, -0.32562), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (6.3853, 0), C # (0, 0.0408), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000555262, -0.0419514), C # (1.98856, 0.0814425), C # (0.0170135, -0.42942), C # (0.000788511, 0.052829), C # (1.21157d-05, -0.00542571), C # (-5.57524, -0.228304), C # (-0.00165939, 0.0442559), C # (0.000264062, 0.0348546), C # (-0.000151981, -0.0201374), C # (-0.000486679, -0.0367698), C # (0, 0), C # (0, 0), C # (0, 0), C # (0.0408, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, -0.0408), C # (6.3853, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000851343, -0.0388733), C # (0.0170135, -0.42942), C # (1.98856, 0.0814422), C # (0.000385964, 0.0570138), C # (5.36068d-05, -0.00585705), C # (-0.00165937, 0.0442557), C # (-5.57524, -0.228304), C # (-0.00051019, -0.0322958), C # (-0.000294071, -0.0186603), C # (-0.000746189, -0.0340719), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0.0408), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (6.3853, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.0127851, 0.315829), C # (0.000788511, 0.052829), C # (0.000385964, 0.0570138), C # (1.98856, 0.0814079), C # (-5.57524, -0.228304), C # (1.29593d-05, -0.00545039), C # (5.443d-05, -0.00588152), C # (-1.02116d-08, 1.06159d-07), C # (0.0120437, -0.303168), C # (-0.0112059, 0.276819), C # (0, 0), C # (-0.0408, 0), C # (0, -0.0408), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.0137295, 0.346532), C # (2.80669d-05, 0.00675756), C # (-2.35543d-05, 0.00729422), C # (5.10626, 0.2091), C # (-3.07027, -0.125729), C # (-0.000372858, -0.00503724), C # (-0.000334615, -0.00543482), C # (-1.15417d-08, 1.19988d-07), C # (-0.0142558, 0.353289), C # (-0.00346891, 0.0875554), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000311465, -0.0460089), C # (5.10626, 0.2091), C # (0.00210793, -0.0550618), C # (2.72943d-05, 0.00678016), C # (-0.000372858, -0.00503724), C # (-3.07027, -0.125726), C # (-0.00191428, 0.0407779), C # (-0.000498236, -0.0406443), C # (0.000288194, 0.02345), C # (-7.86952d-05, -0.0116247), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000636311, -0.0426318), C # (0.00210794, -0.055062), C # (5.10626, 0.2091), C # (-2.43082d-05, 0.00731662), C # (-0.000334615, -0.00543482), C # (-0.00191428, 0.0407779), C # (-3.07027, -0.125726), C # (0.000785054, 0.0376625), C # (0.000453802, 0.0217284), C # (-0.000160771, -0.0107714), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (1.15341d-08, -1.19909d-07), C # (-0.000216307, -0.0299928), C # (0.00042811, 0.0277909), C # (9.35258d-09, -9.72294d-08), C # (-1.15417d-08, 1.19988d-07), C # (-0.000498236, -0.0406443), C # (0.000785054, 0.0376625), C # (3.23091, 0.132321), C # (1.17588d-08, -1.22245d-07), C # (2.91422d-09, -3.02963d-08), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (-5.43468d-07, 1.60246d-05), C # (0.000124449, 0.0173293), C # (0.000246723, 0.0160582), C # (-0.0103514, 0.26089), C # (-0.0142558, 0.353289), C # (0.000288194, 0.02345), C # (0.000453802, 0.0217284), C # (1.17588d-08, -1.22245d-07), C # (3.23091, 0.132288), C # (-1.37313d-07, 4.04881d-06), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (19.422, 0), C # (-6.2442, -0.255698), C # (0.000492274, 0.0371924), C # (0.000754767, 0.0344636), C # (0.0113347, -0.280002), C # (-0.00325794, 0.0822305), C # (-7.39092d-05, -0.0109177), C # (-0.000150994, -0.0101163), C # (2.73699d-09, -2.84538d-08), C # (-1.28963d-07, 3.80257d-06), C # (-14.6692, -0.600698), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-6.45929, 0.264506), C # (-0.000555262, 0.0419514), C # (-0.000851343, 0.0388733), C # (-0.0127851, -0.315829), C # (-0.0137295, -0.346532), C # (-0.000311465, 0.0460089), C # (-0.000636311, 0.0426318), C # (1.15341d-08, 1.19909d-07), C # (-5.43468d-07, -1.60246d-05), C # (-6.2442, 0.255698), C # (-5.9819, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0.000654894, -0.0494788), C # (1.98856, -0.0814425), C # (0.0170135, 0.42942), C # (0.000788511, -0.052829), C # (2.80669d-05, -0.00675756), C # (5.10626, -0.2091), C # (0.00210794, 0.055062), C # (-0.000216307, 0.0299928), C # (0.000124449, -0.0173293), C # (0.000492274, -0.0371924), C # (0, 0), C # (3.74537, 0), C # (0, 0.1824), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0.1824, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0.0010041, -0.0458484), C # (0.0170135, 0.42942), C # (1.98856, -0.0814422), C # (0.000385964, -0.0570138), C # (-2.35543d-05, -0.00729422), C # (0.00210793, 0.0550618), C # (5.10626, -0.2091), C # (0.00042811, -0.0277909), C # (0.000246723, -0.0160582), C # (0.000754767, -0.0344636), C # (0, 0), C # (0, -0.1824), C # (3.74537, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0.1824), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0.0150791, 0.372499), C # (0.000788511, -0.052829), C # (0.000385964, -0.0570138), C # (1.98856, -0.0814079), C # (5.10626, -0.2091), C # (2.72943d-05, -0.00678016), C # (-2.43082d-05, -0.00731662), C # (9.35258d-09, 9.72294d-08), C # (-0.0103514, -0.26089), C # (0.0113347, 0.280002), C # (0, 0), C # (0, 0), C # (0, 0), C # (3.74537, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.1824, 0), C # (0, -0.1824), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.012103, -0.30548), C # (1.21157d-05, 0.00542571), C # (5.36068d-05, 0.00585705), C # (-5.57524, 0.228304), C # (-3.07027, 0.125729), C # (-0.000372858, 0.00503724), C # (-0.000334615, 0.00543482), C # (-1.15417d-08, -1.19988d-07), C # (-0.0142558, -0.353289), C # (-0.00325794, -0.0822305), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000274567, 0.0405584), C # (-5.57524, 0.228304), C # (-0.00165937, -0.0442557), C # (1.29593d-05, 0.00545039), C # (-0.000372858, 0.00503724), C # (-3.07027, 0.125726), C # (-0.00191428, -0.0407779), C # (-0.000498236, 0.0406443), C # (0.000288194, -0.02345), C # (-7.39092d-05, 0.0109177), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.00056093, 0.0375814), C # (-0.00165939, -0.0442559), C # (-5.57524, 0.228304), C # (5.443d-05, 0.00588152), C # (-0.000334615, 0.00543482), C # (-0.00191428, -0.0407779), C # (-3.07027, 0.125726), C # (0.000785054, -0.0376625), C # (0.000453802, -0.0217284), C # (-0.000150994, 0.0101163), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (1.01677d-08, 1.05704d-07), C # (0.000264062, -0.0348546), C # (-0.00051019, 0.0322958), C # (-1.02116d-08, -1.06159d-07), C # (-1.15417d-08, -1.19988d-07), C # (-0.000498236, 0.0406443), C # (0.000785054, -0.0376625), C # (3.23091, -0.132321), C # (1.17588d-08, 1.22245d-07), C # (2.73699d-09, 2.84538d-08), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-4.79086d-07, -1.41263d-05), C # (-0.000151981, 0.0201374), C # (-0.000294071, 0.0186603), C # (0.0120437, 0.303168), C # (-0.0142558, -0.353289), C # (0.000288194, -0.02345), C # (0.000453802, -0.0217284), C # (1.17588d-08, 1.22245d-07), C # (3.23091, -0.132288), C # (-1.28963d-07, -3.80257d-06), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-7.9517, 0.32562), C # (-0.000486679, 0.0367698), C # (-0.000746189, 0.0340719), C # (-0.0112059, -0.276819), C # (-0.00346891, -0.0875554), C # (-7.86952d-05, 0.0116247), C # (-0.000160771, 0.0107714), C # (2.91422d-09, 3.02963d-08), C # (-1.37313d-07, -4.04881d-06), C # (-14.6692, 0.600698), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (19.422, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.370409, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-6.45929, -0.264506), C # (0.000654894, 0.0494788), C # (0.0010041, 0.0458484), C # (0.0150791, -0.372499), C # (-0.012103, 0.30548), C # (-0.000274567, -0.0405584), C # (-0.00056093, -0.0375814), C # (1.01677d-08, -1.05704d-07), C # (-4.79086d-07, 1.41263d-05), C # (-7.9517, -0.32562), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.0408, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (6.3853, 0), C # (0, -0.0408), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000555262, -0.0419514), C # (1.98856, 0.0814425), C # (0.0170135, -0.42942), C # (0.000788511, 0.052829), C # (1.21157d-05, -0.00542571), C # (-5.57524, -0.228304), C # (-0.00165939, 0.0442559), C # (0.000264062, 0.0348546), C # (-0.000151981, -0.0201374), C # (-0.000486679, -0.0367698), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0.0408), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0.0408), C # (6.3853, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000851343, -0.0388733), C # (0.0170135, -0.42942), C # (1.98856, 0.0814422), C # (0.000385964, 0.0570138), C # (5.36068d-05, -0.00585705), C # (-0.00165937, 0.0442557), C # (-5.57524, -0.228304), C # (-0.00051019, -0.0322958), C # (-0.000294071, -0.0186603), C # (-0.000746189, -0.0340719), C # (0, 0), C # (0.0408, 0), C # (0, -0.0408), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (6.3853, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.0127851, 0.315829), C # (0.000788511, 0.052829), C # (0.000385964, 0.0570138), C # (1.98856, 0.0814079), C # (-5.57524, -0.228304), C # (1.29593d-05, -0.00545039), C # (5.443d-05, -0.00588152), C # (-1.02116d-08, 1.06159d-07), C # (0.0120437, -0.303168), C # (-0.0112059, 0.276819), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.0137295, 0.346532), C # (2.80669d-05, 0.00675756), C # (-2.35543d-05, 0.00729422), C # (5.10626, 0.2091), C # (-3.07027, -0.125729), C # (-0.000372858, -0.00503724), C # (-0.000334615, -0.00543482), C # (-1.15417d-08, 1.19988d-07), C # (-0.0142558, 0.353289), C # (-0.00346891, 0.0875554), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000311465, -0.0460089), C # (5.10626, 0.2091), C # (0.00210793, -0.0550618), C # (2.72943d-05, 0.00678016), C # (-0.000372858, -0.00503724), C # (-3.07027, -0.125726), C # (-0.00191428, 0.0407779), C # (-0.000498236, -0.0406443), C # (0.000288194, 0.02345), C # (-7.86952d-05, -0.0116247), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000636311, -0.0426318), C # (0.00210794, -0.055062), C # (5.10626, 0.2091), C # (-2.43082d-05, 0.00731662), C # (-0.000334615, -0.00543482), C # (-0.00191428, 0.0407779), C # (-3.07027, -0.125726), C # (0.000785054, 0.0376625), C # (0.000453802, 0.0217284), C # (-0.000160771, -0.0107714), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (1.15341d-08, -1.19909d-07), C # (-0.000216307, -0.0299928), C # (0.00042811, 0.0277909), C # (9.35258d-09, -9.72294d-08), C # (-1.15417d-08, 1.19988d-07), C # (-0.000498236, -0.0406443), C # (0.000785054, 0.0376625), C # (3.23091, 0.132321), C # (1.17588d-08, -1.22245d-07), C # (2.91422d-09, -3.02963d-08), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (-5.43468d-07, 1.60246d-05), C # (0.000124449, 0.0173293), C # (0.000246723, 0.0160582), C # (-0.0103514, 0.26089), C # (-0.0142558, 0.353289), C # (0.000288194, 0.02345), C # (0.000453802, 0.0217284), C # (1.17588d-08, -1.22245d-07), C # (3.23091, 0.132288), C # (-1.37313d-07, 4.04881d-06), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (19.422, 0), C # (-6.2442, -0.255698), C # (0.000492274, 0.0371924), C # (0.000754767, 0.0344636), C # (0.0113347, -0.280002), C # (-0.00325794, 0.0822305), C # (-7.39092d-05, -0.0109177), C # (-0.000150994, -0.0101163), C # (2.73699d-09, -2.84538d-08), C # (-1.28963d-07, 3.80257d-06), C # (-14.6692, -0.600698), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-6.45929, 0.264506), C # (-0.000555262, 0.0419514), C # (-0.000851343, 0.0388733), C # (-0.0127851, -0.315829), C # (-0.0137295, -0.346532), C # (-0.000311465, 0.0460089), C # (-0.000636311, 0.0426318), C # (1.15341d-08, 1.19909d-07), C # (-5.43468d-07, -1.60246d-05), C # (-6.2442, 0.255698), C # (-5.9819, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.1824, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0.000654894, -0.0494788), C # (1.98856, -0.0814425), C # (0.0170135, 0.42942), C # (0.000788511, -0.052829), C # (2.80669d-05, -0.00675756), C # (5.10626, -0.2091), C # (0.00210794, 0.055062), C # (-0.000216307, 0.0299928), C # (0.000124449, -0.0173293), C # (0.000492274, -0.0371924), C # (0, 0), C # (3.41863, 0), C # (0, -0.1824), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0.1824), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0.0010041, -0.0458484), C # (0.0170135, 0.42942), C # (1.98856, -0.0814422), C # (0.000385964, -0.0570138), C # (-2.35543d-05, -0.00729422), C # (0.00210793, 0.0550618), C # (5.10626, -0.2091), C # (0.00042811, -0.0277909), C # (0.000246723, -0.0160582), C # (0.000754767, -0.0344636), C # (0, 0), C # (0, 0.1824), C # (3.41863, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0.1824, 0), C # (0, -0.1824), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0.0150791, 0.372499), C # (0.000788511, -0.052829), C # (0.000385964, -0.0570138), C # (1.98856, -0.0814079), C # (5.10626, -0.2091), C # (2.72943d-05, -0.00678016), C # (-2.43082d-05, -0.00731662), C # (9.35258d-09, 9.72294d-08), C # (-0.0103514, -0.26089), C # (0.0113347, 0.280002), C # (0, 0), C # (0, 0), C # (0, 0), C # (3.41863, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.012103, -0.30548), C # (1.21157d-05, 0.00542571), C # (5.36068d-05, 0.00585705), C # (-5.57524, 0.228304), C # (-3.07027, 0.125729), C # (-0.000372858, 0.00503724), C # (-0.000334615, 0.00543482), C # (-1.15417d-08, -1.19988d-07), C # (-0.0142558, -0.353289), C # (-0.00325794, -0.0822305), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.000274567, 0.0405584), C # (-5.57524, 0.228304), C # (-0.00165937, -0.0442557), C # (1.29593d-05, 0.00545039), C # (-0.000372858, 0.00503724), C # (-3.07027, 0.125726), C # (-0.00191428, -0.0407779), C # (-0.000498236, 0.0406443), C # (0.000288194, -0.02345), C # (-7.39092d-05, 0.0109177), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-0.00056093, 0.0375814), C # (-0.00165939, -0.0442559), C # (-5.57524, 0.228304), C # (5.443d-05, 0.00588152), C # (-0.000334615, 0.00543482), C # (-0.00191428, -0.0407779), C # (-3.07027, 0.125726), C # (0.000785054, -0.0376625), C # (0.000453802, -0.0217284), C # (-0.000150994, 0.0101163), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (1.01677d-08, 1.05704d-07), C # (0.000264062, -0.0348546), C # (-0.00051019, 0.0322958), C # (-1.02116d-08, -1.06159d-07), C # (-1.15417d-08, -1.19988d-07), C # (-0.000498236, 0.0406443), C # (0.000785054, -0.0376625), C # (3.23091, -0.132321), C # (1.17588d-08, 1.22245d-07), C # (2.73699d-09, 2.84538d-08), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-4.79086d-07, -1.41263d-05), C # (-0.000151981, 0.0201374), C # (-0.000294071, 0.0186603), C # (0.0120437, 0.303168), C # (-0.0142558, -0.353289), C # (0.000288194, -0.02345), C # (0.000453802, -0.0217284), C # (1.17588d-08, 1.22245d-07), C # (3.23091, -0.132288), C # (-1.28963d-07, -3.80257d-06), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (13.1023, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (-7.9517, 0.32562), C # (-0.000486679, 0.0367698), C # (-0.000746189, 0.0340719), C # (-0.0112059, -0.276819), C # (-0.00346891, -0.0875554), C # (-7.86952d-05, 0.0116247), C # (-0.000160771, 0.0107714), C # (2.91422d-09, 3.02963d-08), C # (-1.37313d-07, -4.04881d-06), C # (-14.6692, 0.600698), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (0, 0), C # (19.422, 0) / open(1, file='bb', status='old') read(1,*) ((aa(i, j), i=1, 40), j=1, 40) close(1) call zheevd('V', 'U', 40, aa, 40, w, work, lwork, rwork, lrwork, $ iwork, liwork, info) write(*,*) 'info = ', info end