вычисление интеграла

//Lab 2_2_1
//Integral po formule tropecii
#include <STDIO.h>
#include <MATH.h>
#include <CONIO.h>
#define n 100
int i;
float a, b, h, x, s, sn, f[n];
main()
{
clrscr ();
a=0.5; b=2.5;
h=(b-a)/n;
f[1]=pow(a,2)*sqrt(pow((2+3*a),3));
f[100]=pow(b,2)*sqrt(pow((2+3*b),3));
for(i=1; i<n; i++)
         {
         x=(a+i*h);
         f[i]=pow(x,2)*sqrt(pow((2+3*x),3));
         sn+=f[i];
         }
s=(h/2)*(f[1]+f[100]+2*sn);
printf ("\nZnachenie integrala s=%f", s);
getch ();
}
 
 
 
//Lab 2_2_2
//Integral po formule Simpsona
#include <STDIO.h>
#include <MATH.h>
#include <CONIO.h>
#define n 100
int i;
float a, b, h, x, s, s1, s2, f[n];
main()
{
clrscr ();
a=0.5; b=2.5;
h=(b-a)/n;
f[1]=pow(a,2)*sqrt(pow((2+3*a),3));
f[100]=pow(b,2)*sqrt(pow((2+3*b),3));
for(i=1; i<n-1; i++)
         {
         x=(a+i*h);
         if (i%2==0)
            {
            f[i]=pow(x,2)*sqrt(pow((2+3*x),3));
            s1+=f[i];
            }
         else
            {
            f[i]=pow(x,2)*sqrt(pow((2+3*x),3));
            s2+=f[i];
            }
 
         }
 
s=(h/3)*(f[1]+f[100]+(2*s1)+(4*s2));
printf ("\nZnachenie integrala s=%f", s);
getch ();
}

#include<stdio.h>
#include<math.h>
#define n 100
const double a=2.0, b=10.0;
const double h=(b-a)/n;
double func[n+1];
double f(double x)
{
return(x+1)*2*sqrt(log10(x));
}
double trapecii()
{
double result=0.0;
for(int i=0; i<n; ++i)
result+=(func[i]+func[i+1])*h/2;
return result;
}
double simpson()
{
double result=0.0;
result+=(func[0]+func[n])*h/3;
for(int i=1; i<n; i+=2)
result +=4*h/3*func[i];
for( i=2; i<n; i+=2)
result +=2*h/3*func[i];
return result;
}
void main()
{
for(int i=0; i<=n; ++i)
func[i]=f(a+i*h);
double I1, I2;
I1=trapecii();
I2=simpson();
printf("%f\n%f\n",I1,I2);
scanf("%F%f",I1,I2);
}