Как интегрировать по методу Симпсона?

class proga5 {
    public static void main(String args[]) {
        int n=1000;
        
    
        task123(n,3*(Math.PI)/2.0,0.0);
        
        
    }
 
 
    public static void task123(int n, double b, double a){
        double[] gf = new double[n];
        double[] nj = new double[n];
        double  y=0;
    
        double x=Math.sin(y)*Math.pow(Math.E,y);
        gf[0]=y;
        nj[0]=x;
        
        for(int i=1;i<n;i++) {
            y=y+(b-a)/(n-1);
            x=Math.sin(y)*Math.pow(Math.E,y);
            gf[i]=y;
            nj[i]=x;
            //  System.out.println(y+" grgrwgr");
        }
        
 
            System.out.print("N        "+"l0                    "+"l1                   "+"l2");
            System.out.println();
            System.out.print(n+"     ");
        //  System.out.println(nj[j]);
            
        
        double h=(b-a)/(n-1);
        //System.out.println(h+"fewfewf");
    pryamougolnik(n, h,a,b);
    trapecia(n,h,a,b);
 
    }
    public static void pryamougolnik(int n,double h,double a,double b) {
    double u=h;
    h=0;
    double result=0;
    for(int i=1;i<n;i++) {
    h+=(b-a)/(n-1);
    double q=Math.sin(h)*Math.pow(Math.E,h);    
    result+= (u*q);
        }
    System.out.print(result+"   ");
    }
    
    public static void trapecia(int n,double h,double a,double b) {
    double u=h;
    h=0;
    double result=0;
    double result1=0;
    double a1=Math.sin(0);
    double b1=Math.sin(3*(Math.PI)/2);
    for(int i=1;i<n-1;i++) {
    h+=(b-a)/(n-1);
    double q=Math.sin(h)*Math.pow(Math.E,h);    
     result1=result1+q;
        }
         result=u*(0.5*(a1+b1)+result1);
    System.out.println(result);
    }
    
    public static void simpson(int n,double h,double a,double b){
    
    double u=h;
    h=0;
    double result=0;
    double result1=0;
    double a1=Math.sin(0);
    double b1=Math.sin(3*(Math.PI)/2);
    for(int i=1;i<n-1;i++) {
    h+=(b-a)/(n-1);
    double q=Math.sin(h)*Math.pow(Math.E,h);    
     result1=result1+q;
    }
    
}