/*програма знаходження iнтеграла методом лівих прямокутникiв*/
#include "stdio.h"
#include "math.h"
FILE *stream;
float f(float x)
{return(exp(1/x+sin(x)));};
main ()
{
/*опис змiнних*/
float a,b,h,s=0,x;
int n,i;
stream=fopen("respr.txt","w");
printf("Введiть a ");
scanf("%f",&a);
printf("введiть b ");
scanf("%f",&b);
fprintf(stream,"Iнтегрування методом прямокутникiв\n");
fprintf(stream,"Функцiя f(x)=exp(1/x+sin(x))\n");
fprintf(stream,"Межi a=%4.4f b=%4.4f\n",a,b);
n=10000;
h=(b-a)/n;
/*цикл по вузлах*/
for (i=0;i<=n;i=i++)
{
x=a+i*h;
s=s+f(x);
};
s=s*h;
/*вивiд результату*/
fprintf(stream,"Визначений iнтеграл=%4.3f\n",s);
fclose(stream);
}
Результат роботи програми
Iнтегрування методом прямокутникiв
Функцiя f(x)=exp(1/x+sin(x))
Межi a=5.0000 b=15.0000
Визначений iнтеграл=15.473
Метод трапецій
Замінивши прямокутники розбиття з попереднього розділу трапеціями, одержимо формулу трапецій:
» h( 
Завдання. Методом трапецій знайти 
Програмна реалізація
Здійснена на мові С.
