double integrand(double x)//подынтегральная функция
{
	return x*x*x;//Есть смысл сделать эту функций inline-функцией,или это произойдёт автоматически и не стоит заморачиваться ?
}
double a = 1.0;
double b = 3.0;
int N = 1000000000;//10^9, ~71 секунда
double tau = (b - a) / N;
double res = 0;

int main(int argc,char** argv)
{
	for (double i = a; i < b; i += tau)//правда что раньше инкрементом могло быть только целое число ?
	{
		//printf("%f \n", i);//как сильно он замедляет работу
		res += integrand(i);
	}
	res = res*tau;
	printf("integral calculated by rectangles= %f \n", res);
	return 0;
}

double f(double x)//подынтегральная функция
{
	return x*x*x;//Есть смысл сделать эту функций inline-функцией,или это произойдёт автоматически и не стоит заморачиваться ?
}

double int_by_rec(double(*integrable_function)(double), double a, double b, int nunmber_of_segments)
{
	double res = 0;
	double tau = (b - a) / nunmber_of_segments;
	for (double i = a; i < b; i += tau)
	{
		res += (*integrable_function)(i);
	}
	res *= tau;
	return res;
}

int main(int argc,char** argv)
{
	printf("integral calculated by rectangles= %f \n", int_by_rec(f,1.0,3.0,1000000000));
	return 0;
}

double(*integrable_function)(double),