Puoi usare un puntatore per passare il risultato

codice:
void radq(double x, double *y)
{
      double e = 0.01; // Qui assegna tu una precisione richiesta costante ... ad es. 0.01
      
      *y=1;
      do
      {
          *y = ((*y) + x/(*y))/2;
      }while(fabs(((*y)*(*y)-x)/x)>e);
}
oppure puoi dichiarare la variabile y "globale" prima del main ...