home *** CD-ROM | disk | FTP | other *** search
/ Oakland CPM Archive / oakcpm.iso / sigm / vol292 / scilib.lbr / TRAP3.PAS < prev    next >
Encoding:
Pascal/Delphi Source File  |  1986-12-17  |  512 b   |  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. -r