Обчислити наближено визначений інтеграл за формулою Сімпсона, при n=10

?ecos(x)dx

Результати розв’язку в середовищі Turbo Pascal

program simpson;

uses crt;

const a=0;

b=1;

var h,x,r,yn,yp:real; n,i:integer;

function y (x1:real):real;

begin

y:= exp(cos(x));

end;

begin

clrscr;

n:=10;

yn:=0;

h:=(b-a)/n;

x:=x+h;

for i:=1 to n do

begin

if (i mod 2)=0 then

yn:=yn+cos(x)*y(x);

x:=x+h;

end;

r:=(h/3)*(y(a)+yn+yp+y(b));

writeln ('y(x)=',r:3:4);

readln;

end.

Результати розв’язку в середовищі Mathcad