Скачиваний:
329
Добавлен:
04.03.2014
Размер:
2 Кб
Скачать
program int;
var a,b,ik,is,H,IT:real;
i,p,N:integer ;


procedure kvadr(n:INTEGER;a,H,it:real;var ik:real;var p:integer);
VAR K:INTEGER;X:ARRAY[1..3435] OF REAL;
IH,IHD,E:REAL;
begin
k:=1;p:=0;
repeat
;n:=n*k;
for i:=1 to n do
begin
  x[i]:=a+(i-1/2)*h;ih:=ih+h*(exp(1/x[I])/(x[I]*x[I]));
   end;WRITELN(P);
   for i:=1 to 2*n do
   x[i]:=a+(i-1/2)*h;
   ihd:=ihd+(h/2)*(exp(1/x[I])/(x[I]*x[I]));
   ik:=ihd+(1/3)*(ihd-ih);WRITELN(IK);READLN;e:=abs(it-ik);p:=p+1;h:=h/2;K:=K*2;



   until e<0.0001;
   end;
   procedure sim(n:integer;a,B,it:real;var is:real;
   var p:intEGER);
   var i,K:integer;x,f:array[1..3435] of real;ih,E,IHD,H:real
  ; BEGIN
  h:=(b-a)/(2*n);
   repeat k:=1;n:=n*k;
   for i:=1 to 2*n+1  do
   begin
   x[i]:=a+h*i;f[i]:=(exp(1/x[I])/(x[I]*x[I]));end;
   ih:=h/3*(f[1]+f[2*n+1])
   ;i:=0;repeat i:=i+2;ih:=ih+4/3*h*f[i];until i=2*n;
   i:=1; repeat i:=i+2;ih:=ih+2/3*f[i];until i=2*n-1;
 h:=h/2;  for i:=1 to 4*n+1 do
   begin
   x[i]:=a+h*i;f[i]:=(exp(1/x[I])/(x[I]*x[I]));
   ihd:=h/3*(f[1]+f[4*n+1]);END;
   i:=0;repeat i:=i+2;ihd:=ihd+(4/3)*h*f[i];until i=2*n;
   i:=1;repeat i:=i+2;ihd:=ihd+(2/3)*f[i]*h;until i=2*n-1 ;
   is:=ihd+(ihd-ih)/15
  ; e:=abs(it-is);P:=P+1;K:=K*2;
   until e<0.0001;
   end;
   begin{osn progr}
   writeln('vvedite n,a,b');P:=0;
   readln(n,a,b);h:=(b-a)/n;it:=exp(1/a)-exp(1/b);WRITELN('IT--',IT);
    KVADR(N,A,H,IT,IK,P);
      WRITELN('METOD KVADRATOV-----  ',IK,'PROGONOK----  ', P);
      P:=0 ;SIM(N,A,B,IT,IS,P)  ;
    WRITELN(' METOD SIMPSONA---  ',IS,'PROGONOK--- ',P);
    READLN;READLN
    END.
Соседние файлы в папке Alex