.


           


Procedure Simpson(a, b : real; n : integer; var j : real);

var

dx, c, c1, f : real;

i : integer;

begin

dx := (b - a)/n;

c := a;

c1 := a + dx/2;

f := fx(a) + fx(b) + 4*fx(c1);

for i := 1 to n - 1 do

begin

c := c + dx; c1 := c1 + dx;

f := f + 2*fx(c) + 4*fx(c1)

end;

j := (dx/6)* f

end;

54. -.

Function I(n : longint; a, b : real) : real;

var

x, f : real; k : longint;

begin

randomize;

f := 0;

for k := 1 to n do

begin

x := random*(b - a) + a;

f := f + fx(x)

end;

I := (b - a)*f/n

end;

{----------------------------------------------------------------------------------------}

Procedure Monte_Karlo(a, b, eps : real; var s : real);

var

t, d, f, dd : real; n : longint;

begin

n := 1;

t := I(n, a, b);

s := t;

d :=0;

repeat

n := n + 1;

t := t + I(n, a, b);

s := t/n;

d := d + (n/(n - 1))*sqr(I(n, a, b) - s);

dd := d/(n - 1)

until 5*sqrt(dd/n) < eps

end;



184.
Pi/6 61Pi/360.

185.
dy x = 1 dx = 0.2.

186. , y = ex sinx
y = e-x sinx
( x0 x ).

187. , :

x3 - 9x + 2 = 0, xe = 2, x = 0.538sinx + 1, a = ax a > 1, x2 arctgx = a,


188.

,


189. , ,
.
.