RETURN

ERD


формам И


f.

J-

t I fi £

i] $ J

і












































E9-9O* grPOJ




(pv’wy‘sfA‘оу'ях* ax high iiyo (Г)СЛ~( j'mA Z2

m‘ї=с гг ос wh’Mio тг

ЛМ**( 1 )Х=( м ‘ 1 MX

ЗУ* 1-1 U ос І-У-АУ тг ог 'І-.'І'їШ 02

оМі-зі йг се

ЗИЯУЙПІ.ИЧІОф =s== З

( ( ПГА-( 0'1 )ГА )5<П'П~ГМ ( ) * j

( (ОК )М МчУС-ІМИ 0* Т * ЬЯ * WK ) 4 1 91 імшоо и

‘м< m муч; г? H’- + (zi>a:<+< t і mi-c с )са

М = { Х + ЯІ )дч ?т

М 0-109 ( сто ‘91 > 41 1-SI-91 (91 > Д/А-Г T + SI )Л’л’

01 0109 ( (91 W3!OW 41 Я

( qf );л=^ Т-СС-9І

01 ‘ 2--9Г 9 І ПО

( I ‘ f ).М( 1 Ш1 21 от* мі гт оп ;V’‘ W п 00

A

:=ss=sss КЯ9УЛ09

Yl^'UHO ОООгЧІОУ^Я V'K.HShVHP УМЯЭ<ІМІА0-0 =х== О

ч-сп^сі

1+70=^ї

1+М^Ії

>/Ш = Н 9

91 ГНОЭ

ЯСУ 1ІМоО 9

( 1'Г)Х=(С)СЛ ОУ*І=Г я ОС

9 0X09 ( 1 ‘ЯЧ‘01 ) .41 ?/C1=7l

. HW*W’ЧОІІ/Р.О/'ІОлЧОО

<4А' 'bl * nV * ая * М’/КХУЧ/^Оа'І'ІОО

117' X7W* ОЯ ‘ U!4 ' нч ' ( 01) ДУ 9*1У4Ч

( Т’НлКУ* г І'ЧІ’ЗІАУ*

’ (>» X * Н V*) ОУ* ( 01 *НЧ )9Х' ( 1 * ОУ )я У* ( ЧУ ‘ 0:9 МХУ

* ( I ‘ ON )SC X * ( Ov 1 ХУ9 ‘ ( 0?4 )Г А * ( 01 ’.ЗУ ) А‘ ( О’1’ )Х ЧИС’МЧТС

Ч-УЯ0Я1'1

( 7-0 ' Ч-V )9*1У.ЧУ LlOlldyl

( 75 * ISO ’ 074w?0osq ’ ^У * .MV-*

‘рХ ' ЯХ '5У’ЯХ' XVS'АЧ’йГА ’ СЛ ‘ A * Х)9¥У 1’1 TlfWOlS

****«*««*«** ИХИАОИаііНУїіУХ и.ОКЬОЯОйЙ.АЧУМЛ VJ -h!?¥4 УЧЧУеІJOd’JeOl У*#* з

frS’If W-643W6 10

0


J




OCT 92-4279-80 Л,25



0 иХ£*Х

:= РАСЧЕТ ОСНОВНОЙ ПОГРЕШНОСТИ ===s==s==x=xsss^^3=5;k==x=s======3s:




RYAX=0t ■ ' -

DO 32 J*1,RC

Р5АХ=Й„

pft=xcJ)




00 30 K=!,MR



30

PSA№PSAX+A3(K,1 )*XF( J(K) SAXU)=PSAX




DO 31 1 = 1 ДС



31

RYAX^KYAX-HYU, I )-PSAX)**2



32

CONTINUE

Zf-i-l LC*NC-M**1 *HM1

П (2H,GT,0,1 GOTO 33 £KR = l GOTO 301



33

GAMMAC=100t*KD*DSQHT( RYAX/ZK*SSG)



C =-=

x ЕСЛИ КРп-0 .РАССЧИТАТЬ СЛУЧАЙНУЮ СОСТАВЛЯВШУЮ GSl==™=™=~=



C — SS

x И ПОГРЕШНОСТЬ АППРОКСИМАЦИЯ -===s=====^=:=-==-.=-==^z=z--




If (NP,EQ,0) GOTO 301

HYV=0»

DO 41 J=1,NC



40

P*=Y J( j ) .••-.•

DO 40 1-1,LC

RYYJ=aYYJ+( ( YU, I }~Ptf)**2) ■



41

continue

§



ZA-WMl#ftMl#NC*(ЬС-1)

IE (ZN.GT,0,) GOTO 43



f

£ER=2 GOTO 301'

KJ


43

GStxiJSij^KDfcDSQHTtRYYJ/ZN+SSGJ

*3


.. C ==z

= . РАСЧЕТ ПОГРЕШНОСТИ АППРОКСИМАЦИИ ИАИ СРЕДНЕКВАДР,ЗНАЧЕНИЯ ==»

TO Cl



IE (КОМ,EQ,1) GOTO 55




MAXMJAoS(YJ( 1 )-SAX( 1 ) )




DU 50 3=2,NC

P«=DABS1YJU)-SAXU) )


50

IE (PE,GT,MAX) HAXSPW

GOtiTUb'E


A



GA=KAX/AM1*100,




GOTO 301 . ' ’

"4}

C* .

V

05

UAX=0.

Lu t> *ї і ‘іC

£*>» tj

66

HYXa = £'i*X-HY JU)-SAXU) )**2 .

ZN =KИ1 * E .41 * ( AC — M—l 1

<=Q



If (2N»GT»D,I GOTO 5B

Cj

o,.


ERrt=3

&


GOTO 301

<□


.56

GA=DSQRT(RYAX/2N )

'TQ .D

c

4

301

IP (ERRFEQ,0 ) ■ RETURN

IF {EEK;£q»lJ PRINT 60

IF (ERR,£QB2) PRINT 81


'•' ІЇ

IF (EFA,£Q,3) PRINT 62



60

FOHMATt'. В РАСЧЕТЕ GAMMAO ЗНАМЕНАТЕЛЬ »0Г)

2s

&

CS

61

82

FORMAT!r В РАСЧЕТЕ СЛУЧАЙНОЙ СОСТ.ЗНАМЕНЛТЕЛЬкЭ*)

FORMAT!’’ В РАСЧЕТЕ GAMMA АППРОКСИМАЦИИ ЗН АМЕН АТЕЛ Ь=0 ” )

T

r


EKD

*3>

o


• .... •:

ЇЗО

41











ГОСТ $. (05-68 фовмал H



-

І OCT 92-4279-60 Л.26

4 t

•z •■

Q < ....'■ J

G &4t# ПОДПРОІ PAM/A РАСЧЕТА КОЭФФИЦИЕНТОВ МАТРИЦЫ **************л###;

SUBROUTINE KoEr(X I,XU,XG,YJtAPT,VOUT)

,; DOUBLE i-'Et.3Ir. 1 j:; XI , XB ,X0 , У I , APT , VOUT

DOUBLE ЫіЕС IS І UN- KMKb s ■'

•COmmUA ZRAZM/NM , KD, № ,LC , NP

CCMMCL/US/L.0,4 t e % H

( X

и*

DIMENSION X И <C ,MR) , XB( <*R( N0 ) t.

*XU(MRt;<h), - ' ’ ' .

. *YHKC ,1 ) ,АРТОіІЇ,1 )', VOUmK, 1 ) '

Lh-1

.. GALL THANSrUI ,X3 ,.<G ,РЯ) ' . .

CALL KULTIIXS.XI ,XC,HKДіС-.МЙ)

CALL IRVERSIXC,MR)

CALL MULT I ( Xb , YI ,,A PT , P R, fiC t LR) ’ .

CALL MULT I ( XC , APT r VOUT, k-kj rifc^LR ) .

уа< Псд/шев Зст?О| j „ 1

RETURN '

■ lad . • »

C GO ЦЙ PG Г P A ММ А ТРАНСПОНИРОВАНИЯ МАТРИЦЬ! f****4>*^*****^****^$<t>|c*i

SUBROUTINE THAXSPC A , В , JiR, MR) . ' ■ ■

DOUBLE PRECISION AtB

D IM L H SIG к A ( X ft (i-l R ) f В ( M R, N R )

DO 11 I~1,«ri ' ,

DU. •! J-l ,y;h

B<J,I)~A(I,J)

11 СО.ХТІЛОЕ


RETURN &

■ЧЭІ г

<5

BED

SUBROUTIXL MULT I(A,b,C,LR,«R,RR) C ПОДПРОГРАММА > ХИОлс rl И Я МАТРИЦ

4Q

*4 » д йк 0 <й~ гО «)

О

В- 3

7

DOUBLE PRECISION A,BtC,SUM

DIMENSION A’( A i< , M R ), В ( MR, К R } (0 ( A R t К R )

DO И L'-l ,ЛД .-

DU 1 1 -'J-l , KH

SU/.=B,Z ' '.

DO 13 K=1tMH

13 SUM=SUM+A(t(K

C(1,0)=SUM

IICORTIXUe l.

ReTUHL . I ‘

< END л " ■

•s'-.' f

. •*

о

*

Q

Р «X тГ і

r

о



гс

(

ГПСТOCT 92-4279-30 Л.27

uO Іін PO Г PA У<У: А ИНВЕРТИРОВАНИЯ МАТРИЦ^ *«x« ««♦***i<i»*+*^

SUBROUTINE I AVE RS(A tMR)

uUUBEE PRECISION A,В,C,DET,Y.«,EPS IL

INTEGER P,RfZ

DIMENSION A(MR(NR),Zl12)fB(12),0112)

DET=1,0

00 11 j-l,fiR

ZIJ)=J

DO 12 1=1,NR

л=І

Y-AI ItI }

R=I-1

IE (p,G'IKR) GO To 13

90 13 J=P,NR

-¥-A(I,J)

IE (DAoS(W),GT.DaBS(Y)) GO TO 32

00 TO 1032





DtT=Y*DbT

EHSIL-г,10-26

IE (DАВЕ(Y), 1Д ,EPS IL) GO TO 3

3

х

■3


4 4>

J

0 ro


со'


зз


16


15


17


1 »

І 5

12


FORMAT tlh , 'ЭЛЕМЕНТ МАТРИЦЫ^' ) PM KT It return

Y=1,/Y ‘

DO 17 0=1, Mt

DCJ)=A(J,K)

AlJ,К >-A(J, I )

Al 0,I)-~C(0)*Y

Al I,J)~A(I,J)*Y

BW J=A( 1 , J)

с о ( T 11 0 £■


J=Z(I)

Z<J)=Z(K)

Z(K)=J

DO 16 E = 1 , f. К

IF (K,Ek,I) GO TO 16

D <j 16 v — 1 , И A


А ( К,J)-л(К,J)-&(J )#С(К)






















OCT 92-4279-80 Л,29


1

2

3

4

S

6

7

a

21

22


*****’*r**********JK**#*****>x* + »**xt***********^4c^*ifc#**^;****^4sj([*j|C4

РАСЧеіГ К-ТОВ ФУНКЦИИ ВЯИйНгіЯ и ДОП f ПОГРЕІїНООТИ .

|*'Ж**1і'***'’^ж*чсН^**ї*сі*.*->и3іКХ«******-****>к5фгчс*ч«^-^И4**і*і*ЧЧг*-#*>к**4«4с**4*4»>К*^*^

IMPLICIT K£At*6(A«H,U-Z)

D l.Mt AS J OX KS1 (6 ),YC б , 4 ) , X І б ) , Yh ( 2,4),ТЄТА ( 6,6 ), В ( б ), FVU ) MAL*fe KSI , KSI ii, КЗ II , K51 IU , il, HI

