[C언어 ] 라이프니츠 공식(Leibniz formula)을 50만항을 계산하는 프로그램

 

 

 

안녕하세요 성조입니다.

이전에 연구실에서 풀었던 문제를 갖고 왔습니다!

소스코드, 속도 측정한 소스코드,  sqrt를 사용한 속도 측정 소스코드로 나눠서 코드에 대해 포스팅하겠습니다.

 

(개인)접근 방법

1. (1/n) -> 분모는 계속 n+2값을 갖는다.

2. 500,000번 반복한다.

3. +, -, +, -...가 항이 모두 끝날 때까지 반복한다.

 

조건으로 문제를 풀이했고 다음과 같은 소스코드를 만들었습니다.

#include <stdio.h>
 
int main()
{
    double pi = 0, num = 1;
 
    printf("\n\t");
    for (int i = 1; i <= 500000; i++)
    {
        if (i % 2 == 0)
            pi = pi - (1 / num);
        else
            pi = pi + (1 / num);
        num = num + 2;
 
    }
    pi = pi * 4;
 
    printf("%f", pi);
    printf("\n");
}

속도 측정 값

1~2 m/s이 반복해서 측정됐습니다.

 

이후 sqrt 함수를 사용하여 측정된 값과 속도를 구하면 어떨까? 란 의문이 들어 실험해봤습니다.

 

 

sqrt 함수를 사용해본 소스코드

#include <stdio.h>
#include <math.h>
#include <time.h>
 
int main()
{
    double pi = 0, num = 1;
 
    printf("\n\t");
    clock_t start, end;
    start = clock();
    for (int i = 1; i <= sqrt(500000); i++)
    {
        if (i % 2 == 0)
            pi = pi - (1 / num);
        else
            pi = pi + (1 / num);
        num = num + 2;
 
    }
    pi = pi * 4;
 
    printf("%f", pi);
    printf("\n");
    end = clock();
    printf("\t%d ms", end - start);
}