Значения функции записать в массив и вычислить значение интеграла

#include <STDIO.h>
#include <iostream>
#include <MATH.h>
#include <CONIO.h>
#define n 100
using namespace std;
int i;
float a, b, h, x, s, s1, s2, sn, f[n];
int main()
{
    a= 1; b= 20;
    h=(b-a)/n;
    f[1]=(a*a*log(a))/pow((1+a),2);
    f[100]=(b*b*log(b))/pow((1+b),2);
    for(i=1; i<n; i++)
    {
        x=(a+i*h);
        if(i%2==0)
        {
        f[i]=(x*x*log(x))/pow((1+x),2);
        s1+=f[i];   
        }
        else
        {
        f[i]=(x*x*log(x))/pow((1+x),2);
        s2+=f[i];   
        }
    
    }
    
    s=(h/3)*(f[1]+f[100]+(2*s1)+(4*s2));
    cout<<"Zhach. integrala (Simpson):"<<s<<endl;
    
    for(i=1; i<n; i++)
    {
        x=(a+i*h);
        f[i]=(x*x*log(x))/pow((1+x),2);
        sn+=f[i];
    }
    s=(h/2)*(f[1]+f[100]+2*sn);
    cout<<"Zhach. integrala (Trapecia):"<<s;
    return 0;
}

#include <iostream>
#include <cstdlib>
#include <math.h>
using namespace std;
void Input(float* M,int k);  
void Print(float* M,int k);   
int Work(float* M,int k);
 
 
int main()
{
    int k=101;
    cout<<"Razmer massiva: "<<endl;
 
    float* A=new float[k];
    Print(A,k); 
    int s=Work(A,k);
    delete[] A;
    return 0;
}
void Input(float* M,int k)
{
    for(int i=0;i<k;i++)
    {
            float x=rand()%100;
            float y=(x*x*log(x))/((1.+x)*(1.+x));
            M[i]=y;
    }
}
 
void Print(float* M,int k)
{
    for(int i=0;i<k;i++)
    {
        cout << M[i]<< " ";
    }
    cout <<", "<<endl;
}
int Work(float* M,int k)
{int i;
float h, s, s1, s2, sn; 
char f[1];
double a,b,x;
    a= 1; b= 20;
    h=(b-a)/k;
    f[1]=(a*a*log(a))/((1.+a)*(1.+a));
    f[100]=(b*b*log(b))/((1.+b)*(1.+b));
    for(i=1; i<k; i++)
    {
        x=(a+i*h);
        if(i%2==0)
        {
        f[i]=(x*x*log(x))/((1.+x)*(1.+x));
        s1+=f[i];   
        }
        else
        {
        f[i]=(x*x*log(x))/((1.+x)*(1.+x));
        s2+=f[i];   
        }
    
    }
    
    s=(h/3)*(f[1]+f[100]+(2*s1)+(4*s2));
    cout<<"Zhach. integrala (Simpson):"<<s<<endl;
    
    for(i=1; i<k; i++)
    {
        x=(a+i*h);
        f[i]=(x*x*log(x))/((1.+x)*(1.+x));
        sn+=f[i];
    }
    s=(h/2)*(f[1]+f[100]+2*sn);
    cout<<"Zhach. integrala (Trapecia):"<<s;
    return 0;
}

float h, s, s1 = 0, s2 = 0, sn = 0;

#include <iostream>
#include <cstdlib>
#include <math.h>
using namespace std;
void Input(double* M,int k);  
void Print(double* M,int k);   
int Work(double* M,int k);
 
 
int main()
{
    int k=100;
    cout<<"Razmer massiva: "<<endl;
 
    double* A=new double[k];
    Print(A,k); 
    int s=Work(A,k);
    delete[] A;
    return 0;
}
void Input(double* M,int k)
{
    for(int i=0;i<k;i++)
    {
            double x=rand()%100;
            double y=(x*x*log(x))/((1.+x)*(1.+x));
            M[i]=y;
    }
}
 
void Print(double* M,int k)
{
    for(int i=0;i<k;i++)
    {
        cout.precision(2);
        cout<< M[i]<<", ";
    }
    cout <<"\n"<<endl;
}
int Work(double* M,int k)
{
int i;
double h, s, s1 = 0, s2 = 0, sn = 0;
char f[99];
double a,b,x;
    a= 1; b= 20;
    h=(b-a)/k;
    f[0]=(a*a*log(a))/((1.+a)*(1.+a));
    f[99]=(b*b*log(b))/((1.+b)*(1.+b));
    for(i=1; i<k; i++)
    {
        x=(a+i*h);
        if(i%2==0)
        {
        f[i]=(x*x*log(x))/((1.+x)*(1.+x));
        s1+=f[i];   
        }
        else
        {
        f[i]=(x*x*log(x))/((1.+x)*(1.+x));
        s2+=f[i];   
        }
    
    }
    
    s=(h/3)*(f[0]+f[99]+(2*s1)+(4*s2));
    cout<<"Zhach. integrala (Simpson):"<<s<<endl;
    
    for(i=1; i<k; i++)
    {
        x=(a+i*h);
        f[i]=(x*x*log(x))/((1.+x)*(1.+x));
        sn+=f[i];
    }
    s=(h/2)*(f[0]+f[99]+2*sn);
    cout<<"Zhach. integrala (Trapecia):"<<s;
    return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
 
double F(double x)
{
    return x*x*log(x)/((1.+x)*(1.+x));
}
 
double Trapezium(int n, double a, double b, double h) 
{
    double sum=F(a)/2.;
    for(int i = 1; i < n; i++)    
        sum+=F(a+i*h);
        sum+=F(b)/2.;
    return h*sum;    
}
 
double Simpson(int n, double a, double b, double h)
{
    double s1=0., s2=0.;
    
    for(int i = 1; i < n; i++)
    {
        if(i%2==0) s1+=F(a+i*h);        
        else s2+=F(a+i*h);    
    }    
    return h/3.*(F(a)+F(b)+2.*s1+4.*s2);
}  
 
int main()
{
    double a=1., b=20., n=100, h=(b-a)/n;
    
    cout << Trapezium(n,a,b,h) << "  " << Simpson(n,a,b,h) << "\n";
    
system("pause");
return 0;
}

#include <iostream>
#include <cmath>
using namespace std;
 
double F(double x)
{
    return x*x*log(x)/((1.+x)*(1.+x));
}
 
double Trapezium(double *T, int n, double a, double b, double h) 
{    
    double sum=F(a)/2.;
    T[0]=F(a);
    for(int i = 1; i < n; i++)
        {
        sum+=F(a+i*h);
        T[i]=F(a+i*h);
        }
        sum+=F(b)/2.;
        T[100]=F(b);
    return h*sum;    
}
 
double Simpson(double *S, int n, double a, double b, double h)
{
    double s1=0., s2=0.;
    S[0]=F(a);
    for(int i = 1; i < n; i++)
    {
        if(i%2==0) s1+=F(a+i*h);        
        else s2+=F(a+i*h); 
        S[i]=F(a+i*h);
    } 
    S[100]=F(b);
    return h/3.*(F(a)+F(b)+2.*s1+4.*s2);
}  
 
int main()
{
    int n=100;
    double a=1., b=20., h=(b-a)/n;
    double T[101], S[101];
    
    cout << Trapezium(T,n,a,b,h) << " " << Simpson(S,n,a,b,h) << "\n\n";
    
    for(int i = 0; i <= n; i++)
      cout << T[i] << "   " << S[i] << "\n";
      cout << "\n";      
    
system("pause");
return 0;
}