황금비율(golden ratio) 구하기 with Common Lisp

다음은  이차방정식 x^2 - x - 1  = 0 의 양의 근 즉 황금비율(golden ratio)을 구하는 Common Lisp 애플리케이션 소스이다. 황금비율을 구하는 비례방정식은   1 : x = x : (x+1) 이며, 이를 이차방정식으로 표현한 것이 x^2 - x - 1  = 0 이다.

See:  http://en.wikipedia.org/wiki/Golden_ratio

아래의 소스는 CLisp 으로 테스트되었다.

 

1. ;;  Filename: testGoldenRatio.lsp
2. ;;    황금률(즉, 이차방정식 x^2 - x - 1  = 0 의 양의 근)을 계산한다.
3. ;;
4. ;;   Execute: clisp testGoldenRatio.lsp
5. ;;
6. ;;      Date:  2013. 9. 1.
8. (defun printUsing()
9.     (format t "Using: lisp testGoldenRatio.lsp [-h|-help]~%")
10.     (format t "This calculates the value of the golden ratio.~%")
11. )
13. (defun handle-zero-error(msg)
14.     (catch 'handle-zero-error
15.       (foprmat t msg)
16.       (quit)))
18. ;; 이차방정식 a x^2 + b x + c  = 0 의 근을 구한다.
19. (defun findQuadraticRoot(a b c)
20.     (if (zerop a)
21.         (throw 'fn-a 'done))
22.     (if (< (- (* b b) (* 4 a c)) 0.0)
23.         (throw 'fn-a 'done))
24.     (setf x1 (/ (+ (- b) (sqrt (- (* b b) (* 4 a c)))) (* 2.0 a)))
25.     (setf x2 (/ (- (- b) (sqrt (- (* b b) (* 4 a c)))) (* 2.0 a)))
26.     (list x1 x2)
27. )
29. ;; 실행 시작 지점
30. (if (and (> (length ext:*args*) 0) (or (equal (nth 0 ext:*args*) "-h") (equal (nth 0 ext:*args*) "-help"))) (progn
31.     (printUsing)
32.     (quit) ))
34. (setf values (findQuadraticRoot 1.0 (- 1.0) (- 1.0)))
35. (setf x1 (nth 0 values))
36. (setf x2 (nth 1 values))
37. (if (>= x1 x2) (progn
38.     (format t "The bigger root is ~F, ~%" x1)
39.     (format t "and the less root is ~F.~%" x2) )
40.        (progn
41.     (format t "The bigger root is ~F, ~%" x2)
42.     (format t "and the less root is ~F.~%" x1) )
43. )

실행> clisp testGoldenRatio.lsp
The bigger root is 1.618034,
and the less root is -0.618034.