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

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

See:  Golden ratio - Sajun.org

1. /*
2.  *  Filename: testGoldenRatioCPP.cpp
3.  *    황금률(즉, 이차방정식 x^2 - x - 1  = 0 의 양의 근)을 계산한다.
4.  *
5.  *   Compile: cl -GX testGoldenRatioCPP.cpp
6.  *
7.  *   Execute: testGoldenRatioCPP
8.  *
9.  *      Date:  2008/03/24
10.  *    Author:  PH Kim   [ pkim (AT) scripts.pe.kr ]
11.  */
13. #include <iostream>
14. #include <string>
15. #include <cmath>
17. using namespace std;
19. typedef struct _PAIR {
20.     double x1;
21.     double x2;
22. } PAIR;
24. void printUsing() {
25.     cout << "Using: testGoldenRatioCPP [-h|-help]" << endl;
26.     cout << "This calculates the value of the golden ratio." << endl;
27. }
29. // 이차방정식 a x^2 + b x + c  = 0 의 근을 구한다.
30. PAIR *findQuadraticRoot(double a, double b, double c) {
31.     static PAIR zeros;
32.     if (a == 0.0) {
33.         cerr << "Since the highest coefficient is zero, the given equation is not a quadratic equation." << endl;
34.         exit(1);
35.     }
36.     else if (b*b - 4*a*c < 0.0) {
37.         cerr << "Since the discriminant " << (b*b - 4*a*c) << " is negative, the given equation has no real root." << endl;
38.         exit(1);
39.     }
41.     zeros.x1 = (-b + sqrt(b*b - 4*a*c)) / (2.0 * a);
42.     zeros.x2 = (-b - sqrt(b*b - 4*a*c)) / (2.0 * a);
43.     return (PAIR *) &zeros;
44. }
46. void main(int argc, char *argv[]) {
47.     double x1, x2;
48.     PAIR *values = (PAIR *) findQuadraticRoot(1.0, -1.0, -1.0);
49.     if (argc > 1 && (strcmp(argv[1], "-h") == 0 || strcmp(argv[1], "-help") == 0)) {
50.         printUsing();
51.         exit(1);
52.     }
54.     values = findQuadraticRoot(1.0, -1.0, -1.0);
55.     x1 = values->x1;
56.     x2 = values->x2;
57.     if (x1 >= x2) {
58.         cout << "The bigger root is " << x1 << endl;
59.         cout << "and the less root is " << x2 << endl;
60.     }
61.     else {
62.         cout << "The bigger root is " << x2 << endl;
63.         cout << "and the less root is " << x1 << endl;
64.     }
65. }

컴파일> cl -GX testGoldenRatioCPP.cpp

실행> testGoldenRatioCPP
The bigger root is 1.618034,
and the less root is -0.618034.




이 저작물은 크리에이티브 커먼즈 코리아 저작자표시-비영리-변경금지 2.0 대한민국 라이센스에 따라 이용하실 수 있습니다.