Текст програми. /*програма знаходження iнтеграла методом лівих прямокутникiв*/

 

/*програма знаходження 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(

Завдання. Методом трапецій знайти

 

Програмна реалізація

 

Здійснена на мові С.