00030 :
00031 n=len(d)
00032 theta=2*math.pi*fdt*(Numeric.array(range(1,n+1), Numeric.Float))
00033 zr=1.0-Numeric.dot(Numeric.cos(theta), d)
00034 zi=Numeric.dot(Numeric.sin(theta), d)
00035 return xms/(zr*zr+zi*zi)
00036
00037
if __name__=="__main__":
