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

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

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

1. (*
2.  *  Filename: testGoldenRatio.ml
3.  *    황금률(즉, 이차방정식 x^2 - x - 1  = 0 의 양의 근)을 계산한다.
4.  *
5.  *   Execute: ocaml testGoldenRatio.ml
6.  *
7.  *   Compile: ocamlc -o testGoldenRatio.exe testGoldenRatio.ml
8.  *   Execute: testGoldenRatio
9.  *
10.  *      Date:  2013/01/27
11.  *    Author: P. Kim   [ pkim _AT_ scripts.pe.kr ]
12.  *)
14. open Printf;;
16. let printUsing() =
17.     printf "Using: TestGoldenRatio [-h|-help]\n" ;
18.     printf "This calculates the value of the golden ratio.\n" ;;
20. exception ZeroCoefficientException;;
21. exception NegativeDiscriminantException;;
23. (* 이차방정식 a x^2 + b x + c  = 0 의 근을 구한다. *)
24. let findQuadraticRoot a b c =
25.     if a = 0.0 then begin
26.         raise ZeroCoefficientException
27.     end
28.     else if (b*.b -. 4.0*.a*.c < 0.0) then begin
29.         raise NegativeDiscriminantException
30.     end;
32.     let x1 = ((-.b) +. (sqrt (b*.b -. 4.0*.a*.c))) /. (2.0 *. a) in
33.     let x2 = ((-.b) -. (sqrt (b*.b -. 4.0*.a*.c))) /. (2.0 *. a) in
34.     x1, x2  ;;  (* 리턴 겂 두 개 *)
36. (* 실행 시작 지점 *)
37. let cmdArgs = Sys.argv;;
38. if (Array.length cmdArgs > 1) && (cmdArgs.(1) = "-h" || cmdArgs.(1) = "-help") then
39. begin
40.     printUsing();
41.     exit 1
42. end ;;
44. let x1, x2  = findQuadraticRoot 1.0  (-1.0)  (-1.0) ;;
46. if x1 >= x2 then begin
47.     printf "The bigger root is %g " x1 ;
48.     printf "and the less root is %g.\n" x2
49. end
50. else begin
51.     printf "The bigger root is %g, " x2 ;
52.     printf "and the less root is %g.\n" x1
53. end ;;

 

컴파일> ocamlc -o testGoldenRatio.exe testGoldenRatio.ml

실행> testGoldenRatio
The bigger root is 1.61803 and the less root is -0.618034.