Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

la2

.pas
Скачиваний:
21
Добавлен:
09.12.2013
Размер:
2 Кб
Скачать
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('bez bankirovania');
 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;
  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;
 writeln('bankirovanie');
 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.
Соседние файлы в предмете Вычислительные методы