READ m.IFID

PRINT 92e,C(Ae(I,J),J=l,N),l=l,N)

PRINT 921,CCBCI,J),J=1,IR) ,I=1..N) PRINT 922,C(FieCI,J),J=l,N),I=1,N) PRINT 923,UPSI(I,J), J=1,IR) ,I=1,IR) Vim FORMAT (2Х,Ав=,йЕ12.5)

921 FORMAT (2X,BaME12.5)

922 FORMAT C2X,Fie=,i!.E12.5)

923 FQFSiAT (2X,FSI=ME12.5) PRINT 905,(FIVAK(I),I=1,5) PRINT Vie,IFID

CALL ARI-iAYC2,N,N,Nl,Nl,Aevi,Ae)

CALL ARRAYC2,N,IR,N1,IR1,BV,B)

CALL MTRACAeV,AeTV,N,N,e)

CALL Mn«*(BV,BTV,N,IR,e)

CMX ARl.AYCl,N,N,Nl,Nl,A»Ty,fteT)

CALL ARRAY(1,IR,N,IRl,N1,ВТV,ВТ)

CALL MCPY(P3I,P3Il,IRl,IRl,e)

CMX ARRAY С 2, IR, IR, IRl, IRl. PSIIV, PSIl)

CALL MIhW(PSIlV,IR,D,L;M)

CALL ARFWYС1,IR,IR,IRl,IRl,PSIlV.P3I1)

PRINT 666,Ш,BV,AdTV,BTV,PSI1V 666 F0RMATC2X,H>",6E12.5) С ВВОД НАЧАЛЫ*1Х УСЛОВИЙ ДЛЯ РЕШЕНИЯ СИСТЕМЫ А,Ч-ЕБРАИЧЕСКИХ УРАВНЕНИЙ

FiEAD 9ei,C)CS(I),l=l,NXS)

REAB 917,EFS

READ 9в6,IT 9в7 FQiaiAT (2X,X HA4.=,6E12.5) 908 FORMAT C2X,IT=M5)

PRINT 9e7,CXS(I),I=l,NX3)

916 F0RMAT(2X,EF3=,E12.5) PRINT 916,EPS

PRINT 9вЗ,1Т READ 917,P0PPR1 PRINT 913,Р0РРТ:1 913 FORMATC2X,nonp.B PRl ДЛЯ РЕШ.АЛГ.УР.,C12.5) READ 917,CPRMT(I),I=1,4) READ 917,HTIME

917 FOIdlAT (E12.5)

909 FORMAT (2X,PFSir=,4E12.5,HTIME=,E12.5)

PRINT 909, С FRMT CI),I=1,4),HTI«E 911 FORMAT (9F6.3)

READ 911,CDERY(I),I=1,NDIMRK)

flETURN

С ВНЕШНЯЯ ГГОДПРСГРАММА, ВЫЧИСЛЯЮЩАЯ ЗНАЧЕНИЯ ФУНКЦИИ Ш РЕШЕНИЯ СИСТЕМ С НЕЛИ1-ЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ С ВИД УР-НИЛ FI(X1,X2, ...XN)=e, 1=1,2 N

SUBFOUTINE FUN(C. ХЗ)

COMMON /А1/ Ав, и. Fie, F3I

COMMON /А2/ АвТ, ВТ, PSIl, FI

COMMON /А4/ N,ir-:,NXS,NDIMiy<.,Nl,IRl,NX.Sl,№iIMr<l , COMMON /А5/ AeV,AeTV,E»,BTV,FSIlV FE/iL Ав(6,6),В(6,2),Р1в(6,6),Р31(2,2) FiEAL AdTC6,6),BT(2,6),PSIlC2,2),FIC6,6) REAL A(6,6),DC2,6) FEAL XS(21),C(21) F<EAL K(6,6),CM(6,6),KV(6,6) FdlAL РР:1(6,6),Р1<2С6,6),РТ!21(6,2),РР22(6,2), * PF<23(6,6),PR24{6,6)

REAL 31(6,6),32(6,6) С ВЕКТОРЫ ДЛЯ 1ГРЕДСТАБЛЕНИЯ МАТРИЦ ВЕКТОРНЫМ СПОСОВОМ-REAL AeV(36),AeTV.36),BV(12) ,BTV(12),P3I1V(4) F£AL KKv436),FIV(36),CMV(36)

