#include <stdio.h>
#include <conio.h>
#include <math.h>

int n,p;
float table[10],matrix[20][20];
float a,b,h,y;

void getfx()
{
	int i;
	printf("Enter the degree of equation : ");
	scanf("%d",&n);
	printf("Enter the coefficients in decreasing order of degree :\n");
	for(i=n;i>=0;i--) scanf("%f",&table[i]);
}

void displayfx()
{
	int i;
	printf("The entered equation is :");
	for (i=n;i>=1;i--) printf("%.2fX^%d + ",table[i],i);
	printf("%.2f",table[0]);
}

void getlimit()
{
	int i;
	printf("\nLower limit = ");
	scanf("%f",&a);
	printf("Upper limit = ");
	scanf("%f",&b);
	printf("\nEnter the number of intervals required : ");
	scanf("%d",&p);
	matrix[0][0]=a;
	h=(b-a)/p ;
	for (i=0;i<=p;i++) matrix[i+1][0]=matrix[i][0]+h;
}

float cal2fx(float j)
{
	int i;
	float sum=table[0];
	for(i=n;i>=1;i--)  sum=sum+pow(j,i)*table[i];
	return exp(sum);
}

void cal1fx()
{
	int i;
	for(i=0;i<=p;i++) matrix[i][1]=cal2fx(matrix[i][0]);
}

void integrate()
{
	int i;
	y = matrix[0][1]+matrix[p][1];
	for(i=1;i<=p-1;i++)
	{
		if (i%3 !=0) y+=3*matrix[i][1];
		else y+=2*matrix[i][1];
	}
	y*=3*h/8;
}

void main()
{
	int i;
	clrscr();
	getfx();
	displayfx();
	getlimit();
	cal1fx();
	integrate();
	printf("%.1f\n",b);
	printf("%c  e^(",244);
	for (i=n;i>=1;i--) printf("%.2fX^%d + ",table[i],i);
	printf("%.2f",table[0]);
	printf(")dx= %.4f\n",y);
	printf("%c\n",245);
	printf("%.1f",a);
	getch();
}