Como calcular integral em Java usando a Regra dos Trapézios - Java para Métodos Numéricos

package estudos;

public class Estudos {
  public static void main(String[] args) {
    // vamos definir os limites de integração
    double inferior = 0;
    double superior = 2;
    
    // quantidade de trapézios. Quanto 
    // mais trapézios, maior a precisão do cálculo
    int quant_trapezios = 1000;
    
    // vamos mostrar o resultado
    System.out.println("Valor da integral com " + quant_trapezios +
      " trapézios: " + calcularIntegralTrapezio(inferior,
      superior, quant_trapezios));
  }
  
  // função que recebe o valor de x e retorna o
  // y corresponente
  static double funcao(double x) {
    return 4.0 * Math.pow(x, 3.0);
  }
  
  // função usada para calcular a integral pela soma das
  // áreas dos trapézios
  static double calcularIntegralTrapezio(double a, 
    double b, double n) {
    
    // vamos calcular a distância entre os dois intervalos
    double h = (b - a) / (n - 1);
    
    // calculamos a soma inicial usando função(a) e função(b)
    double soma = 0.5 * (funcao(a) + funcao(b));
    
    // e agora somamos os termos intermediários
    for(int i = 1; i < n - 1; i++) {
      soma = soma + funcao(a + (i * h));
    }
    
    // e retornamos o resultado da integração
    return (h * soma);
  }
}