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


procedure kvadr(n:integer,a,it:real;var ik:real;var p:integer);

begin
k:=1/2;p:=0;
repeat
k:=k*2;n:=n*k;
for i:=1 to n do
begin
  x[i]:=a+(i-1/2)*h;ih:=ih+h*(exp(1/x)/(x*x));
   end;
   for i:=1 to 2*n do
   x[i]:=a+(i-1/2)*h;
   ihd:=ihd+(h/2)(exp(1/x)/(x*x));
   ik:=ihd+1/3(ihd-ih);e:=abs(it-ik);p:=p+1;h:=h/2;
   until e<0.0001;
   end;
   procedure sim(n:integer,a,it:real;var is:real;
   var p:int);
   var i:integer;x,f:array[1..70]of real,k,n,ih:real
   h:=(b-a)/(2*n);
   repeat k:=k*2;n:=n*k;
   for i:=1 to 2*n+1  do
   begin
   x[i]:=a+h*i;f[i]:=(exp(1/x)/(x*x));end;
   ih:=h/3*(f[1]+f[2n+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 4n+1 do
   begin
   x[i]:=a+h*i;f[i]:=(exp(1/x)/(x*x));
   ihd:=h/3*(f[1]+f[4n+1];
   i:=0;repeat i:=i+2;ihd:=ihd+(4/3)*h*f[i];until i=2n;
   i:=1;repeat i:=i+2;ihd:=ihd+(2/3)*f[i]*h;until i=2n-1 ;
   is:=ihd+(ihd-ih)/15
   e:=abs(it-is);
   until e<0.0001;
   end;
   begin{osn progr}
   writeln('vvedite n,a,b');
   readln(n,a,b);h:=(b-a)/n;it:=exp(1/a)-exp(1/b);



Соседние файлы в папке Alex