Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- program Int_17;
- {$APPTYPE CONSOLE}
- //uses
- // SysUtils;
- label
- m1;
- var
- a,b,i,j,n,iter:integer;
- h,eps,s1,s2,x,intN:real;
- f:array[1..100000] of real;
- begin
- s1 := 0; s2 := 0; iter := 1;
- a := 0; b := 2; // predel
- n := 4;
- eps := 0.001; // tochnost'
- h := ( b - a ) / n;
- x := a + h / 2;
- i := 1;
- while i <= n do
- begin
- f[i] := exp( sqrt( 1 + x + 0.2*sqr(x) ) );
- s1 := s2 + abs( f[i] ) * h;
- x := x + h;
- Inc(i);
- end;
- m1: n := n * 2;
- h := ( b - a ) / n;
- x := a + h * 2;
- j := 1;
- while j <= n do
- begin
- f[j] := exp( sqrt( 1 + x + 0.2*sqr(x) ) );
- s2 := s2 + abs( f[j] ) *h;
- x := x + h;
- Inc(j);
- end;
- if abs( s2 - s1 ) < eps then
- begin
- intN := s2;
- writeln('Integral: ', intN:10:3 );
- end
- else
- begin
- s1 := s2;
- s2 := 0;
- iter := iter + 1;
- goto m1;
- end;
- writeln( 'Summa iteraciy ', iter:10 );
- readln;
- end.
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement