% The computation of a0 leads to a real mess of surds.  It's not clear
% a further simplication exists.
% on rounded;  % Worse with this on!
% on fort;
p := -(a5*a9+a6*a8)/(2*(a6*a9-a7^2));
q := (p^2-a5*a8/(a6*a9-a7^2))^(1/2);
g0 := p+q;
g1 := p-q;
r0 := (a5+a6*g0)^2;
r1 := (a5+a6*g1)^2;
let a0 = a5*g0^2*a7^/r0+a8+a6*g0^3*a7^2/r0;
end;
b0 := a0-2a7^2*g0^2/sqrt r0+a9*g0;
a1 := a5*g1*a7^2/r1+a8+a6*g1^3*a7^2/r1;
b1 := a1-2a7^2*g1^2/sqrt r1+a9*g1;
end;
