Voorspellingen van het gedrag van de proefopstelling

Voorspellingen van het gedrag van de proefopstelling

Citation for published version (APA):

van Heck, J. G. A. M. (1982). Voorspellingen van het gedrag van de proefopstelling. (DCT rapporten; Vol.

1982.020). Technische Hogeschool Eindhoven.

Document status and date:

Gepubliceerd: 01/01/1982

Document Version:

Uitgevers PDF, ook bekend als Version of Record

Please check the document version of this publication:

• A submitted manuscript is the version of the article upon submission and before peer-review. There can be

important differences between the submitted version and the official published version of record. People

interested in the research are advised to contact the author for the final version of the publication, or visit the

DOI to the publisher's website.

• The final author version and the galley proof are versions of the publication after peer review.

• The final published version features the final layout of the paper including the volume, issue and page

numbers.

Link to publication

General rights

Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights. • Users may download and print one copy of any publication from the public portal for the purpose of private study or research. • You may not further distribute the material or use it for any profit-making activity or commercial gain

• You may freely distribute the URL identifying the publication in the public portal.

If the publication is distributed under the terms of Article 25fa of the Dutch Copyright Act, indicated by the “Taverne” license above, please follow below link for the End User Agreement:

www.tue.nl/taverne

Take down policy

If you believe that this document breaches copyright please contact us at:

openaccess@tue.nl

providing details and we will investigate your claim.

- & -

I

f

-f

s f

ip

i

;L I

i

J =

t

Y

"

i X 3 Y

s

2 FUNKTIE Z.t r , L STeMDWWT 5 . 5 % -3.58 Z.58 R t B B E K i J K U L a K 3.88 VERBORGEN L I J M E N CESTIPPELD RCHTERULhKKEN RLS D R M D F I C U U S GETEKEND P E R S P E k T I E F TEKEHINC X (DEflPiMC

B )

ff1 n RdX

D E R P I N C S U E R H O U D I ~ C ZETB RODE 1 Y (STIJFHEID K ) x v

z

