home *** CD-ROM | disk | FTP | other *** search
/ Simtel MSDOS 1992 September / Simtel20_Sept92.cdr / msdos / turbopas / pas_sci.arc / TRAP3.PAS < prev    next >
Pascal/Delphi Source File  |  1985-09-06  |  504b  |  30 lines
  1. program trap3;          { -> 270 }
  2. { integration by the trapezoidal rule }
  3.  
  4. const   tol             = 1.0E-6;
  5. var     sum,upper,lower : real;
  6.  
  7. function fx(x: real): real;
  8. { find f(x)=1/x }
  9. { watch out for x=0 ! }
  10. begin
  11.   fx:=1.0/x
  12. end;
  13.  
  14. function dfx(x: real): real;
  15. begin
  16.   dfx:=-1.0/sqr(x)
  17. end;
  18.  
  19. {$I trapez.lib}
  20.  
  21. begin           { main program }
  22.   ClrScr;
  23.   lower:=1.0;
  24.   upper:=9.0;
  25.   writeln;
  26.   trapez(lower,upper,tol,sum);
  27.   writeln;
  28.   writeln(chr(7),'area=',sum)
  29. end.
  30.