CMMf./U/HI ,MXFV, JT, JPR, IUt IV.ZIV

I-TtGbR ZIV,MU(6 ),MV(6)

DATA Y Ю/З ,3,3,4.6,6/

DATA MIV/1,2,3,3,4,5/

FORMAT!' ",/////)


fb

£f rj J

CO

<>


rGrtXATI'

FOiUuiTHIl ) РОЙУ.АТІ r ',411)

F0 Rv.AT (6F6 , Г) FORMAT!' '.6F8.3) FJRMAT(3A6)

FURMATi' ',3Ab ) FORMAT(11,F6 * 1) FORMAT!' 'tIl»F6t3) pRIhT 1 read г

PM XT 2

READ 2

PMXT 2

MAD 3, NFV, JT, JPRJIP


at сил Э ,*

S1’E>* O-es- <n

t-Q

Q ’


31


PHIAT 4(XFV,JT,JPR.NP

l'r UF\bT,l,0R.£iFV,GT,6) GOTO 5Ь6 . lV=y.IV(NFV}

IU=4IU(AFV)

READ 7,rV-

PKINT b,FV

MM 5,KSl

HW 6,(KSI(I),I=1,IU) дЬП sKSI ( 1 ) KbI IU-КЫ ( ГО ) J? ( 0FV,LE3) GOTO 30 HtAO b,(X(I ),1 = 1,IU > hint a, (x (і), i=i, IU > DO 31 01=1,JT

MAD 5, ( Y( I , JI J ( 1=1 , IU)

PAINT 5,(Y(1,Л ),1=1,IU)


О

32

MAD 5, ( Ш I , M , 1=1,2)

PKIRT o((YU(1 , J ) , 1 = 1,2)

<s>

1


fc=DAbS(YNt2,1)-YN(1,1)1'

GOTO 24

o



м3 £3r>



ад


[r (JPR,X£,1) GOTO 33

DU 32 J=1,JT















































21


22


31

23


DV 23 1=1 tiiH K=Z( I )

П <К,еЛ» I ) GO ТО 31

DO 22 0-1 , ftй iV=A( I tJ) Al ItJ)-A(KtJI A ( К, J ) -Vi COjiTLIJUE P=Z(I) г (і)-z(к) ZU)=r ' OtT=-DET GO TO 21 GOXTIXuE COATI fiUE RETURN END


OCT 92-4279-80 Л,28


§ OJ

3 о tj П с:


-“М


G rO C2

Q


*5»


ftoi
















































33

V

23 '

23

24

34

С- s=x

ОСТ 92-4279-60 Л. 30

•f

READ b.f ( Yi( I ,1 ) , 1-1,2) ' • ■ •*■

PKIUT. ;> , ( Yi.U , 1 ) . i-1 ,2 )

КЧ)*86(П(г,1)-и(1*1))‘ - ■

IE(nfv»gts3) goto йЬ

A~jr.3so::t i (i)) .2. '

ЇМ 2,1 )~YM1 ,1 ) f /

