File: EDYCHI.FC of Tape: Various/ETH/fc1

(.FC,.FO Focal format converted to 'WRITE' listing) 

---

C-OS/8 FOCAL, 1972

01.01 C EDYCHI
01.02 C AC SUSZEPTIBILITAET IN ZYLINDER
01.03 C 19.9.74
01.10 O O EDY
01.20 F X=1,.2,8.8;D 2
01.30 S S=0;T S
01.40 O C
01.50 O O TTY:,E
01.90 Q

02.01 C CHI=M1+I*M2=RR*COS(PH)+I*RR*FSIN(PH)
02.10 S N=2;D 3.1;S N=0;D 3.2
02.20 S PH=PH(1)-PH(3);S RR=RR(3)/RR(1)
02.30 S M1=RR*FCOS(PH);S M2=RR*FSIN(PH)
02.32 S Y=1/(FATN(PH))
02.33 S X1=1.46*Y+.36
02.34 S X2=X1-.01/(Y-1)
02.35 S PH=PH*180/PI
02.40 I (PH)2.45,2.5,2.5
02.45 S PH=PH+360
02.50 T %5.03,X,PH,!

03.01 C HILFSROUTINE
03.10 S A=1;S B=-1;D 7;S BR(N)=S;S A=0;S B=2;D 7;S BI(N)=-S;D 4
03.20 S A=0;S B=0;D 7;S BR(N)=-S;S A=1;S B=-1;D 7;S BI(N)=S;D 4
03.90 S PH(N+1)=FATN(BI(N)/BR(N));S RR(N+1)=FSQT(BR(N)^2+BI(N)^2)

04.10 S P=FATN(BI(N)/BR(N));S RR(N+1)=FSQT(BR(N)^2+BI(N)^2)
04.20 I (BR(N))4.3;S PH(N+1)=P;R
04.30 I (BI(N))4.4;S PH(N+1)=P+PI;R
04.40 S PH(N+1)=P-PI;R

07.01 C KELVIN FUNKTIONEN
07.02 C BER(X)+I*BEI(X), ORDNUNG N
07.03 C REIHENDARSTELLUNG NACH ABRAMOWITZ FORMEL 9.9.9
07.04 C GEGEBEN: X,N,A,B
07.10 S J=0;S S=0
07.20 S K=2*J+A;S Z=J+(A+B)/2;S Y=N+2*K;S G=(X/2)^Y*(-1)^Z
07.30 D 7.8;D 7.9;S G=G/KF/NF
07.40 S S=S+G;I (FABS(G)-FABS(S)/1000)7.5;S J=J+1;G 7.2
07.50 R
07.80 S KF=1;F I=1,K;S KF=KF*I
07.90 S NF=1;F I=1,N+K;S NF=NF*I

09.10 A ?N,A,B,X?;D 7;T %5.04,S,!
09.20 G 9.1

---