c - 使用中点法积分

Question

我需要编写一个 C 程序，使用中点法计算三次多项式的积分。

基本上积分:f(x) = mx^3 + nx^2 + px + q使用 block /中点方法。

到目前为止我有这个:

#include <stdio.h>
#include <math.h>
int main(){
    char option;
    float m, n, p, q, a, b, N, S=0,T,h;//S for midpoint, T for Trap, h for interval length
    do{
        printf("\n Select one of following options:\n");
        printf("a)Enter new coefficients (m,n,p, and q)\n");
        printf("b)Enter the value for the range from a to b, and N, the number of intervals\n");
        printf("c)Calculate the integral from a to b\n");
        printf("d)Exit\n");

        option=getc(stdin);
        printf("You have selected option: \"%c\" \n", option);

        switch (option){
        case 'a':
            printf("Enter m,n,p,q (in this order):\n");
            scanf("%f%f%f%f", &m,&n,&p,&q);
            break;
        case 'b':
            printf("Enter the interval from a to b and the number of intervals N:");
            scanf("%f%f%f", &a, &b, &N);
            break;
        case 'c':
            h = (b - a) / N;
            for (float i = a; i < b - h; i + h){
                S =+(m*pow(i, 3)) + n*pow(i, 2) + p*i + q;
            }
            S = S*h;
            printf("The integral using midpoint rule: %f", S);
            break;
        case 'd':
            printf("Exit.\n");
            break;
        default:
            printf("Invalid entry, try again... \n");
            break;
        }
    } while (option != 'd');
    return 0;
}

这是我尝试使用的主要算法:

h = (b - a) / N;
            for (float i = a; i < b - h; i + h){
                S =+(m*pow(i, 3)) + n*pow(i, 2) + p*i + q;
            }
            S = S*h;
            printf("The integral using midpoint rule: %f", S);
            break;

但它不起作用。它卡在循环部分。因为在我输入所有其他变量并选择 c 选项后，它什么也没做。

该算法首先查找间隔的长度。然后循环根据中点公式添加函数，这就是我认为我做错了什么的地方。

Answer 1

通常，将浮点变量用于循环计数器被认为是一种不好的做法，因为比较浮点值是否相等很困难。解决方法是使用 int 作为计数器并将其乘以浮点增量:

double delta = (b - a) / N;
double start = a + delta / 2.0;

for (int i = 0; i < N; i++) {
    double x = start + delta * i;

    // do calculations with x
}