Завдання:
Обчислити площу криволінійної трапеції методом прямокутників, трапецій і Сімпсона.
X
0.5
2.0
3.5
5.0
6.5
8.0
9.5
11.0
12.5
14.0

y
1.3
-1.3
-2.5
-1.5
0.4
1.8
1.9
2.2
2.1
4.1


#include <stdio.h>
#include <conio.h>
#include <math.h>
#define N 10
main()
{float x[N]={0.5 2.0 3.5 5.0 6.5 8.0 9.5 11.0 12.5 14.0},
y[N]={1.3 -1.3 -2.5 -1.5 0.4 1.8 1.9 2.2 2.1 4.1 },
h=0.25,sL=0,sR=0,sLR,sTr=0,sSim=0;
short i;
clrscr();
for(i=0;i<=N-1;i++)
y[i]=y[i]+25;
getch();
for(i=0;i<=N-1;i++)
printf(" %.2f ",x[i]);
printf("\n");
for(i=0;i<=N-1;i++)
printf(" %.2f ",y[i]);
for(i=0;i<=N-2;i++)
sL=sL+h*y[i];
printf("\n sL=%.2f ",sL);
for(i=1;i<=N-1;i++)
sR=sR+h*y[i];
printf("\n sR=%.2f ",sR);
sLR=(sL+sR)/2;
printf("\n sLR=%.3f ",sLR);
for(i=1;i<=N-2;i++)
sTr=sTr+h*y[i];
sTr=sTr+h*(y[0]+y[N-1])/2;
printf("\n sTr=%.3f ",sTr);
for(i=1;i<=N-2;i+=2)
sSim=sSim+4*h*y[i]/3;
for(i=2;i<=N-2;i+=2)
sSim=sSim+2*h*y[i]/3;
sSim=sSim+h*(y[0]+y[N-1])/3;
printf("\n sSim=%.4f ",sSim);
getchar();
return 0;
}
sL = 11.856
sR = 12.899
sLR = 12.342
sTr = 13.657
sSim=12.949