C언어 / 미적분 공식 소스

#include <stdio.h>
 
// 미분 함수
double differentiate(double (*f)(double), double x, double h) {
    return (f(x + h) - f(x)) / h;
}
 
// 적분 함수
double integrate(double (*f)(double), double a, double b, int n) {
    double integral = 0.0;
    double dx = (b - a) / n;
    
    for (int i = 0; i < n; ++i) {
        double x = a + i * dx;
        integral += f(x) * dx;
    }
    
    return integral;
}
 
// 예시 함수: f(x) = x^2
double exampleFunction(double x) {
    return x * x;
}
 
int main() {
    double x = 2.0;
    double h = 0.0001;
    int n = 100000;
    double a = 0.0;
    double b = 3.0;
 
    printf("함수 f(x) = x^2 에서 x = %.2f 에서의 미분값: %.6f\n", x, differentiate(exampleFunction, x, h));
    printf("함수 f(x) = x^2 에서 [%.2f, %.2f] 구간에서의 적분값: %.6f\n", a, b, integrate(exampleFunction, a, b, n));
 
    return 0;
}