STWWUN? -5.59 4.89 2.89 RIBBE O(XJKULf3K 3.98

ACHféRULIPKKEM &LS Q f f ~ ~ Q F I G ~ U R GETEUEND PERSPEKTIEF TEKEWING

- 1 8 -

A R C (PHItlIl I I ' Y (STXJFHEXD R )

,

I

I

' Y (STXJFHEXD R ) X Y 2 STAFIWUOIT 5.58 -3.58 2.S8 RIBBE KIJKULAK 3.863

VERBORGEN LIJNEW QEGTXPPELD

&CHTERULhKKEN @LS D R ~ A D F I Q U ~ ~ OPEKEND PERSPEKTIEF tEKENlN0 PWASE tPHXC13) RODE I Y V X Y FUPIKTIE PIIN Raw 7 ST4t4DPUNT 5.së - 3 . 5 8 2 . G RIBBE K I J K U L A K 3.88

UERBORCEti LIJNER GESTIPPELD

ACHTERWLMKEH ALS DRAQDFICUUR CETEKEPID

PERSF'EKTIEF TEKEr(1YG

X (DERPING B I

:<

f

4 9 -

I

I I V (STIJFHEID K ) X (DEAPINO €3) X Y t STANDCOWT 6.58 - 3 . 5 8 2 . s ~ RIBBE KIJKULIIK 3 . 0 8 'JERBORQEFI LIJNEN O E S t X P P L L D

ACHTERULAKKEt4 ALS DRIIADFICUUW OEPEKEND

P E R S P E K T I E F TEKENING R*X x Y FUPIKfIE 8. X Y

z

STANDPUFIT 5 - 5 8 4 . 0 8 2 - 0 8 RIBBE K I J K U L R K 3 . 0 8

VERBORGEN LIJNEN CESTIPPEL3

ACHTERULAKKEH 4Lc3 DRAADFIGUUR GETEKEND

MAC ( P H I C 3 3 ) RODE ! Y (STIJFWEID K ) X Y

z

SPANDPWT 5.5% -3.m 2.58 R I B B E KIJKUL@K 3.88

VERBORGEN LIJNEN OBSTIPCELD

ACHtERULûKKEN AL8 D f f A A D ~ I 6 U ~ ~ GETEKEND

P E R S P E K T I E F 'IEKENINQ V (StIJFHEID Bo X (DERQINQ I)) X Y

z

STWOPUti? 5.58 3.48 2.98 RIBBE KIJKULkK 3 . 8 8

UERBOROEW LIJNEW CESTIPPEL3

RCHTERULAKKEN c)LS DRRADFICUUR CETEKEVD PERSPEKTIEF TEKENING

-22.-

______-_--

-

_t

- 2 3 -

FREOUENTKE RODE 2 (STIJFHEID K x Y STmDPU)(T RIBBE KIJUULBiK 3.88 U E U ~ O ~ G € M LIJREM OEOTIPPELD PERSPEKTIEF EN

ILKEWI#Q

kLS ~ B i ~ D F I ~ U U V RIW ANLX X I .

a

6. 3 Y i. 7 i. 8 FUNKTIE 7. i?

a.

3

6STIJFHEiD

to

f

X (DENPING 8 ) X Y

z

STFINDPUNT 5.58 -3.58 2.56 R I B B E KXJKWLhK 3.88

UERBORCEN LIJNEN GESTIPPELD

ACHTERULOCKFN A L S DPaQDFTCIJIJR CETEICENP

DERPXt466UERWOUDINQ RODE Z Y o ( 9 f4XN RA)( X 1

*e

8. Y 1.8 i. CUNKTIE 7.1 t .

-sz

I

330U 2

íCI3IHd) 35UHd

- 2 1 -

9 (STIJFHEID K > X v 2 5.58 -3.55 2.58 WLM 1 . 9 5 LIJ?%ñ CES'TXPCELD YEP( A ~ ~ 1 6 U U R QE

act

ANDEREH, Z!*KIJKRICHTIN6 U E U ~ ~ ~ E R E N 3.PLOTTEW TE :KPHD M f M fl&# X 1.90 8 5. Y 1 . m 7 i. FUNKTXE B). i. Y (STIJFHEID K ) X Y 2 S? FIN WüN? 5.66 4.88 2.89 RIBBE KIJKVLAK 3.88 UERBOUGEN L I J N E N GESTIPPELD

ACHTERULFIKKEN ALS DRAfiDFIOUUR CETEKEND

RODE 2 X (DERPIMC

B)

x Y STACIWUNT UERBaRGEf4 LIJNEW

~

~

$

~

X

~

~

~

L

D

RCHTERVLAKKEN RLS ~ R ~ ~GETEKEND ~ ~ I Q ~ U ~ PERSPEKTIEF TEKENING RIBBE KXJKWLAK 3.eQ

A

X Y FUNKTIE v (STIJFMEID i[) X 1 2 STANDPUNT 5 . 5 0 - 3 . 5 8 2 . 5 0 P E R S P E K T I E F T E < E N l N C 1s RIBBE KIJKULAK 3 . 8 0 UERBORCEN LIJNEN G E S T I P P E L D RCHTERULALKEtí 6 L S DReADFICUUR C E Ï E K E B D

-39-

X Y N LXJNEN-OESTSPPELD hKKhN ALS D ~ ~ ~ ~ I Q U U ~ QETEKEND PERSPEKTXEF TEKEWXM n RIP( X 1. Y i. FUNCTIE 4. X (DEWIN6 B I x Y

z

STFiNDPUNf 5.58 -3.59 2 . 5 8 R I B B E KIJKULAK 3.08 VERBORGEN L I J N E N G E S T I P P E L D

6KHTERULAKKEN ALS DR6tFiDFICUUR GETEKEND

X V f

JKULM 3.88

t4 LIJPEN GESTIPPELD

~ ~ E f f U ~ K K ~ ~ ALS DRfiADFIOUUR GETEKEND

P € ~ 5 ~ € ~ T I ~ F PEKENINO

DI

X Y

z

StûNDBDWT -5.88 i . 7 ~

1.ee

RIBBE KIJYULM 3.88 E R B O R G E N LIJWEN GESTIPPELD

MHTERULPIKKEN ALS DRRI~DFICU~R GETEKEND

PMFILLEL PROJEKTIE

-33-

RAC

RüDE

X Y

U ~ R ~ LXJNEN ~ G ~QESTIPPELD N

ACHTEBULAKKEN AL9 DR6ADFIGUUR GETEKEND

PERSBWTICF TEKENING PITH X I . Y 1. FUNKTIE X v

z

STAHDPUWT 5 . 5 8 - 3 . 5 8 2 . 5 8 RIBBE KIJKULAK 3 . 8 8

UERBORCEN LIJNEN GESTIPPELD

ACHTERüLMKE!4 aLS DRfi&DFItUUR GETEKEND P é R S P E K T P E F TEKENING

X (DERPINP

B )

ffXN AAX X Y CUMKTIE

B.

( , PERSPEKTXEF TEKENING Y ( S T I J F H E I D K I X Y 2 STúffDPUNT -5.09 1 . 7 8 2.08 R I B B E YIJKVI6tK 3.88

VERBORGEN LAJNEN GESTIPPELD

GCHTERULAKYEtí fiLS DSAúDFiCUüR GETEKEND s A R & l L E L PROJEKTIE

-35-

V (STIJFHEID Io x BCHTERWLPIKKEN A i 9 ~ R ~ D F I ~ U U ~ GETEKEND PERSPEKTIEF TCKENXNO

fi

R I N X 1. v 11. FUNUTXE O. X Y 2 STfitiDPUNT

-s.ee

7 . 8 8 2.56 RIBBE K I J K U L ~ K 3 . 9 8

UERBORCEM LIJHEN GESTIPPELD

I K x ~ T E R U L M K E N RLS D R A A U F I G U U R GETEKEND

DERSPEKYIEF TEKENI%C

-6

3

7038

/o

m

I

-

- B i -

$SET INSTWATICN

BEGIN FILE IN(KI~REM3TE),UIT(KIND;DISK,~CSIZE=22), F77FO ( KINPDISK,MAXRECSIZl5 ) , MRXFO (KIND=DISK, p.w(REcsIzE=22) ,F77~TA(KI~DIS<,MRXRECSIZE;15),

MAxZETA(KIhD=DISK,MAxT(ECSIzE=22) ,F77F'HIlM(KImDIS(, r.lAxRECSIZE=15 ) , W H I l M ( KINWDISK,MWRECIZE=22), P77 PHIZH(KIND=DISK, bnx~~cSIZE=15 )

,

MAXPHI2M (KIND=DISK, MAXREcCIZE=22) ,P77~113M(KIND=DISK,MAXRECSIZFrl5),

W H I 3 M (KINWDISK ,MAXREG1 ZE=22 ) ,OUT (KIND=REM3TE)

,

F77 PHI

1

P ( KINWDISK ,MWRECSIZE=15 )

,

MAXPHI1 P( KIND=DISK, Iyu<HECSI;sE=22 ) ,F77PHI2P(KINBDISK,MN<RECSIZE=15), MIU<PHI2P(KIND=DISK,MAXRECSIZG22 ) ,F77PHI3P(KIND=DISK, b u u t ~ c S I ~ 1 5 ) ,MAxpHI3P(KINBDISK,MAxIIEcsIzE=22) ;

INTEXER I ,I1,12, I3,IB, IK,KP, J,W,POOLPLoT,'i"E;

REAL AE,FE,K,iU,Ml,M2.M3,Bi

FO[ 1 : 61 ,BPIM'[ 1 : 51 ,KE'IM'[ 1 : 51

,

FoPIM'[ 1 : 5 I 1 : 51 I

ZFXAPLI3P[1: 5,1:5] ,PHIlPLh[ 1: 5,1:5], PHILPLOTN 1:5,1:51, PH13PIM'M[1:5,1:5] ,P€IIlPLDW[

1

:5,1 :SI, PHI2PLDW[1: 5,l :SI, PHI3 PLoTP[ 1 : 5 , l : 5J i

u)MpLEX ARRAY PI: 1 :6] ,A[ 1 :6,1:3,1:3] ,PHI[ 1 :6,1:31;

BXILEAN AURAY AcC[1:61; DEFILE TO= STEP 1 UNTIL #i

ARRAY (;DEF[ O : 61

,

DC O: 4,l: 3,l: 31 ,WO[

1

: 61, ZETX

1

: 61 t

PK(~CEWRE RESIDUES( P , D , M ~ M ~ M ~ ,A,NPOLEN) ; cokP.lf>LM ARRAY P[*l,Ac*,*,*I:

€UP& MlM2M3;

mAL ARRAY D[*,*,*li INTEGER NPûLEN;

NGIN 1NTM;ER I,J,K,L,MODE;

FüR WDE:=l 'IQ NPOLEN Do BEGIN OOMPLM ARRAY PPP[O:43;

CDMPLEX W,PIP; PPPCO] :=l: PPPC11 :=PcM)DEI ; PPPC21:=PPPC11*PPPC11i PPPC3.j :=PPPC21*PPPClI ; PPP[4j :=PPPC3J*PPPCll i EüR 1:=1 1'0 3 W EDR J:=1 TO 3 D3 =IN DP:=~MPLEX(o,o) ; FDR L:=O TO 4 Do PIP:=(XIMPLIX(MlM2M3 I O ) ;

EUR L:=l 'IO NPOLEN Do IE L NED m D E THDi PIP:=PIP*(PCM)DEI-P[L1) i A[ MIDI:, I ,J] :=W/PIP : DP : =DP+OMPLEX( DEL, I, 51, O ) *PPP[ LI i END END END RESIDUES: PIXKEDURE EIGEN(A,P€II,NPOLEN) ;

CDPIPLEX ARRAY A[*,

*,

*I, PHI[*,

"1

;

INTGER NPOLEN: EGIN I t a u m WDE,I;

03MPLEx

tax;

FOR WDE:=l 'IO NPOLEN Do BEGIN M4X:=cDMpLM(O,O);

FüR I:=l 'IO 3 Do

FOR I:=l 'It.) 3 Do

IF CABS(A[M3DE.,l,I])>(aBS(MAX) 'RIFSJ EW:=ACM)DE, 1,Il;

-

8 2 -

OMPm

c:

EGIN

REAL

CPH;

C!PH:=AFKTAN( MiG( C )/REAL(C ) ) *45/ARC.TNJ (

1

) i

IF REAL(C) < O WMEN CPH:=BH+l80: IF CPH>180 W E N BH:=BH-360; B€WE:-qH:

END @€WE:

WRITE((xpp, </,"???MAX ITERATIES,ABS NAuhM,REL NAUW">) :

READ(IN,/,MAX,AE,RE):

Ml:=o.211;

M2 ~2.316 ;

M3:=1.094:

Kl:=1.77@8,

WITE ( UIT, <"Ml="

,

E12 .5,

/,

"MZ-"

,

E12 .5,

/ ,

"MS"

,

E12 .5,

/,

"Kl=" ,E12.5,

///>

,M1 ,M2 ,M3 ,Kl ;

MITE(CUT, <"WELKE POOL M)lT G W m WRDEB?" > ) ;

WITE(ouT, <"WELK TYPE PLOT WILT U?",

/,

$INCLUDE "I%X;RAPH ON USERS"

READ(IN,/,WoLPLOT);

"O=GEXN PIM'FILE AANWiKlW',

/,

"l==GRAFIEK OPIKIJWEN UIT PRISm'S",/,

"%GRAFIEK OPIKIJWEN UIT ONDER EN I 3 0 V ~ V L R K " , /,

"3=GRAFIEX OPECiJhEN UIT IIOVJBVLRK" > ) ;

READ(IN,/,TYPE);

EDR IB:=1 'IO 5 Co

EDR IK:=1 TD 5 Co PEGIN B:=lOW(I&1)*1225; K:=17 : EPIDJ?[ IB] :=B; KE'LM'[IK1 :=K; OEF[ O] :=Ml*M2%3 : CDEF[lJ:=2*Ml*M2*B+Ml*M3*B; 03EF[2J:=2W*M2*K

+

Ml*B*B

+

Kl*Ml*M3

+

03EF[3]:=2*B*PM1

+

2*B*Ml*KL

+

2*M2*B*Kl + 03EF[4]:=Ml*K*K

+

2*Ml*K*Kl

+

2*M2*K*Kl

+

CDEFC 51 :=2*B*PKl; Om6]:=K*K*Kl;

ZEFOPQL( 6 ,COEF ,AF: ,RE ,MAX, P ,ACC,KP) ;

FOR I:=l 'ID KP W IF NOT ACCC11 'IMEN P[I~:=Q)kPm(o,o)i

EDR 1 1 ~ 0 To 4 Co JBR IZ:=l Tü 3 Co JBR 1 3 ~ 1 ?o 3 Co D[Il,I2,131:=0; DC0,i ,i1 :=(K+2*Kl )*K; N0,1,2]

:=NO,

2 ,i] :=2*K*Kl: DCO, 1,3] :ID[ O, 3,1] :=K*Kl; D[O, 2,21:=2*I(aKl; N0,2,31:=Dc0,3,21:=@K1;

NO,

3,3] :=@KI ;

Ni,

1

,i] :=2*B*(K+KL i

M1,i

,2]

:=N

1,2 ,i] :=2*B*Kl; D[ 1,1,3]

:=M

1,3,13 :=B*Kl i K*Ml*M3

+

Kl%2*M3; M3*B*Kl: B*B*Kl

+

M3*@Kl:

-

6 3 -

-

8 4 -

D[1,2,2] :=2*B*K1;

a

1,2,3] :=D[ 1,3,2j :=B*m ; D[1,3,3] :=B*Kl ; D[ 2,1,1] :=2*M2*K+M3*ICl+M3*K+B*B: U[ 2,1,21 :=D[2,2 -13 :=M3*H; D[ 2,2,2] :=2*Ml*K+M3*a; D[2,2,3] :=D[2,3,2] :=Ml*K: D[ 2 , 3 , 3

1

: =Ml* ( K+Kl )+M2*Kl; DL 3,1,1] :=( 2*M2+M3 )*B; a3,2,2]:=2*Ml*U; ~€3,2,31:=~€3,3,21 :=Ml*B; Dl 3,3,3] :=MI *B; D[4,1,1] :=Mî*M3 : D[4,2,21:=Ml*M3; D[ 4,3,3

1

:=Ml*M2 ; RESIDJES( P, D,M1 *MZ*M3 ,A, KP) ; EüR I:=l 10 KP Do DEGIN W[ 11 :=cAI3s (PC I1 ;

EIGEN (A, PHI, KP ) ;

ZETa 11 :=-REAL ( P[ I1 )/WO[ I1 : FOCI] :=wOCIl/(8*ARCI'AN(

1)

)

ENU;

FOPLOT[ IB, 1K1 :=FO[poOLp~1; ZETAPIDT[ IB, IK] :=ZETA[FO3LPlDTI :

PHIlPLOTM[ IB. IK] :=CABS ( PHI[ wOLPL13T, 11 ;

PۓI2PILYI'M[ LB, IK] :=cABc (PtlI[wOLP~,2]) ;

PHI3PLoTM[ IB. IK] :=CABS ( ]?flIC pcx)LpLup, 3

1

i

PHIlPLoTPC IB, IK] :=CPHASE ( WI[ WOLPLDT,

11

;

PHI2PLI3TP[ IB, IK] :=CPHASE(WI[F03LPLOTt21) ;

PHI3PLoTP[ IB, IK] :=CPHASE( PHI[ PcXILpID', 31 i

WITE( UIT, <"K" ,El, "B"

,m,

"RE( Pc 11 ) " ,X2, "IM(PC 11 ) " ,X2, "FO[ I]",X5, "ZETM I]" ,X3, "ffi(PHI1)" ,X2,

"EI( PHI1 ) " ,X2, "MG( FHI2

1''

,

X2 I "PH( PHI2

1''

,X2,

"MG( Ai13 )" ,X2, "PH( PHI3

1''

I X2 /e

*

( 12

(m.

2 A 1 I / ) a

/

,

"RE(Al1)" ,X3, "IM(A11)"

,U,

"RE(A12

1''

,X3, "IM(A12)",X3 ,"RE(A13 )" ,X3 ,"IM(AI3)" ,X3, "IM(A22 )" ,X3, "RE (A23

1''

,X3, "IM(A23 )"

,U,

"RE (A33 ) " , X3, "IM(A33 )" ,

/

*(12(E9.2,Xi),/)>,KP,FOR I:=l To KF' W [K,B,P[II, FO[I],'ZETA[I],FûR J:=1 TO 3 Ly) [CABS(PHI[I,J1), CPHASE( PHI[ 1,

JI

) 3

1,

KP,FOR 11:=1 TO KP M) [FOR I2:=1 TO 3 Do [EDR 1 3 ~ 1 2 TO 3 Do [ A ~ I 1 , 1 2 , 1 3 ~ ~ ~ ~ ~ ;

"RE(A22

1''

,X3,

EI\ID;

~ ( U I T , C R U N C H ) ;

~I?e(OUT,<7(5(E12.5,X1)/)>,EUR IB:=l TO 5 Ly)[EiPwr[IB11,

8

%

% EUR IK:=l TO 5 Ly) [KPLOTCMJJ,FOR IB:=l TO 51 M)

[FOR IK:=1 To 5 Co [EY)PLlYKIB,IK111); IE' TYPE>O AND TY!?E<4 ?Hili

&GIN GWM(BPLOT,WLOT, FOPLOT,5,5,F77EO ,MpXED ,TYPE, GRAPH( BPLOT

,

KPIx>T, ZETAPLOT, 5,5 ,F77ZEI'A, MhxZETA,TYPE, GRAPII( BPLOT,KPLOT,PHIlPLOTM,5,5 ,F77P€lIlM,MPXPHIlM,TYF'E,

GRAPH(BPLOT,KPLOT,PHI2PLOTM,5,5,E'77PHI2M,MAXPHI2M,TYPE, GRAPH(BPLOT,KPLOT,PHI3PLOTM, 5 , s ,F77PHI3M,MTUPHI3M,TYPE,

PRIS~,WRITEFILE,O,O,O,~,O,O,FALSE) ; PRISM,WRITEFILE,O,O,O,O,O,O,FALSE);

PRISM,WRITEFILE,O,O,O,O,O,l ,TRUE) ;

PRlS~,wRITEFILE,O ,O,O ,O ,O, 1 ,TlüJE) ; PRISM,WRITEFILE,O,O,~,~,~,~ ,TRUE) ;.

END;

m.

GRAPH( BpL13T, KPIDT, PHI1 PI13Tp, 5,5 ,F77P€lIl P ,MAXPHI1 P,TYL'E I

GRAPH( BPUX ,KPLDT, PHILPLOTP ,5,5 ,F77PIlI2P,MAXPHI2P ,TYPE, GRAPH( BPLOT ,KPLOT, PHI3PIDIF, 5,s ,F77PHI3 P ,MNCPHI3P ,TYPE,

PRIst.R,wRITEF1LE,0,0,0

,o,o,o

,FAT.JSE) :

P R I ~ , W R I T J % F I L E , O , O , O , O , O , O , ~ ~ ~ ) ;

PROCEWW PRI~(X,Y,A,MIM<,MINY,MINA,MA)(X,MAXY.MAXA, IX,IY,FL,TYPE,WRITEFILE) : REAL MINX,MINY,MINArMIV(X,MAXY,MAXA; PIXxtXUm WRITEFILE; FILE EZ; INTEGEI IX,IY,TYPE; AKMY XL*I,Yc* J,AC*,

*I;

BEGIN ARRAY uX)RDC1:6,1:31;

INTEGEX 1,J;

FEAL CX,CY,CA,PX,PY,PA; I F Ivpyo[ NEQ MINX 'EEN

C X : = ~ / ( M - M I N X ) ELSE

0<:=1@24,

IF W X Y NlXl MINY TPIEN

I F PAXA NDJ MINA THEN

Cri:=l/(MAXA-MINA) ELSE cA:=1@24;

1~

:=-MINx*Q(; PY:=-MIWCy; PA: =-MINA*cA; WRD[

1 , l l

:=X[ M J*CX+PX ; a X R D [ l , 2 J :=YCIYl*CY+PY; W W [ 2,13

:=xc

I x + 1 j

* a + m

; ~ R D C 2 . 2 3 :=(xxIRDc 1,21;

mmc

3,11 :=ux>RDc2,11: W R D [ 3 , 2

1

: =yc IY+l1 *cY+PY;

LY:=~/(MAXY-MINY) ELSE

m:=l

a 2 4 ,

iDR I:=1

Si'm

1 W I L 3 Co E I N 030RD[ I + 3 , 1 1 :=OORD[ I,

11

;

mmc

1 + 3 , 2 1 :=oxIRDc I, 2 1 ;

WRDC I,

31

:=O; END;

uxIRD[ 4,31

:=x

Ix, IYI*CA+PA;

@3ORD[5,3J:=AlIx+l,IY]*Cn+PA;

WRJl[ 6,3] :=A[ M+1, IY+l]*CA+PA; WITEFILE

( m m ,

FL,TYPE ) ; mRDI: 1 , 2 1 :=Yi IY+11 *CY+PY;

WOWC 2 , l l

:=xc

MI*CX+Ex ;

W R D [ 2 , 2 J :=YCIY+l]*E+E>Y;

ux>RD[ 3 , l ~ : = W R D [ 2 , 1 J ;

W R D C 3 , 2 1 :=YCIYI*CY+PY;

EDR I:=l STEP 1 W I L 3 D;>

BEGIN W R D [ I+3, l i :=WRtf I ,

11

;

ux)m[

1 I 1 J

:=x[

=+11*(3<+m

; W R D C 1+3,21 :=mRDc I. 21 ; 030RDC I,

31

:=O; END: ~RD[4,3]:=A[M+l,IY+lI*CA+PA; W R D [ ~ , ~ ] :=ALM ,N+l]*cA+PA; WRDL6,31:=A[Ix ,N I*CA+PA; MWEFLLE(C~~KD,FL,TYPE) ;

mNL,

P N m ; PKKEWRE WR~~TEFILE(WRD,FL,TYPE) ; FILK EL; ARRAY 03oRD[*,*]: IMI'IXEX TYPE; E G I N IMFGER l , J ; CASE TYPE OF BEGIN

1: WIN WWENT VOLLEDIG PRISW;

W I T E ( F L , <6( I10,3 (F12.8 , X l ) / )

,

"-1"

,

/,

"5 1 3 2 O"/

"5 4 5 6

o"/

" 5 2 3 6 5 O"/

"5

"-2">,FOR 3 1 4 6 I:=l O"/ STEP 1 UNI'IL 6 Co

"5 1 2 5 4 O"/

CI,FOR J:=I STEP i u N r I L 3 DO

[cDoRD[I, J l l 1 ) :

END VûLLEDIG PRISm;

2: í33GIN COiNEN" ALLEEN E3C34 EN CNDER; hRITE(EL, <6( I10J (F12.8 ,X1)/ )

,

"-1"

,

/,

"5

1

3 2 O"/

" 5 4 5 6 O"/

[I,FOR J:=l STEP 1 UM'IL 3 Co

cQ)oRD[I,JI1J);

"-2">,FOR I:=l STEP 1 UNTIL 6 Co

END BOVEN EN QM>ERviAK;

3 : BEGIN

cxx+mT

ALLEFN B O V E N V M ;

W I T E ( F L , <6( 1 1 0 , 3 ( F 1 2 . 8 , X l ) / ) , "-l",

/,

" 5 4 5 6 O"/,

"-2">,EDR I:=l STEP 1 UNTIL 6 Co

[I,FOR J:=1 STEP 1 UM'IL 3 Co

C O R D C I . J J J J ) ;

END ALLEEN B3vEMILAK:

END W E TYPE OF; END WRITEFILE;

PKKZWRE GRAPH ( X, Y ,A, NX ,NY FL^ ,FU ,TYPE, PRISM ,WRITEFILE,

XMIN,YMIN,AMIN,X,YT%X,W,EXPAND) ;

DOOLEAN EXPAND;

RE% XMIN,YMIN,AMIN,XMW,YI"RX,M; ARRAY XC*I,YC*I ,A[*, *I;

INEGER NX,NY,TYPE;

FILE FU,FL2:

PRcaDuRE PRIrn,wRITEFILE;

BEGIN INiBZER IXJY;

REAL CA,PA,MINX,MINY,MINA,MAXX,MAXY,MAXA;

I F EXPAND WEN

BEGIN I F X M I N NEQ O OR XWX NEQ O ?HEN

FüR IX:=l STEP 1 W I L NX Do

BEGIN I F X[ M]>XMX ?HEN X[

Ixl

:=M:

IF X[M]<XMIN 'R1EN X[Ixl:=XMIN;

END ;

IFYMIN NEQ O O R W NEQ O liE%

FüR ïY:=l STEP 1 W I L NY Do

BEGIN I F Y [ N 1 > Y N U ?HEN Y[Iï1:='iMU;

I F Y[IY]<YMIN WEN Y[IYI:=YMIN;

END;

I F AMIN NEQ O OR W NEQ O 'WEN

FDR

IX:=1

STEP

1

W I L NX Do

FDR IY:=1 STEP 1 W I L NY Do

BEGIN I F A[M,IY]>AWX RlEN dIX,IYJ:=AMRx;

I F A[M,IY]6AMIN dIX,IYI:=AMIN: END; END: M I N x : = M : = X [ l ] ;

MINY:=M:=yc11;

MINA:=MAxA:=Ac

1,il:

KIR I X : = ~

srm

i WIL NX

w

BEGIN I F X [ M ] > M ?ZIEN EIRxx:=aIX];

I F X[M]<MINX THEN MINX:=X[IXI;

- 83

-

I;DK IY:=l STEP 1 W I L M!

BE(;IN IF Y [ N I > M "HEN

m : = y [ : m 3 ;

IF Y[IY]<MINï lliE?i MINY:=Y[IYl;

kul? IX:=l STEP 1 UNTIL Nx a,

BEGIN IF A[M,IYI>E.PXA THDi W = A [ lX,IYI;

IF A[IX,IY]<MINA THEN MINA:=A[.M,IYI;

END;

END;

IF EXPAND

nm

BEGIN

IF XMIN NfU O OR XWAX NEQ O ?HEN

EFGIN M I N K : = X M I N ; W : = m END;

IF YNIN NEU O OR Y N U

Nw

O ?HEN BEGIN MINY:=YMIN;M:=W END: IF AMIN NFQ O OR A W U Na2 O ?HIPI

I3EcIN MINA:=AMIN;bM:=M END;

mu;

EDR ix:=i

srw

1 UNTIL NX-i DO

FOR IY:=~

srw

1 UNTIL NY-1 D3

PRISM ( X I Y ,A, MINX ,MINY ,MINA, bw[x ,W ,MAxA,

Ik' WìXA NEQ MINA THEN CA:=l/(MAXA-MINA) HLSE

pE\:=-MINA*CA; I X

,

IY ,FLi ,"PE,WRITEFiLE) ;

cA:=181q;

mITE(FLl,<"l

o o

O"/ "2

o

1 O"/ "3 1 1 O"/ "4 1

o

O"/ " 5

o o

i"/ "7 1 1 1"/ " -1"

/

"3 1 4 3 2 O"/ " 3

1

5

o"/

1' -2"

/

"1 1

o

o"/

"-i"/ "3 1 2 @'/ '1 -2"

/

"1

o

1 O"/ "-l"/ " 3 1 2 O"/ "-2"/ "1 1 1 O"/ "-1"/ " 3 1 2 O"/ " -2"

/

" 6 O 1 1"/ " 8 1 O i"/ "2 1 O ",F19.4/ "2 O 1 ",F19.4/ " 2 1 1 " ,F19.4/ "-3"/ > ,A[NX ,11*CA+PA, AC1 ,NYI*CA+PA, A[NX,NY]*CA+PA) ;

I~~ITE(EU, ~ ( ~ 1 9 .12,ii, ,MINX,MINY,MINA,MAXX,MAXY,MPXA) : Lcxx(FLl ,CRUNCH) ;

LCXX(FL2,CRUNQ-I);

- R I -

-

$33"

I N S ~ I C L \ I

BEGIN

FILE IN(KIND=/I%X'E)

, U I T ( K I B l 3 = D I S K , P C S I Z E + 2 2 ) ,

F 7 7 F O ( K I N D = D I S K , M S I ~ l 5 )

,pIIAxFo(KIM)=DISK,

MAXZETA

( K I D D I S K ,

MAXRECS

I

!222=22

)

,

F77 PHIlM

(

KIkJDISK,

MmRECSIZE=l5)

, M A X P H I l N ( K I ~ D I S K , M A X R E C S I Z E = 2 2 ) ,

F77PHI2M(KIND=DISK,

E4AxREcsIzE=15

)

,

MAXPHï2M

(KIND-DISs(,

N?W?Ec=SIZE=22

)

, F 7 7 P H I 3 I . I ( K I N D = D I ~ , ~ ~ 1 ~ = 1 5

)

,

N;W;PHI3M(KIND=DISK,MAXRECSIZE=22)

,OUT(KIDREF.lE),

F77PK11P

( K D J P D I S K ,

MAxRECSIzE=15

)

,

MAXF'HIlP

( KlND=DISK,

~ ~ C S I Z E = 2 2 ) , F 7 7 P H I 2 P ( K I M D = D I S K , M A X R E S I ~ 1 5 ) ,

MAXPHI2P(

KINI)=DISK,MAXEECSIZE=22)

,F77PI-II3P( KIDDISK,

Mz1xTzEGIzIE=15)

,P.IAXPHI3P(KINC=DISK,PIAXFE~I~22)

;

INTEG=

1 , 1 1 , 1 2 , 1 3 , 1 B , I K , K P , J , M A X , ~ L E ' ~ , ~ E ;

REAL

AE,RE,K,Kl,Ml,M2,M3,B;

REAL

ARFUiY

CXIEF[0:6],D[O:4,1:3,1:3],WO[1:6],ZET~1:6],

FO[l:6] ,BPLMT[1:5] ,KPLCYT[1:5] ,FOPLOTCI:5,1:51,

Z F 2 A P l M f

1

:

5 , l :

51,

PHIlPI&$l:

5,l:

51

,

PHI2PlXYTP$1:5,1:

52,

P H I 3 P L O T M ~ 1 : 5 , 1 : 5 ] , P H I 1 P W T P ~ 1 : 5 , 1 : 5 ] , P H I 2 P ~ ~ 1 : 5 , 1 : 5 ~ ,

cyIp/sipLM

XFRY

P[1:6],A[1:6,1:3,1:3],PHI[l:6,1:31;

E93oLEAw

kfaRAY

ACCc1:61;

DEFIT33

Tc)=

STEP

1

W I L

#;

PROCEDURE ESIDUES(P,D,MIMZG

,A,NPOLE;?J)

;

CDPIPLEX ARRAY

P[*],A[*,*,*];

REAL

Mlprsm;

F

G

m

L

ARRAY

D[*,*,*I;

WTEGER

NPom;

BEGIN INTEGEIR

I,J,K,L,MODE;

p/IAxREcL>IZE+22)

, F 7 7 Z E T A ( K I N D = D I S K , ~ ~ ~ I ~ - l 5 )

PHI3PwfPC

1

:5,1:

51

;

FOR

PDDE:=1

To

NPûLEl54

Dc

BEGIN o0MPL;Ex

ARRAY

PPP[O:43;

ODMPLM

DP,PIP;

PPPcl]:=PcE!DDEI;

PPPCO1

:=l;

PPPC2f :=PPPCi~*PPP~lI:

?P?[3>

:=PPPf2]*PPPL1!

;

P P P [ ~ ~

: = P P P [ ~ ~ * P P P ~ ~ ] ;

FOR Ls=l

m

a

Do

FOR

J:=l

TO

3

Do

=GIN

DP:=CDMPLEX(O,O);

FDR

L:=O

m

4

Do

PIP

:

=a34PLFx

(Mlmm, o

DP

:=DP+û3WLEX (

DLL, I ,

JI

,

O

)

*PPP[L]

:

;

FUR

L:=l

To E3POLEN

Do

I F L

NEQ

M9DE

T!3N

PIP:=PIP*(PDoDE~-PCLI);

A

[

PDDE

,

I

,J

J

:

=DP/PIP

;

-

8 2 -

PE€t[Y'DDE,

I]

:=ACEODE,

1

,I]/p.py(;

E m

END EIGEE;

REAL

PROCEWRE

B H A S E ( C )

;

m w m

c;

BEGIN

R"0AL

B H ;

BH:=MCL%X( IPaG (C)/FEXL (C ) )

*45/ARCTAN

(

1

)

;

IE'

REAL(C)

< O THE? CPil:=BH+180;

IF

CPH>180

TBEN

CPN:=CPH-360;

B W E

:=BH;

$JNCUIDE

"DDGRlipH

ON

U S E S "

END

BHASE;

kRITE(OUT,</,"???1\IPhX

ITEEATIES,A€S

N A W i , R E L

NAUWK">);

READ(m,/,m,AE,RE);

Ml:=0.211;

FL?

~ 2 . 3 1 6

;

M3

:=1.094;

HRITE(UIT,

<''i'fl=''

,E12.5

/,

"M2="

,E12.5,

/,

"N3="

,E12.5,

/*

"KL=",E12.5,///>,Ml,P42,P13,Wf

:

WRITE(

m,

<''WELKE

POOL

P P X ~

GEPLOT

homm?"

> )

;

rn(IN,/,PCOLPIMT)

:

WITE(CYIT,

("WELK

TYPE

PLMT

WILT

U?",

/,

"O=GEEN

PLOTFILE

A A l ~ " ,

/,

'I

l=GRAFIEK

OPBXJWEN

U I T

PRISER

' S"

,

/ ,

" 2 - G W I E K

0PBc;Uc.dEM U I T

ONDER

EX

BWE8VM",

/,

" 3 = G W I E K

OPlXUWEN

U I T

BovENvLs1K"

>

1

;

EIFWD(IN,/,WE);

FOR

IB:=l

'ID

5

Do

FOR

IK:=1

rn

5

co

E G I N

B:=lOW(IB-1)*1225;

K:=17 :

BPLOTC IB] :=B;

KPLMf[

IK]

:=K;

CDEFL

11

:=2*MI*E2*B+I'4d*M3*B;

WEF[21:=2*r;P1?42*K

4-

Hl*B*B

-i-

K l % l l * H 3

+

uIEF[ 31

:=2*B*K%l

i-

2*B*Ml*Kl

i-

2*&%2*B*kM.

+

CX1EF[43:=Ml*K*K

+

2*Ml**KL

+

2*M2*K*KL

+

CDWI

51

:=2*B*K*Kl;

CDWC61 :=K*K*Kl;

FüR

I:=l

TO

Kp

Co

I F . N O T ACC[I]

33EN

P[I~:=UX~!?LIEX(O,O);

FOR

Il:=O

TO

4

Co

EUR

I2:=1

TO

3

Do

FOR

I3:=1

TO 3

Do

D [ I l

,

12,132

:=O;

D[ O, 1,1]

:=(

K+2*K1.

) *K;

D[

0,1,23

:=D[

O,

2,1] :=2*K*Kl;

D[O, 1,3]

:=D[O, 3,1] :=K*KL

;

D[0,2,21:=2*K*Kl:

DE0,2,31 :=D[0,3,2]

:=K*Kl;

D l O, 3,3

1

:=K*Kl;

D[ 1,1,11:=2*B*(K+Kl)

:

D[

1 , 1 , 2

1

:=D[ 1,2,11:=2*B*I(I

;

D[l,

1,33 : = D n ,

3,lI

:=S*Kl;

aEF[v!

:=Pa142?4?3;

K*-Ml*M3

+

KL*m*l43;

m*B*m;

B*B*EU

+

M3*K*KL;

GRAPEI(BPWT,KPWT,PHIlPWI1P,5,5,F77PHIlP,~~lP,TYPE,

GRApH(BPLMT,KpwT,PHI2PlDTP,

5,s

,F77PE3ISP,MAXPHI2P,'lTPE,

GRAPH (

BPLMT

,

KPWT,

PHI3 PïL3TP, 5

,5

,F77PHI3

P

,t.ITuu?HI3

P

,

TYPE,

PRIs~,wRImFILE,o,o,o,o,o,o,FALsE)

:

PRIS~,WRITEFILE,O,O

,O,O

,O,O,FALSE)

: