home *** CD-ROM | disk | FTP | other *** search
/ CP/M / CPM_CDROM.iso / mbug / mbug003.arc / JRTPAS-1.LBR / JRT24B.PAS < prev    next >
Pascal/Delphi Source File  |  1979-12-31  |  512b  |  22 lines
  1. {Program 4.5
  2. compute the cosine using the expansion:
  3. cos( x ) = 1 - x**2/( 2 * 1 ) + x**4/( 4 * 3 * 2 * 1 ) - ...}
  4.  
  5. Program cosine;        { * ( input, output ) is implicit in Pascal/Z * }
  6.  
  7. const   eps = 1.0E-14;
  8. var    x, sx, s, t : real;
  9.     i, k, n : integer;
  10.  
  11. begin    
  12.     for i := 1 to n    do
  13.     begin
  14.            read( x ); t := 1; k := 0; s := 1; sx := sqr( x );
  15.         while abs( t ) > eps * abs( s ) do
  16.         begin
  17.             k := k + 2; t := -t * sx/( k * ( k-1 ) ); s := s + t
  18.         end;
  19.         writeln( x, s, k div 2 )
  20.     end
  21. end.
  22.