Advertisement
Guest User

Untitled

a guest
Mar 19th, 2018
78
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Pascal 0.99 KB | None | 0 0
  1. program Int_17;
  2.  
  3. {$APPTYPE CONSOLE}
  4.  
  5. //uses
  6. //  SysUtils;
  7. label
  8.   m1;
  9. var
  10.   a,b,i,j,n,iter:integer;
  11.   h,eps,s1,s2,x,intN:real;
  12.   f:array[1..100000] of real;
  13.  
  14. begin
  15.  
  16.   s1 := 0; s2 := 0; iter := 1;
  17.   a := 0; b := 2; // predel
  18.   n := 4;
  19.   eps := 0.001; // tochnost'
  20.   h := ( b - a ) / n;
  21.   x := a + h / 2;
  22.  
  23.   i := 1;
  24.   while i <= n do
  25.     begin
  26.       f[i] := exp( sqrt( 1 + x + 0.2*sqr(x) ) );
  27.       s1 := s2 + abs( f[i] ) * h;
  28.       x := x + h;
  29.       Inc(i);
  30.     end;
  31.  
  32.   m1:   n := n * 2;
  33.   h := ( b - a ) / n;
  34.   x := a + h * 2;
  35.  
  36.   j := 1;
  37.   while j <= n do
  38.     begin
  39.       f[j] := exp( sqrt( 1 + x + 0.2*sqr(x) ) );
  40.       s2 := s2 + abs( f[j] ) *h;
  41.       x := x + h;
  42.       Inc(j);
  43.     end;
  44.  
  45.   if abs( s2 - s1 ) < eps then
  46.     begin
  47.       intN := s2;
  48.       writeln('Integral: ', intN:10:3 );
  49.     end
  50.   else
  51.     begin
  52.       s1 := s2;
  53.       s2 := 0;
  54.       iter := iter + 1;
  55.       goto m1;
  56.     end;
  57.  
  58.   writeln( 'Summa iteraciy ', iter:10 );
  59.   readln;
  60. end.
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement