dclog.f

UNIX v6源代码 这几乎是最经典的unix版本 unix操作系统设计和莱昂氏unix源代码分析都是用的该版

	function dclog(z1)
	double complex z1, z2, dclog
	double precision a,b,c
	double precision dummy(2)
	equivalence (a,z2,dummy(1)), (b,dummy(2))

	z2=z1
	c=dcabs(z2)
	b=datan2(b,a)
	a=dlog(c)
	dclog=z2
	return
	end