Разработать программу для вычисления определенного интеграла методом средних прямоугольников. Проверка курсовой

#include "stdafx.h"
#include <math.h>
#include <iostream>
 
using namespace std;
 
//---------------------------------------------------------------------------
//Подынтегральная функция
double func(double,double,double,double);
 
//---------------------------------------------------------------------------
// Метод средних прямоугольников
double Int(double,double,double);
 
//---------------------------------------------------------------------------
// Проверка количество точек разбиения отрезка интегрирования (n) 
int test(int n);
 
void main()
{
setlocale(LC_ALL,"Russian" );
 
double x1,x2,a,b,c;
unsigned long int n;
cout<<"Задание:"<<endl<<"Разработать программу для вычисления определенного интеграла"<<endl<<"методом средних прямоугольников для функции типа:"<<endl<<"y=a*cos(b*x)*x^c"<<endl<<"Коэффициенты a, b, c и границы интервала интегрирования считываются с экрана."<<endl;
 
cout<<"-----------------------------------------------------------------------------"<<endl;
cout<<"\nОпределите границы интервала:"<<endl<<"\nx1:\t"; 
cin>>x1;
cout<<"x2:\t";
cin>>x2;
 
 
while(true)
{
int cs;
cout<<"\nВведите количество точек разбиения отрезка интегрирования n."<<endl<<"\nn:\t";
cin>>n;
cs=test(n);
if (cs==0)
cout<<"\nЗначение n меньше единицы. Введите заново.\n";
else
break;
}
 
//--------------------------------------------------------------------------- 
//Вывод результата
 
cout << endl <<"\nРезультат интегрирования= " << Int(x1,x2,n)<<endl;
 
}
 
//---------------------------------------------------------------------------
// Метод средних прямоугольников
double Int(double x1,double x2,double n)
{
double a,b,c;
cout<<"\nВведите коэффициенты подынтегральной функции{y=a * cos (b * x )+ x^c}:"<<endl<<"\na\t";
cin>>a;
cout<<"b:\t";
cin>>b;
cout<<"c:\t";
cin>>c;
    double dx=fabs(x1-x2)/n;//приращение
    double Integ=0;
    for ( int i=0 ; i<n ; i++)
    {
        Integ += func((x1+dx*i),a,b,c);
    }
    Integ *= dx;
    return Integ ;
}
//---------------------------------------------------------------------------
//Подынтегральная функция
 
double func(double x,double a,double b,double c)
{
return (a * cos (b * x )+ pow (x,c));
}
 
//---------------------------------------------------------------------------
// проверка количество точек разбиения отрезка интегрирования (n) 
int test(int n)
{
if (n<1)
return 0; 
else
return 1;
}