
function [n1,n2]=g_sqrt(n)
  tab=factor_list(n)
  f1 = floor(tab(:,2)/2)
  f2 = tab(:,2)-2*f1
  n1 = prod(tab(:,1).^f1)
  n2 = prod(tab(:,1).^f2)
  if n<>n1^2*n2
    error('bug in mysqrt')
  end
endfunction

