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

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

  1. /*
  2.  *  Filename: testGoldenRatio.scala
  3.  *    황금률(즉, 이차방정식 x^2 - x - 1  = 0 의 양의 근)을 계산한다. 
  4.  *
  5.  *   Execute: scala testGoldenRatio.scala
  6.  *
  7.  *      Date:  2009/03/09
  8.  *    Author:  PH Kim   [ pkim (AT) scripts.pe.kr ]
  9.  */
  10. def printUsing() {
  11.     println("Using: scala testGoldenRatio.scala [-h|-help]")
  12.     println("This calculates the value of the golden ratio.")
  13. }
  14. // 이차방정식 a x^2 + b x + c  = 0 의 근을 구한다.
  15. def findQuadraticRoot(a: Double, b: Double, c: Double) : Array[Double] = {
  16.     if (a == 0.0) {
  17.         throw new RuntimeException("Since the highest coefficient is zero, the given equation is not a quadratic equation.")
  18.     }
  19.     else if (b*b - 4*a*c < 0.0) {
  20.         throw new RuntimeException("Since the discriminant " + (b*b - 4*a*c) + " is negative, the given equation has no real root.")
  21.     }
  22.     var x1 : Double = (-b + Math.sqrt(b*b - 4*a*c)) / (2.0 * a)
  23.     var x2 : Double = (-b - Math.sqrt(b*b - 4*a*c)) / (2.0 * a)
  24.     var array : Array[Double] = Array(x1, x2)
  25.     return array
  26. }
  27. // 실행 시작 지점
  28. if (args.length > 0 && (args(0).equals("-h") || args(0).equals("-help"))) {
  29.     printUsing()
  30.     System.exit(1)
  31. }
  32. var values : Array[Double] = findQuadraticRoot(1.0, -1.0, -1.0)
  33. var x1 : Double = values(0)
  34. var x2 : Double = values(1)
  35. if (x1 >= x2) {
  36.     println("The bigger root is " + x1 + ", ")
  37.     println("and the less root is " + x2 + ".")
  38. }
  39. else {
  40.     println("The bigger root is " + x2 + ", ")
  41.     println("and the less root is " + x1 + ".")
  42. }



실행> scala testGoldenRatio.scala
The bigger root is 1.618033988749895,
and the less root is -0.6180339887498949.



크리에이티브 커먼즈 라이선스
Creative Commons License


 

Posted by Scripter
,