1 /SQUARE ROOT SUBROUTINE OS8 FORTRAN II LIBRARY LIBRARY
11 /COPYRIGHT (C) 1974,1977 BY DIGITAL EQUIPMENT CORPORATION
22 /THE INFORMATION IN THIS DOCUMENT IS SUBJECT TO CHANGE WITHOUT NOTICE
23 /AND SHOULD NOT BE CONSTRUED AS A COMMITMENT BY DIGITAL EQUIPMENT
24 /CORPORATION. DIGITAL EQUIPMENT CORPORATION ASSUMES NO RESPONSIBILITY
25 /FOR ANY ERRORS THAT MAY APPEAR IN THIS MANUAL.
27 /THE SOFTWARE DESCRIBED IN THIS DOCUMENT IS FURNISHED TO THE PURCHASER
28 /UNDER A LICENSE FOR USE ON A SINGLE COMPUTER SYSTEM AND CAN BE COPIED
29 /(WITH INCLUSION OF DIGITAL'S COPYRIGHT NOTICE) ONLY FOR USE IN SUCH
30 /SYSTEM, EXCEPT AS MAY OTHERWISE BE PROVIDED IN WRITING BY DIGITAL.
32 /DIGITAL EQUIPMENT CORPORATION ASSUMES NO RESPONSIBILITY FOR THE USE
33 /OR RELIABILITY OF ITS SOFTWARE ON EQUIPMENT THAT IS NOT SUPPLIED BY
46 / VERSION NUMBER IS AVAILABLE AT ENTRY POINTS
50 / IF X0 IS AN APPROXIMATION FOR Y**(1/2)
51 / THEN (X0+(Y/X0))/2 IS A BETTER APPROXIMATION
59 SQER, 6321 /"SQRT" ERROR FROM LOC XXXXX
71 CALL 1,STO /Y=ARGUMENT
75 TAD ACH /IF Y NEGATIVE THEN ERROR
82 POS, SZA /IF Y=0 THEN ROOT=0
85 NONZ, RAR CLL /FORM INITIAL APPROXIMATION
100 TAD (-3 /DO 3 ITERATIONS
108 CALL 1,FAD /(Y/X0)+X0
110 TAD ACH /((Y/X0)+X0)/2