fc TAPE 131-238-923 ABRAMOWITZ MODIFIED BESSEL FUNCTION K0(X)
(24,6)
a1,+0.1	a2,+0.0000001	a3,+0.1	a4,+0.05	a5,+10.0	a6,+0.0	a7,+0.0	a8,+0.0	a9,+0.0
a10,+1.0	a11,+0.0	a12,+0.0	a13,+0.0	a14,+0.0	a15,+0.0	a17,+0.5	a18,+0.0	b2,icaa1	isua5
icpb1	STOP	b1,icaa7	ispm1+5	imra1	ispk1	itsa11	icaa10	idva11	itsa12
iada8	itsa8	icaa2	isua7	icpa13	icaa12	imra2	itsa6	icsa12	imra17
iada8	itsa8	a13,icsa12	iada6	icpa14	icaa7	iada4	itsa7	ispb1
a14,icaa1	iTOAp28.123t	icaa8	imra4	iTOAp28.1234567c	icaa1	iada3	itsa1
ispb2
k1,
LSR FU1 EXPONENTIAL SUBROUTINE
||||||||||
za1=a14
0r,ita38r	OUT	cs2042	su39r	ts54r	cp7r	sp9r	cs39r	ts2042	IN
10r,its+za1	idv52r	iad40r	imr+za1	idv50r	iad40r	imr+za1	idv48r	iad40r	imr+za1
20r,imr46r	iad40r	imr+za1	idv44r	iad40r	imr+za1	imr42r	iad40r	imr+za1	iad40r
30r,isp33r	IN	imr+za1	its+za1	OUT	ao54r	cp31r	IN	isp	+2
40r,+1.		+.5		+3.		+.25		+5.	
50r,+6.		+7.		+0
||||||||||
END OF SUBROUTINE
{61}m1,
LSR FU5 HYPERBOLIC FUNCTIONS
za1=a13
0r,ita53r	OUT	ca87r	sp8r	+0	ita53r	OUT	ca86r
ts51r	cs2042	10r,su88r	ts4r	cp17r	cs51r	ad87r	cp54r
sp19r	cs88r	ts2042	IN	20r,its+za1	idv84r	iad72r	imr+za1
idv82r	iad72r	imr+za1	idv80r	iad72r	imr+za1	30r,imr78r	iad72r
imr+za1	idv76r	iad72r	imr+za1	imr74r	iad72r	imr+za1	iad72r	40r,isp43r
IN	imr+za1	its+za1	OUT	ao4r	cp41r	IN	ica72r	idv+za1
50r,iex+za1	+0	imr74r	isp	IN	its+za1	imr+za1
idv82r	idv84r	iad72r	60r,imr+za1	imr+za1	imr78r	idv80r	iad72r
imr+za1	imr+za1	imr74r	idv76r	iad72r	70r,imr+za1	isp53r	+1.
+.5	+3.	+.25	80r,+5.	+6.	+7.	iad+za1	isu+za1	+2
END OF SUBROUTINE
||||||||||
{61}i START AT b2
v