int main()
{
printf("\tSimpson 16\tSimpson 32\tPonto Medio 16\tPonto Medio 32\n");
printf("F1\t%10f\t%10f\t%10f\t%10f\n", simpson(0,4, 16,f1), simpson(0,4, 32,f1),ponto_medio(0,4,16,f1), ponto_medio(0,4,32,f1) );
printf("F2\t%10f\t%10f\t%10f\t%10f\n", simpson(1,3, 16,f2), simpson(1,3, 32,f2),ponto_medio(1,3,16,f2), ponto_medio(1,3,32,f2) );
printf("F3\t%10f\t%10f\t%10f\t%10f\n", simpson(0,M_PI, 16,f3), simpson(0,M_PI, 32,f3),ponto_medio(0,M_PI,16,f3), ponto_medio(0,M_PI,32,f3) );
}
int main (void)
{
// primeira funcao ------------------------------------------------------
printf("\n\tRegra de Simpson\n\t f(x) = int(xdx/sqrt(x^2+9))\n ");
printf("\tSm = 16\n\n \t S(fx) = %g\n", simpson(0, 4, fx1,16));
printf("\tSm = 32\n\n \t S(fx) = %g\n", simpson(0, 4, fx1,32));
printf("\n\tRegra do Ponto Medio\n\t f(x) = int(xdx/sqrt(x^2+9))\n ");
printf("\tMm = 16\n\n \t M(fx) = %g\n", ponto_medio(0, 4, fx1,16));
printf("\tMm = 32\n\n \t M(fx) = %g\n", ponto_medio(0, 4, fx1,32));
printf("\n\n\n\n------------------------------------------------------------------------\n");
//------------------------------------------------------------------------
// segunda funcao --------------------------------------------------------
printf("\n\tRegra de Simpson\n\t f(x) = int(x^2*log(x)dx)\n ");
printf("\tSm = 16\n\n \t S(fx) = %g\n", simpson(1, 3, fx2,16));
printf("\tSm = 32\n\n \t S(fx) = %g\n", simpson(1, 3, fx2,64));
printf("\n\tRegra do Ponto Medio\n\t f(x) = int(x^2*log(x)dx)\n ");
printf("\tMm = 16\n\n \t M(fx) = %g\n", ponto_medio(1, 3, fx2,16));
printf("\tMm = 32\n\n \t M(fx) = %g\n", ponto_medio(1, 3, fx2,32));
printf("\n\n\n\n------------------------------------------------------------------------\n");
//------------------------------------------------------------------------
// terceira funcao -------------------------------------------------------
printf("\n\tRegra de Simpson\n\t f(x) = int(x^2*sin(x)dx)\n ");
printf("\tSm = 16\n\n \t S(fx) = %g\n", simpson(0, M_PI, fx3,16));
printf("\tSm = 32\n\n \t S(fx) = %g\n", simpson(0, M_PI, fx3,32));
printf("\n\tRegra do Ponto Medio\n\t f(x) = int(x^2*sin(x)dx)\n ");
printf("\tMm = 16\n\n \t M(fx) = %g\n", ponto_medio(0, M_PI, fx3,16));
printf("\tMm = 32\n\n \t M(fx) = %g\n", ponto_medio(0, M_PI, fx3,32));
printf("\n\n\n\n------------------------------------------------------------------------\n");
//------------------------------------------------------------------------
system("PAUSE");
return 0;
}
