Исходные тексты программы

Исходные тексты программы. REAL A3,3, U3,3, V3,3, SIGMA3, WORK3,Y3,C3,Y03 INTEGER I,IERR, J, M, N, NM OPEN 6,FILESVD.OUT,STATUSUNKNOWN,FORMFORMATTE D OPEN 5,FILE SVD.IN,STATUSUNKNOWN,FORMFORMATTED 140 FORMAT3I5 150 FORMAT4E15.7 READ5,140 NM,M,N DO 131 I1,M READ5,150 AI,J,J1,N 131 CONTINUE READ 5,150 YI,I1,M CALL SVDNM,M,N,A,SIGMA TRUE U TRUE V,IERR,WORK IFIERR.NE.0 WRITE 6,2 IERR 2 FORMAT15H TROUBLE.IERR,I4 WRITE6,120 120 FORMATМАТРИЦА А DO 121 I1,M WRITE6,130 AI,J,J1,N 130 FORMAT8E15.7 121 CONTINUE WRITE 6,160 YI,I1,N 160 FORMATПРАВЫЕ ЧАСТИ8E15.7 210 FORMATСИНГУЛЯРНЫЕ ЧИСЛА WRITE6,210 DO 3 J1,N WRITE6,6 SIGMAJ 3 CONTINUE SMASIGMA1 SMISIGMA1 DO 211 J2,N IFSIGMAJ.GT.SMA SMASIGMAJ IFSIGMAJ.LT.SMI.AND.SIGMAJ.GT.0. SMISIGMAJ 211 CONTINUE OBUSMASMI 230 FORMATЧИСЛО ОБУСЛОВЛЕННОСТИ,E15.7 WRITE6,230 OBU SIGMA10. DO 30 J1,N IFSIGMAJ .GT. SIGMA1 SIGMA1SIGMAJ CJ0. 30 CONTINUE TAUSIGMA10.1E-6 DO 60 J1,N IFSIGMAJ.LE.TAU GO TO 60 S0. DO 40 I1,N SSUI,JYI 40 CONTINUE SSSIGMAJ DO 50 I1,N CICI SVI,J 50 CONTINUE 60 CONTINUE write 6,560 WRITE 6,6 CI,I1,3 DO 322 J1,N SS0. DO 321 I1,M 321 SSAJ,ICISS 322 Y0JSS write 6,570 WRITE 6,6 Y0I,I1,3 C WRITE6,7 C DO 4 I1,M C WRITE6,6 UI,J,J1,N C4 CONTINUE C WRITE6,7 C DO 5 I1,N C WRITE6,6 VI,J,J1,N C5 CONTINUE 6 FORMAT3E15.7 560 format2x, roots 570 format2x, right 7 FORMAT1H STOP E N D SUBROUTINE SVDNM,M,N,A,W,MATU,U, MATV,V,IERR,RV1 REAL ANM,N,WN,UNM,N,VNM,N,RV1N LOGICAL MATU,MATV IERR0 DO 100 I1,M DO 100 J1,N UI,JAI,J 100 CONTINUE G0.0 SCALE0.0 ANORM0.0 DO 300 I1,N LI1 RV1ISCALEG G0.0 S0.0 SCALE0.0 IFI.GT.M GO TO 210 DO 120 KI,M 120 SCALESCALEABSUK,I IFSCALE.EQ.0.0 GO TO 210 DO 130 KI, M UK,IUK,ISCALE SSUK,I2 130 CONTINUE FUI,I G-SIGNSQRTS,F HFG-S UI,IF-G IFI.EQ.N GO TO 190 DO 150 JL,N S0.0 DO 140 KI,M 140 SSUK,IUK,J FSH DO 150 KI,M UK,JUK,JFUK,I 150 CONTINUE 190 DO 200 KI,M 200 UK,ISCALEUK,I 210 WISCALEG G0.0 S0.0 SCALE0.0 IFI.GT.M.OR.I.EQ.N GO TO 290 DO 220 KL,N 220 SCALESCALEABSUI,K IFSCALE.EQ.0.0 GO TO 290 DO 230 KL,N UI,KUI,KSCALE SSUI,K2 230 CONTINUE FUI,L G-SIGNSQRTS,F HFG-S UI,LF-G DO 240 KL,N 240 RV1KUI,KH IFI.EQ.M GO TO 270 DO 260 JL,M S0.0 DO 250 KL,N 250 SSUJ,KUI,K DO 260 KL,N UJ,KUJ,KSRV1K 260 CONTINUE 270 DO 280 KL,N 280 UI,KSCALEUI,K 290 ANORMAMAX1ANORM,ABSWIABSRV1I 300 CONTINUE IF.NOT.MATV GO TO 410 DO 400 II1,N IN1-II IFI.EQ.N GO TO 390 IFG.EQ.0.0 GO TO 360 DO 320 JL,N 320 VJ,IUI,JUI,LG DO 350 JL,N S0.0 DO 340 KL,N 340 SSUI,KVK,J DO 350 KL,N VK,JVK,JSVK,I 350 CONTINUE 360 DO 380 JL,N VI,J0.0 VJ,I0.0 380 CONTINUE 390 VI,I1.0 GRV1I LI 400 CONTINUE 410 IF.NOT.MATU GO TO 510 MNN IFM.LT.N MNM DO 500 II1,MN IMN1-II LI1 GWI IFI.EQ.N GO TO 430 DO 420 JL,N 420 UI,J0.0 430 IFG.EQ.0.0 GO TO 475 IFI.EQ.MN GO TO 460 DO 450 JL,N S0.0 DO 440 KL,M 440 SSUK, IUK,J FSUI,IG DO 450 KI,M UK,JUK,JFUK,I 450 CONTINUE 460 DO 470 JI,M 470 UJ,IUJ,IG GO TO 490 475 DO 480 JI,M 480 UJ,I0.0 490 UI,IUI,I1.0 500 CONTINUE 510 DO 700 KK1,N K1N-KK KK11 ITS0 520 DO 530 LL1,K L1K-LL LL11 IFABSRV1LANORM.EQ.ANORM GO TO 565 IFABSWL1ANORM.EQ.ANORM GO TO 540 530 CONTINUE 540 C0.0 S1.0 DO 560 IL,K FSRV1I RV1ICRV1I IFABSFANORM.EQ.ANORM GO TO 565 GWI HSQRTFFGG WIH CGH S-FH IF.NOT.MATU GO TO 560 DO 550 J1,M YUJ,L1 ZUJ,I UJ,L1YCZS UJ,I-YSZC 550 CONTINUE 560 CONTINUE 565 ZWK IFL.EQ.K GO TO 650 IFITS.EQ.30 GO TO 1000 ITSITS1 XWL YWK1 GRV1K1 HRV1K FY-ZYZG-HGH2.0HY GSQRTFF1.0 FX-ZXZHYFSIGNG,F-HX C1.0 S1.0 DO 600 I1L,K1 II11 GRV1I YWI HSG GCG ZSQRTFFHH RV1I1Z CFZ SHZ FXCGS G-XSGC HYS YYC IF.NOT.MATV GO TO 575 DO 570 J1,N XVJ,I1 ZVJ,I VJ,I1XCZS VJ,I-XSZC 570 CONTINUE 575 ZSQRTFFHH WI1Z IFZ.EQ.0.0 GO TO 580 CFZ SHZ 580 FCGSY X-SGCY IF.NOT.MATU GO TO 600 DO 590 J1,M YUJ,I1 ZUJ,I UJ,I1YCZS UJ,I-YSZC 590 CONTINUE 600 CONTINUE RV1L0.0 RV1KF WKX GO TO 520 650 IFZ.GE.0.0 GO TO 700 WK-Z IF.NOT.MATV GO TO 700 DO 690 J1,N 690 VJ,K-VJ,K 700 CONTINUE GO TO 1001 1000 IERRK 1001 RETURN E N D ПРИЛОЖЕНИЕ 2. контрольный пример Входные данные матрица изначально сингулярна первая строка равна сумме второй и третьей с обратным знаком 3 3 3 .320E 02 .140E 02 .740E 02 -0.240E 02 -0.10E 02 -0.570E 02 -0.80E 01 -0.40E 01 -0.170E 02 -0.140E 02 0.130E 02 0.10E 01 Полученный результат МАТРИЦА А .320E02 .140E02 .740E02 240E02 10E02 570E02 80E01 40E01 170E02 ПРАВЫЕ ЧАСТИ 140E02 .130E02 .10E01 СИНГУЛЯРНЫЕ ЧИСЛА .1048255E03 .7310871E-06 .1271749E01 ЧИСЛО ОБУСЛОВЛЕННОСТИ .1433830E09 Корни .1215394E01 .1821742E01 1059419E01 Правые корни после проверки 140E02 .130E02 .101E01 Видно, что правые части соответствуют начальным данным.

Решение верно.