Расчет правила Симпсонов < 10^-6

Question

Я пытаюсь вычислить правило Симпсонов, чтобы получить точность < 10^-6. Но всякий раз, когда функция запускается, я получаю ошибку NaN или DNE из-за 0/sin(0). Как мне исправить ошибку, когда я =1? ** Отредактировал мою функцию, чтобы включить, когда я == 1.

public static double simpsonsRuleFunction1(double valueN, double valueA, double valueB, double valueDx) {

    double e = 0.0;
    double simpsonsRule = 0.0;
    double valueHolder = 0.0;

    valueN = 2;
    valueA = 0;
    valueB = (Math.PI)/2; 

    for(int i = 1; i<=valueN+1 ; i++){
        valueDx = (valueB-valueA)/valueN;
        e = valueA + ((i-1)*valueDx);

        if (i==1) {
        // Limit as x -> 0
        simpsonsRule += Math.pow(10,-10);

        }
        else if ((i % 2 == 0) && ( i > 1) && (i < valueN+1 ))  {
        simpsonsRule += 4*(e/((Math.sin(e))));
        }
        else if ((i % 2 != 0) && ( i > 1) && (i < valueN+1 )) {
        simpsonsRule += 2*(e/((Math.sin(e))));   
        }
        else if (i == valueN+1 ) {
            simpsonsRule += (e/((Math.sin(e))));

        }

    }
    simpsonsRule = simpsonsRule *((valueDx)/3);



     while(Math.abs(valueHolder - simpsonsRule) > Math.pow(10,-6)) {
           System.out.println("\nValueHolder" + valueHolder);
           valueHolder = simpsonsRule;
           valueN +=2;
           valueDx = (valueB-valueA)/valueN;
           simpsonsRule = 0;
    for(int i = 1; i<=valueN + 1; i++){
        e = valueA + ((i-1)*valueDx);

        if (i==1) {
        // Limit as x -> 0
        simpsonsRule += Math.pow(10,-10);

        }
        else if (i % 2 == 0) {
        simpsonsRule += 4*(e/((Math.sin(e))));
        }
        else if ((i % 2 != 0) && ( i > 1) && (i < valueN + 1)) {
        simpsonsRule += 2*(e/((Math.sin(e))));   
        }
        else if (i == valueN + 1) {
            simpsonsRule += (e/((Math.sin(e))));

        }

    }
    simpsonsRule = simpsonsRule *((valueDx)/3);


    }
    return valueN;
}

Answer 1

Вам нужно будет обработать случай, когда e==0 в явном виде.

Я бы порекомендовал вытащить функцию, чтобы вам не пришлось повторять логику:

double f(e) {
  return e==0 ? 1 : e/Math.sin(e);
}

А затем просто используйте это везде, где вам нужно вычислить функцию, например,

simpsonsRule += 4*(e/((Math.sin(e))));

становится

simpsonsRule += 4*f(e);

И, конечно же, вам также может понадобиться +/- N piв зависимости от входных значений, которые вы разрешаете.

Правка для предложения pcarter: вы можете использовать ненулевой порог, ниже которого использовать 1 в качестве значения функции, для обработки числовых неточностей около нуля, например

return (Math.abs(e) < 1e-10) ? 1 : e/Math.sin(e);

Где порог (в данном случае 1e-10) выбран, чтобы дать достаточно точный результат.