Метод Симпсона

#include<iostream>
#include<math.h>
using namespace std;
void metodSympsona(double a,double b)
{
        double x[41],y[41],iks[41],igrik[41],m,h,integral,s,sum=0;int i,j,k;
    cout<<"vvedite kolvo x m";
    cin>>k;
        switch(k)
    {
    case 1:
        m=10;
        break;
    case 2:
        m=20;
        break;
    case 3:
        m=40;
        break;
    default:
        cout<<"Vibrano nedopustimoe zna4enie";
        }
        h=(b-a)/m;
        for(i=1;i<m;i++)
        {
            x[i]=a+(i-1)*h;
            y[i]=0.1*pow(x[i],3)+pow(x[i],2)-10*sin(x[i]);
            iks[i]=(x[i]+x[i+1])/2;
            igrik[i]=0.1*pow(iks[i],3)+pow(iks[i],2)-10*sin(iks[i]);
        }
        for(i=1;i<m;i++)
        {
            s=y[i]+4*igrik[i]+y[i+1];
                sum=sum+s;
        }
        integral=h/6*s;
        cout<<integral;
}
int main (void)
{
    double a=-4,b=2;
    metodSympsona(a,b);
    system("pause");
    return(0);
}

float function(float x){
    return sin(x);
}
 
float simpson(float left_gr, float right_gr, int n){
    float ksi = 0.0;
    float delta_x = 0.0;
    float simpson_integral = 0.0;
    int m = 2 * n;
 
    delta_x = (right_gr - left_gr) / m;
    simpson_integral = function(left_gr) + function(right_gr);
    ksi = left_gr;
 
    for (int i = 1; i < m; ++i){
        ksi += delta_x;
        if ((i % 2) == 0){
            simpson_integral += 2 * function(ksi);
        }
        else{
            simpson_integral += 4 * function(ksi);
        }
    }
 
    return simpson_integral * delta_x / 3;
}