LSR DE 3 KUTTA-GILL  EXTRA PRECISION
zf1={61}	zn1={61}
0r,ita50r	ica16b+zn1+zn1	imr4b	its2b	ica16b+zn1	imr1b	imr5b	its3b
imr6b	iad2b	10r,its2b	imr108r	isu3b	iad16b+zn1+zn1	its16b+zn1+zn1
ica16b	iad2b	its16b	OUT	ao1r	20r,ao4r	ao13r	ao14r	ao15r	ao17r
su55r	cp57r	ca51r	clh16	ab1r	30r,ab4r	ab13r	ab14r	ab15r	ab17r
ca52r	clh16	ab2r	ab6r	ab8r	40r,su53r	cp59r	ca0	ts56r	cs54r
clh16	ab2r	ab6r	ab8r	IN	50r,isp	-zn1	+3	imr15b	+12
its15b+zn1	+0	IN	isp1r	ao56r	60r,clh15	cp64r	IN
ispzf1	IN	ica1b	imr112r	iadb	itsb	ispzf1	70r,its15b+zn1+zn1+zn1
IN	isp77r	+.	75r,ita107r	ica73r	its16b+zn1+zn1	OUT	ao77r
80r,su70r	cp71r	IN	ica110r	its6b	its9b	its12b	ica112r	its5b
its14b	90r,ica116r	iad118r	its15b	ics15b	its13b	ics9b	its4b
ica120r	iad122r	its8b	100r,ics8b	its7b	ica124r	iad114r	its11b	ics11b
its10b	isp	+3.	110r,+1.	+.5	0.74632	1.57777	0.52525	1.76777
0.52525	1.56777	120r,0.45373	1.76777	0.63143	1.53777	0.66501	0.01000
END OF SUBROUTINE
{61}m