code 34340;
real procedure COMABS(XR,XI);value XR,XI;real XR,XI;
begin XR:= ABS(XR); XI:= ABS(XI);
COMABS:= if XI > XR then SQRT((XR/XI)**2+1)*XI
else if XI= 0 then XR else SQRT((XI/XR)**2+1)*XR
end COMABS;
        eop