/* @(#)msqrt.c 4.1 %G% */
msqrt(a
,b
,r
) MINT
*a
,*b
,*r
;
if(a
->len
<0) fatal("msqrt: neg arg");
if(a
->len
%2==1) x
.len
=(1+a
->len
)/2;
x
.val
=xalloc(x
.len
,"msqrt");
for(j
=0;j
<x
.len
;x
.val
[j
++]=0);
if(a
->len
%2==1) x
.val
[x
.len
-1]=0400;
sdiv(&y
,2,&y
,(short *)&j
);