comment hulpprogramma R 764, MLP 090762 / 1929;
begin	real V, u; array f[1: 2];
	real procedure ZERO (x,a,b,fx,e);
	value a,b; real x,a,b,fx; array e;
	begin	real c,fa,fb,fc,re,ae;
		re:= e[1]; ae:= e[2];
		x:= a; fa:= fx; x:= b; fb:= fx; goto entry;
	begin:  1f abs (a - b) <= fa then a:= b + sign (c - b) * fa;
		if sign (a - x) = sign (b - a) then x:= a;
		a:= b; fa:= fb; b:= x; fb:= fx;
		if sign (fc) = sign (fb) then
	entry:  begin c:= a; fc:= fa end;
		if abs (fb) > abs (fc) then
		begin a:= b; fa:= fb; b:= c; fb:= fc; c:= a; fc:= fa end;
		a:= (a * fb - b * fa) / (fb - fa);
		x:= (c + b) / 2; fa:= abs (b * re) + ae;
		if abs (x - b) > fa then goto begin;
		ZERO:= x
	end ZERO;
	f[1]:= 0; f[2]:= &apos;-6;
	for u:= .15891, .10805, .02351, .14542, .09780, .02014 do
	begin	NLCR; FIXT (0, 5, u);
		FIXT (0, 6, ZERO (V, 0, 1, .0008404 * V ** 5 - .021009 * V ** 4 + .23951 * V ** 3 - .23111 * V ** 2 + 1.2219 * V - u, f))
	end
end