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

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

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

1. #  Filename: testGoldenRatio.jl
2. #    황금률(즉, 이차방정식 x^2 - x - 1  = 0 의 양의 근)을 계산한다.
3. #
4. #   Execute: julia testGoldenRatio.jl
5. #
6. #      Date:  2013. 3. 5.
7. #    Author:  pkim __AT__ scripts ((DOT)) pe ((DOT)) kr
9. function printUsing()
10.     println("Using: python testGoldenRatio.py [-h|-help]")
11.     println("This calculates the value of the golden ratio.")
12. end
14. # 이차방정식 a x^2 + b x + c  = 0 의 근을 구한다.
15. function findQuadraticRoot(a, b, c)
16.     if a == 0.0
17.         throw("Since the highest coefficient is zero, the given equation is not a quadratic equation.")
18.     elseif b*b - 4*a*c < 0.0
19.         throw("Since the discriminant " * "$(b*b - 4*a*c)" * " is negative, the given equation has no real root.")
20.     end
22.     x1 = (-b + sqrt(b*b - 4*a*c)) / (2.0 * a)
23.     x2 = (-b - sqrt(b*b - 4*a*c)) / (2.0 * a)
24.     return (x1, x2)
25. end
28. # 실행 시작 지점
29. if length(ARGS) > 0 && (ARGS[1] == "-h" || ARGS[1] == "-help")
30.     printUsing()
31.     exit(1)
32. end
34. try
35.     # values = findQuadraticRoot(0.0, -1.0, -1.0)
36.     # values = findQuadraticRoot(1.0, -1.0, 3.0)
37.     values = findQuadraticRoot(1.0, -1.0, -1.0)
38.     x1 = values[1]
39.     x2 = values[2]
40.     if x1 >= x2
41.         println("The bigger root is " * string(x1) * ", ")
42.         println("and the less root is " * string(x2) * ".")
43.     else
44.         println("The bigger root is " * string(x2) * ", ")
45.         println("and the less root is " * string(x1) * ".")
46.     end
47. catch x
48.     println(x)
49. end

실행> julia testGoldenRatio.jl
The bigger root is 1.61803398875,
and the less root is -0.61803398875.