home *** CD-ROM | disk | FTP | other *** search
/ ftp.barnyard.co.uk / 2015.02.ftp.barnyard.co.uk.tar / ftp.barnyard.co.uk / cpm / walnut-creek-CDROM / MBUG / MBUG003.ARC / JRTPAS-1.LBR / JRT25.PAS < prev    next >
Pascal/Delphi Source File  |  1979-12-31  |  1KB  |  34 lines
  1. {Program 4.6
  2. NOTE: Have patience with this one. I thought it was broke
  3.     until I looked at the code carefully. It takes
  4.     awhile. JRT took ten minutes,(I didn't use a stop
  5.     watch) while Pascal/Z took 5 minutes. But JRT
  6.     had more digit accuracy. 
  7. computer 1 - 1/2 + 1/3 - ... + 1/9999 - 1/10000, 4 ways:
  8.     1) left to right, in succession
  9.     2) left to right, all pos and neg terms, then subtract
  10.     3) right to left in succession
  11.     4) right to left, all pos and neg terms, then subtract}
  12.  
  13. Program summing;    
  14.  
  15. var    s1, s2p, s2n, s3, s4p, s4n, lrp, lrn, rlp, rln : real;
  16.     i : integer;
  17.  
  18. begin
  19.     s1 := 0; s2p := 0; s2n := 0; s3 := 0; s4p := 0; s4n := 0;
  20.     for i := 1 to 5000 do
  21.     begin
  22.         lrp := 1/( 2 * i - 1 );        {pos terms, left to right}
  23.         lrn := 1/( 2 * i );        {neg terms, left to right}
  24.         rlp := 1/( 10001 - 2 * i );    {pos terms, right to left}
  25.         rln := 1/( 10002 - 2 * i );    {neg terms, right to left}
  26.         s1 := s1 + lrp - lrn;
  27.         s2p := s2p + lrp; s2n := s2n + lrn;
  28.         s3 := s3 + rlp - rln;
  29.         s4p := s4p + rlp; s4n := s4n + rln;
  30.     end;
  31.     writeln( s1, s2p - s2n );
  32.     writeln( s3, s4p - s4n )
  33. end.
  34.