REAL FT<1V(36) ,РР2У{36) ,РГч21У(12) ,PR22V(12) ,РГ<23У(36) ,F«24V(36), S1V(36),32У(36)

Ю 1ввЗ 11-1, Ml DO leey J1=1,N1 КУ(11, JD-e.

loets C0(4TiMUt:

iiU 101© I1-1,H DQ 1в2в Jl-1,«

CALl. L0C(Il.Jl,iy,(4,U,l) KV(I1,J1)=X3(IV) leze США IrtUE 1в1й CONTINUE

UmU. rtCPY(KV,K,Nl,Nl,e) С ПОЛУЧЕНИЕ МТРИиЫ CM(N1,M1) ,ГДЕ ХГАНЯГСЯ ИТиГОМ-Е ЗНАЧЕНИИ С ЛЕВЫХ HACIBl УШЖ1Ш, ВСЕ ДЕЙСТйИ)! НАД МАТГШАЧИ. 111-ТДСТАВ.ИННЫЧИ

С В ш-;то1. ФОРМЕ

CALL т<АЧ(.?,ti,n,i4,m,\асу,i;>

CALL AKR.WiZ.d.H.m.MbnV.n)

CALL 0НРГЧ11 ((ЗчУ, АвУ Д 1<1У, N, N, И)

CALL шт(АвТУ.IачУ,гГСУ,Н,а,Ы1

CliLL OilI-RD 11\1\У, ВУ, F -Id У, N, N ,1R)

JALL 0П1гаь;пЧ21У,Р311У,ГРа2У,М,1КДР;)

CALL WRD (1ТС2У, БТ У, РГСЗ., N Д R, N)

CAUl OrtlRLi(r-f:23y, 1ч1\У,(ТСчУ,,N,N)

CALL GiWiIi!PRlУ,IFiJl.31 У,П,П)

CALL MSUh. 31У, ГКгУ, 32У, N. N)

С-i-L G1IADD (52У. F1У, CflV, N, t)) С QEPATilOE ПГЕС)БГА30а-М)ЯЕ К «АГРИЧНОМ> И-ЩУ

CiiLL А1л1 iAY (1, tl, Н. N1, N1, ему, Crt) С IIOAIDUt ИЗ liniTMlW СП ГСКТОРА CiXS) ,С00ТГсЕТС1П/«Щ1Г0 F£ С BETXlCt ТРОГОЛЬКЧ!*.;

DO U~1,N

UhLL LOCCll, ЛДУ.М.МД ) СДУХ.гиИ, JI)

10-iO coNii«:

1330 l.0f4IINUL RLTURN UID

L ПОДНГОГРАММА ВЫЧ1СЛЕЛИИ ПРАВЫХ ЧАСШН CHCTUibI ДЦ(ЛфЕГи-Щ11А.1КН1К

С УГАВННШ

GUBROUritlE ГСГ(Х, ГДСРГ)

[iCAL Y.:36),t€Rf,3i) С МмСИСЫ ДЛЯ ВВОДА НАЧА.ЛЬНЫЧ

REAL Ав(й,й),В.6,2),Г1в,й,й) ,ГЗ[.2,2)

REy. АвТ(6,6),В1(2,й),Р311(2,2),Р1(6,6)

RLAL К<6,6) С iVnCCHIiU ДЛЯ 0rr.M-M3(W« СЧЕ1А ПРАВЫХ ЧАСТЕЙ

rCAL В.2,6),А(6,6),АТ(6,6), , « L(u,u),

-< А1\(3й,3й) .АККЗй.Зй) ,DMINy,L«INyC3u) ,МИ1МУ(Зй)

REAL Г1РК(ЗйЛ) ,riP(i,6) ,РИ\Г<3(Зй, 1) С й11Ди-;С ИСС.1ЕДУЕМ01-0 3JB1QHTA MATPtmU FI

INTEGER IFID

I4EAL АвУ(Зо),АвТУ(Зй) ,ВУ(12) ДЛУДг) ,Г311У{4)

1-£AL АУСЗЛ),ЬУа2)

REAL АТУ(ЗД),AhCyt 1296) .AKiyC 1296)

I4LAL рга*1У(129й) ,PRKFCVC129U)

COMMON /А1/ Ав,Е,Г1в,Р31

CQMi-iON /А2/ АвГД-:Т,Р311,Г1

C0Mt\3N /iriD/ IFID

