metodo simpson 2

/*PRIMER DEPARTAMENTAL MÉTODOS NUMÉRICOS
MÉTODO DE SIMPSON 1/3 EJERCICIO 2*/
#include <stdio.h>
#include <stdlib.h> 
#include <conio.h> 
#include <math.h>
float linf,lsup,a,b,x1,fa,fb,fx1,Ea,Aa,An,a1,a2,b1,b2,linf1,lsup2,lsup4;
main()
{
	system("color 3f");
	gotoxy(27,9);
	printf("INSTITUTO POLIT%cCNICO NACIONAL",144);	
	gotoxy(16,11);
	printf("x");
	gotoxy(26,13);
	printf("x");
	gotoxy(33,15);
	printf("M%cTODOS NUM%cRICOS",144,144);
	gotoxy(29,17);
	printf("xxx",181,224);
	gotoxy(6,24);
	printf("PRESIONE CUALQUIER TECLA PARA CONTINUAR");	
	getch();
	system("cls");
	system("color 5f");
	gotoxy(5,5);
	printf("Funci%cn: %c(x)=(3x%c-18X+5)/5 dx",162,159,253);
	gotoxy(5,6);
	printf("Inserte el l%cmite inferior\t",161);
	scanf("%f",&linf);
	gotoxy(5,7);
	printf("Inserte el l%cmite superior\t",161);
	scanf("%f",&lsup);
	gotoxy(5,8);
	system("cls");
	system("color 4f");
	gotoxy(5,5);
	printf("Funci%cn: %c(x)=(3x%c-18X+5)/5 dx	Lim inf=%.0f   Lim sup=%.0f",162,159,253,linf,lsup);
	//MÉTODO NUMÉRICO
	a=linf;
	b=lsup;
	x1=(b-a)/2;
	gotoxy(5,7);
	printf("a=%.4f	b=%.4f	x1=%.4f",a,b,x1);
	fa=(3*(pow(a,2))-18*(a)+5)/5;
	fb=(3*(pow(b,2))-18*(b)+5)/5;
	fx1=(3*(pow(x1,2))-18*(x1)+5)/5;
	//An= Área numérica y Aa= Área analítica
	An=((x1/3)*(fa+4*fx1+fb));
	gotoxy(5,9);
	printf("f(a)=%.4f f(b)=%.4f f(x1)=%.4f",fa,fb,fx1);
	//MÉTODO ANALÍTICO
	/* Según las fórmulas analíticas, la Integral
	A=(3x²-18X+5)/5 dx Se resuelve separándola, lo que nos deja
	1/5x³-9/5x²+x y para los límites usamos la siguiente fórmula
	lim=(lsup=(1/5(lsup)³-9/5(lsup)²+(lsup)) - (linf=(1/5(linf)³-9/5(linf)²+(linf))
	*/
	a1=(pow(linf,3));
	a2=(pow(linf,2));
	b1=(pow(lsup,3));
	b2=(pow(lsup,2));
	linf1=(((1/5)*(a1))-((9/5)*(a2))+linf);
	lsup2=(0.2*b1)-(1.8*b2)+lsup;
	/*
	PRUEBA PARA COMPROBAR FUNCIONAMIENTO DE FÓRMULAS
	gotoxy(5,11);
	printf("Linf= %.4f Lsup = %.4f b1= %.4f b2= %.4f",linf1,lsup4,b1,b2);
	*/
	Aa=lsup2-linf1;
	Ea=Aa-An;
	gotoxy(5,11);
	printf("An= %.4f u%c",An,253);
	gotoxy(5,13);
	printf("Aa= %.4f u%c",Aa,253);
	gotoxy(5,16);
	printf("Error de Aproximaci%cn %fu%c",162,Ea,253);
	getch();
}