Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
25
Добавлен:
09.12.2013
Размер:
1.12 Кб
Скачать
program lab2;
uses crt;

function simp(a,b:real):real;

function f(x:real):real;
begin
f:=x*x*x*cos(3*ln(x));
end;

begin
simp:=(b-a)/6*(f(a)+4*f((b-a)/2)+f(b));
end;

var J1,J2,J,hi,ei,ai,bi,ri,R,h,e0:real;
i:integer;

const a=0.5;
b=2;
e=0.001;

begin
clrscr;
writeln('ЪДДДВДДДДДДДДДДДДДВДДДДДДДДДДДДДВДДДДДДДДДДДДДї');
writeln('і ь і  Ј hi і €­вҐЈа « Ji і Џ®ЈаҐи­. ei і');
writeln('ГДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДґ');
i:=1;
ai:=a;
bi:=b;
h:=0;
e0:=0;
R:=0;
repeat
hi:=bi-ai;
ei:=e/(b-a)*hi+e0;
J1:=simp(ai,bi);
J2:=simp(ai,ai+hi/2)+simp(ai+hi/2,bi);
ri:=abs(J2-J1)/15;
if ri<=ei
then
begin
writeln('і',i:2,' ','і',h:13:10,'і',J1:13:10,'і',ri:13:10,'і');
i:=i+1;
J:=J+J1;
ai:=bi;
bi:=b;
e0:=ei-ri;
R:=R+ei;
h:=h+hi
end
else
begin
bi:=ai+hi/2;
e0:=0
end;
until ai=b;
writeln('ГДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДЕДДДДДДДДДДДДДґ');
writeln('ібг¬і',h:13:10,'і',J:13:10,'і',R:13:10,'і');
writeln('АДДДБДДДДДДДДДДДДДБДДДДДДДДДДДДДБДДДДДДДДДДДДДЩ');
readkey;
end.
Соседние файлы в папке laba22