COMMON /А4/ МДР,№<3,М111Мга\,М1Дк1,№(31,М111МК1 COMMON /A3/ АвУ,АвТУ,ВУ,ВТУ,Г311У С nCPECTPUW<A ВЕ1Л0-А ПЕ11-Ж1НЫХ Y.NIIIMFUO В МАССИВ K(N,N) DO 1210 IY1=1,N

DO 1220 IY2=1,N

IY3=lr2t(IYl Dwri l.(IYl,IY2)==r(IY3) 122« CQNriNUL 12ie COMIINUt;

CiUi. iVWDD (A,e, » fSIl,BT,K,

* Ae,B,M,IR,Nl,IRl,AV,DV)

DQ 1230 IC1-1,N1 DO 1230 IE2=1,N1 Eaa,IE2)=e.

IF (IE1.EQ.1E2) EaEl,lE2)=l. 1230 CONTINUE N2=N»»2 N12=N1»»2 N22=N2»»2 N122-N12»»2

CALL mRA(ftV,ATV,N,N,0)

CALL ARF:AY(1,N,N,N1,N1,ATI.,AT) € ВЫЧИСЛЕ)-ЬС HATPVWI AK(N2,N2> COCTOillfW rt3 СУММЫ №OHEEPOBUX С ПГ0ИЗВЕДЕ1М1

CALL aON---K(AT, E, N, N2, N22, N1, N12, N122, FRKRlV)

CALL □iQNPR(t,AT,N,(C,N22,Nl,N12,N122,PRl<R2V)

CALL GMADD (PFMlU.fWCy,Ab;V,N2C) С ПОЛУЧЕНИЕ М,ОТШ1 ПРОИЗГОДНЫК D FI/D FKI.J) CN»»2,1) С IFID - ИНДЕКС ЗЛЕ«01ГА МАТРИЦЫ Р1,В.ПИЯНИЕ KOTOf-OIO ИСС:ЛЕДУЕТСЯ С П0ЛУЧЕИ41Е МАТРИЦЫ 1[Р0ИЗВ0ДНЫХ riP(N,N)

DO 1010 1=1,N

DQ 1010 J=1,N FIP(I,J)=e. 1010 CONTINUE

DO 1030 1=1,0

1030 iFd.m.iriD) ripa,ii-i.

с получше МсТРИЦЫ FIPK(N)<«2,1) ИЗ FIP(N,N) DO 1040 11,N DQ lOiO J-1,N L=-Jia 1)»N FIPK(L, l)=rH4I, J) 1050 CONTINUE 1040 CONTINUE

С OE;i-r=illEIML МАТРИЦЫ Al\IV(42•N2) .ЗАДАННОЙ В РИДЕ ВНчТОРА ВО 1060 1=1,N2Z 1060 AKiyCD-iWd) 1112 FORM. I (1X, AKIV, 10Е11.4)

CALL MINV,AKiy,N2,DHINV,LMINU,MMINV) URITE (6,1112) DMINvi,AKlV С ПЕРЕМНОШ-МЕ МАТРИЦ AK1(N2,N2)»FIPK(N2, l)==F«sJ<3(N2,l)

CALL GMI-RD (А1Ч1V,FTPK, PFja"<3,N2, G, 1) С ВЫЧИСЛЕНИЕ 1ПАИЫХ ЧАСТЕЙ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНВМН С DERY=D (1\К) /В С FI (I, J) = - PRKR3

НО 1070 1=1, h42 1070 BERY(I)= РГМа{1,1) RETUFi4

С подпр-ал-АммА ПЕЧАТИ оигр

ЗиВГЧОи ГINE ОиТР(X, Y, UERY, IHLF, NDIM, РКМТ) F£AL Y(36),IC1Y(36),PIWC5) F£AL К(6,6)

REAL DC2,6),А(6,6),AV(36),DVC12)

FCAL АвТ(6,6)Д!Т(2,6),Р311С2,2),Р1(6,6)

F£hL Ав(6,6),B!6,2),FieC6,6),PSIC2.2)

FiEAL EVr-t(6) ,EVI.6) ,VECK(6,6) ,i,tCI(6,6> ,XARF0L(7) ,UC6) ,У(6)

rCAL FACT(6) ,SUBBIA(6) ,mRK(6) ,UQRK1(6) ,U0I4<2(6)

INFEGER II.JQ1"K(6),LDCAL(6)

INTEGD;- lER(A)

COMMON /А1/ AO,B,riO,PSl