Advertisement
Guest User

Untitled

a guest
Nov 15th, 2018
111
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Pascal 0.75 KB | None | 0 0
  1. program abc;
  2.  
  3.  
  4. {$mode objfpc}{$H+}
  5. uses Classes, SysUtils,math;
  6.  
  7.  
  8. const k=2000;
  9. var a,b,h,s1,s2,s:real;
  10.  i,n,ntest:integer;
  11.  y:array[1..k+1] of real;
  12.  
  13. function f(x:real):real;
  14. begin
  15.  result:=power(cos(3*x+2),6)/sqrt(1+x*x);
  16. end;
  17.  
  18. begin
  19.  
  20.  writeln('a = ');
  21.  readln(a);
  22.  writeln('b = ');
  23.  readln(b);
  24.  ntest:=2;
  25.  while ntest>0 do
  26.  begin
  27.   Writeln('Enter nmber of integerals, less 2000');
  28.   writeln('n = ');
  29.   readln(n);
  30.   ntest:=n mod 2;
  31.  end;
  32.  h:=(b-a)/n;
  33.  y[1]:=f(a);
  34.  y[n+1]:=f(b);
  35.  for i:=2 to n do
  36.  begin
  37.   y[i]:=f(a+h*(i-1));
  38.  end;
  39.  s1:=0;
  40.  s2:=0;
  41.  for i:=2 to n do
  42.  begin
  43.   if (i mod 2)=0 then
  44.    s2:=s2+y[i]
  45.   else
  46.    s1:=s1+y[i];
  47.   end;
  48.  s:=(y[1]+y[n+1]+s1*2+s2*4)*h/3;
  49.  writeln('Integral = ',s);
  50.  readln;
  51. end.
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement