Приложение 2 Листинг программы вычисления значения Qiпо формуле (2.22)

(let ((x 0))

(defun fact(n)

(cond

((< n 2) 1)

(t (* n (fact (- n 1))))

)

⁾

(defun C(n i)

(cond

((< n i) 1)

(t (/ (fact n) (* (fact i) (fact (- n i)))))

)

⁾

(defun deg(n p)

(cond

((zerop p) 1)

(t (* n (deg n (- p 1))) )

)

⁾

(defun Q2(n_ l p n)

(let ((outer_sum 0) (inner_sum 0) (inner_term 0) (term 0) (i_max 0) (prod_C 0) (r_max 0) (i 0)) (cond

((< n_ (+ p (- l 1))) 0)

((> n_ (* p l)) 0)

((eq l 1) (deg n p))

(t

(setq outer_sum 0)

(setq i 1)

(if (< i (\ n_ p)) (setq i (\ n_ p)) nil)

(setq i_max (\ n_ l))

(loop

(if (> i i_max) (return outer_sum) nil)

(if (< (- n_ (* i l)) (- p i)) (return outer_sum) nil)

(setq prod_C 1)

(dotimes (j i prod_C)

(setq prod_C (* prod_C (C n l)))

)

(setq term (* prod_C (C p i)))

(setq r_max (- l 1))

(if (> r_max (- n_ (* i l))) (setq r_max (- n_ (* i l))) nil)

(setq inner_sum 0)

(dotimes (r r_max inner_sum)

(setq inner_sum (+ (Q2 (- n_ (* i l)) (add1 r) (- p i) n) inner_sum)) )

(setq term (* term inner_sum))

(setq outer_sum (+ outer_sum term))

(setq i (add1 i))

)

)

))

⁾

(defun sum_Q2(n_ n)

(let ((i 0) (s 0))

(setq s 0)

(dotimes (i n s)

(setq s (+ s (Q2 n_ (add1 i) n n)))

))

⁾

;;(defun P(n_ n l s) (/ (Q2 n_ l n n) s))

(defun T2(n n_ s)

(setq tt2 0)

(setq l 1)

(loop

(if (> l n) (return (/ tt2 n_)) nil)

(setq tt2 (+ tt2 (* l (/ (Q2 n_ l n n) s))))

(setq l (add1 l))

)

⁾

(defun main(n)

(let ((n_ 0) (tt2 0) (s 0))

(setq n_ n)

(loop

(if (> n_ (* 2 n)) (return 'done) nil)

(setq s (sum_Q2 n_ n))

(print n) (prints "") (print n_) (prints "")

(setq tt2 (T2 n n_ s))

(printline tt2)

(setq n_ (add1 n_))

))

⁾

⁾