#include <iostream>
#define a 0
#define b 1
using namespace std;

double f(double x)
{
	return (0.5-x);
}

int main()
{
	int i;
	double max,min,area,n,temp;
	area=0.0;
	max=f(a);
	min=f(a);
	n=20.0;
	for (i=0;i <= n-1;i++)
	{
		temp=(f(a+(i/n)*(b-a))+f(a+((i+1)/n)*(b-a)))*0.5*(1/n)*(b-a);
		area+=temp;
		if (f(a+(i/n)*(b-a))>max)
			max=(f(a+(i/n)*(b-a)));
		if (f(a+((i+1)/n)*(b-a))<min)
			min=(f(a+((i+1)/n)*(b-a)));
		//cout<<max <<"   -  "<<min<<"  -  "<<area<<"  -  "<<temp<<endl;
		cout<<f(a+(i/n)*(b-a))<<"     "<<f(a+((i+1)/n)*(b-a))<<"    "<<0.5*(1/n)*(b-a)<<"   "<<temp<<endl;
	}
	cout<<"The integral of the given function f over (0,1) is "<<area<<endl;
	cout<<"Its maximum value on this interval is "<<max <<endl;
	cout<<"Its minimum value on this interval is "<<min <<endl;
	return 0;
}