[unix-history] / tests / t08.p

program summing(output);
var
	s1, s2p, s2n, s3, s4p, s4n,
	lrp, lrn, rlp, rln: real;
	i: integer;
begin
	s1 := 0;
	s2p := 0;
	s2n := 0;
	s3 := 0;
	s4p := 0;
	s4n := 0;
	for i := 1 to 500 do
	begin
		lrp := 1/(2*i-1);
		lrn := 1/(2*i);
		rlp := 1/(10001-2*i);
		rln := 1/(10002-2*i);
		s1 := s1 + lrp - lrn;
		s2p := s2p + lrp;
		s2n := s2n + lrn;
		s3 := s3 + rlp - rln;
		s4p := s4p + rlp;
		s4n := s4n + rln;
	end;
	writeln(s1, s2p-s2n);
	writeln(s3, s4p-s4n);
end.
Commit	Line	Data
fde98a2d BJ	1	program summing(output);
	2	var
	3	s1, s2p, s2n, s3, s4p, s4n,
	4	lrp, lrn, rlp, rln: real;
	5	i: integer;
	6	begin
	7	s1 := 0;
	8	s2p := 0;
	9	s2n := 0;
	10	s3 := 0;
	11	s4p := 0;
	12	s4n := 0;
	13	for i := 1 to 500 do
	14	begin
	15	lrp := 1/(2*i-1);
	16	lrn := 1/(2*i);
	17	rlp := 1/(10001-2*i);
	18	rln := 1/(10002-2*i);
	19	s1 := s1 + lrp - lrn;
	20	s2p := s2p + lrp;
	21	s2n := s2n + lrn;
	22	s3 := s3 + rlp - rln;
	23	s4p := s4p + rlp;
	24	s4n := s4n + rln;
	25	end;
	26	writeln(s1, s2p-s2n);
	27	writeln(s3, s4p-s4n);
	28	end.