DO 23 i-2,3T • ♦ • •

YM 1 > 1 }-¥&<! ,1 )

У ii ( 2 , I )~YH2,1 )

MMIMt L - ■ Л

REaD 5,hSlii '

PMM b,KM:. ••.<••

IE (KSIU.EQ.0,) K5IU»20, t

READ 21 ,ZIV,HI

PRINT 22.ZJV.riI ■ '

DO 34 1-І,ID . '

KM.! 1 і I )-KSIN

CALL FVGSTC Ka£ tX ,¥,Y;i ,TETAtGiiD,GG )GO,G)

sstsss D&iBO Ji РЕЗ.У llbT ATOB ssss - = “s“ = ::;r=:==2'xzs“ss*=x“--




400

500

417

403

404

405

406 •

407

406

409

333

■ 41

51

61 ■

71

61

DATA В/'Ц0','Nil ' ,'м2','CO','Q4'/ -

FORJ-AT(' f . 2 OX . f ФУ ri KU HR В Jr И Я ri И a * ) ■

FORMATC './Ж )) ■ - ' ,

FORMAT ( ' r ,4X , ' BjU'GUAa В ЕЛ И 40 A f ,2 іX , 3 A Q ), ji*

FOR’iATC* ' ,4X , 'Л2А0АЗОН И 3 4 E К О H BOOOO)

F0RMAT(Z ' ,4X ЛвЕЛйЧПОі’’ , 31X.F10.5 , r .’ * t E10»b)

FORMA'IM*' ФУНКЦИИ ЙЛИЯНИМ*’ ) •

FORhATO * ,4л,'’КОЭФФИЦИЕНТЫ ФУНКЦИИ ОИОИЯ')

