24.04.2009, 23:42
общий
это ответ
Здравствуйте, Лилечка Васина!
Код программы в приложении.
Удачи!
Приложение:
program z_165744;
const
a = 0.4;{нижний предел}
b = 1.2;{верхний предел}
n = 100;{точность}
function f(x: real): real;
begin
f := 1/(x*sqrt(x)+3);
end;
function IntegralTrap(a, b: real): real;
{метод трапеций}
var
h, k, s: real;
i: Integer;
begin
h := (b - a)/N;
s := (f(a) + f(b))/2;
k := a;
for i := 1 to N - 1 do
begin
k := k + h;
s := s + f(k);
end;
s := s*h;
IntegralTrap := s;
end;
function IntegralSimpson(a, b: real): real;
{метод симпсона}
var
h, k, s: real;
i: Integer;
begin
h := (b - a)/N;
s := 0;
k := a;
for i := 1 to N - 1 do
begin
k := k + h;
if i mod 2 = 0 then
s := s + 2*f(k)
else
s := s + 4*f(k);
end;
s := (s + f(a) + f(b))*h/3;
IntegralSimpson := s;
end;
begin
writeln('Metodom Trapecij = ', IntegralTrap(a, b):0:5);
writeln('Metodom Simpsona = ', IntegralSimpson(a, b):0:5);
readln
end.