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

다음은  이차방정식 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: testGoldenRatio.c
3.  *    황금률(즉, 이차방정식 x^2 - x - 1  = 0 의 양의 근)을 계산한다.
4.  *
5.  *   Compile: cl testGoldenRatio.c
6.  *
7.  *   Execute: testGoldenRatio
8.  *
9.  *      Date:  2008/03/24
10.  *    Author:  PH Kim   [ pkim (AT) scripts.pe.kr ]
11.  */
13. #include <stdio.h>
14. #include <string.h>
15. #include <math.h>
17. typedef struct _PAIR {
18.     double x1;
19.     double x2;
20. } PAIR;
22. void printUsing() {
23.     printf("Using: testGoldenRatio [-h|-help]\n");
24.     printf("This calculates the value of the golden ratio.\n");
25. }
27. // 이차방정식 a x^2 + b x + c  = 0 의 근을 구한다.
28. PAIR *findQuadraticRoot(double a, double b, double c) {
29.     static PAIR zeros;
30.     if (a == 0.0) {
31.         fprintf(stderr, "Since the highest coefficient is zero, the given equation is not a quadratic equation.\n");
32.         exit(1);
33.     }
34.     else if (b*b - 4*a*c < 0.0) {
35.         fprintf(stderr, "Since the discriminant %f is negative, the given equation has no real root.\b", b*b - 4*a*c);
36.     exit(1);
37.     }
39.     zeros.x1 = (-b + sqrt(b*b - 4*a*c)) / (2.0 * a);
40.     zeros.x2 = (-b - sqrt(b*b - 4*a*c)) / (2.0 * a);
41.     return (PAIR *) &zeros;
42. }
44. void main(int argc, char *argv[]) {
45.     double x1, x2;
46.     PAIR *values = (PAIR *) findQuadraticRoot(1.0, -1.0, -1.0);
47.     if (argc > 1 && (strcmp(argv[1], "-h") == 0 || strcmp(argv[1], "-help") == 0)) {
            printUsing();
48.         exit(1);
49.     }
51.     values = findQuadraticRoot(1.0, -1.0, -1.0);
52.     x1 = values->x1;
53.     x2 = values->x2;
54.     if (x1 >= x2) {
55.         printf("The bigger root is %lf, \n", x1);
56.         printf("and the less root is %lf.\n", x2);
57.     }
58.     else {
59.         printf("The bigger root is %lf, \n", x2);
60.         printf("and the less root is %lf.\n", x1);
61.     }
62. }

컴파일> cl testGoldenRatio.c

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


Ch를 이용하면 C 언어 소스 코드를 (컴파일하지 않고) 직접 실행시킬 수 있다.

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