Viết hàm tính tích phân của f(x)=x^n trên một đoạn [a, b] với n, a và b là tham số truyền vào.

float nguyenham(int x, int n)
{
    n = n + 1;
    float s = x;
    for (int i = 1; i < n; i++)
        s *= x;
    s = (float)s / n;
    return s;
}

float tichphan(int a, int b, int n)
{
    float s;
    s = nguyenham(a, n) - nguyenham(b, n);
    return s;
}