FORMATl” '.ЗоХ,A2,5X,К14*5) ; >

FORMATf' ',4X('ПРЕДЕЛЬНОЕ ЗНАЧЕНИЕ Н£ИСКДЮЧЕННОЙ')

ОRMAT(',4Х , 'ДОПОЛНИТЕЛЬНОЙ ПОГРйЬНОСТИЛ'.1ЗХ,F6,3)

PRINT 500 •••••’ .

PRINT 400 ’ ;

PRINT 003

PRINT 417.FV

PRINT 423

PRINT 404.KSI1,KSIIU

PRINT 400 .. .

FORMAT! ' ' ,45X , 'Q0' )

FORMAT! ' ' ,.45Л ,'Ў0+Ц1Е' )

FORMAT!' ',45X,'40+Q1E+Q2b**2') •

FORMAT!' ' ,>45X ,'Q0+Q1 E+Q2EX' ) -

FORMAT!' ',45X,'ЦЯ+Ц1E+Q2E**2*X+t3EX')

■ FO tiMAT! ' ' ,45X , 'Q0 + QlE + O.2t/**2*X+C3t;X + O4*E**2*X' ) GOTO(3b,36.37,36,39,40) ,7!FV ,

PRINT 333

>';

GOTO 42.

PRINT 41 ■

GOTO 42 1 -


8ъам_иньл) Псдлиса и 5^/па

і я-м


! /Ыгт и дата 1

'Чл'

35

36


X

от

>

ед



>


ГОСТ 2. ' 68' CP qстm f I

' .■ . , ..Avi.. - . x мх-^ЪЙЧ—. Лі-І : . .Г*. '■ * _ > «w -■ WWV . I || | ,x ...37 PRINT 51