Интегрирование функции методами симпсона(порабол) и трапеций

    Приложение  А

    Листинг файла integrals.h

    #include<stdio.h>

    #include<conio.h>

    #include<math.h>

    double rectangle(double n,double a,double b)

    {

    double y=0,x=a,c=(b-a)/n;

    for(int i=0;i<=n;i++)

    {

    y+=pow((1+3*x+2*pow(x,2)),(-0.5));

    x+=c;

    }

    y*=c;

    return y;

    }

    double simpson(double n,double a,double b)

    {

    double y=0,x=a,c=(b-a)/n;

    for(int i=0;i<=n;i++)

    {

    if ((i==0)||(i==n))

    y+=pow((1+3*x+2*pow(x,2)),(-0.5));

    else

    if(((i%2)==0)&&(i!=n))

    y+=2*(pow((1+3*x+2*pow(x,2)),(-0.5)));

    else

    y+=4*(pow((1+3*x+2*pow(x,2)),(-0.5)));

    x+=c;

    }

    y*=c/3;

    return y;

    }

    double proverka(double a,double b)

    {

    double y;

    y=pow(2,(-0.5))*(log(4*b+3+pow((pow((4*b+3),2)-1),0.5))-log(4*a+3+pow((pow((4*a+3),2)-1),0.5)));

    return y;

    }