poly_pow.pro

来自「IDL语言编写的用于天文自适应光学仿真的软件CAOS V6.0的第一部分。」· PRO 代码 · 共 24 行

; $Id: poly_pow.pro,v 1.2 2002/03/14 11:49:11 riccardi Exp $

function poly_pow, p, power, DOUBLE=double
;+
; POLY_POW
;
; coeff_out = poly_pow(coeff_in, e [, DOUBLE=double])
;
; returns the coeffs of the polynomial given by the e-th power
; of the polynomial of coeffs coeff_in
;
;-

if power eq 0 then begin
    if keyword_set(double) then return, [1d0] else return, [1]+0B*p[0]
endif

p_out = p

for k=2,power do p_out = poly_mult(p, p_out, DOUBLE=double)

return, p_out
end