Sì ... guardavo solamente la prima colonna ...

Secondo me, il codice dovrebbe essere

codice:
void radice_newton(int &n,int v1[],double E,double rq2[])
{     
    int i;

    for(i=0;i<n;i++)
    {
        rq2[i]=1;
        do
        {
            rq2[i]=(rq2[i] + v1[i] / rq2[i])/2;
        } while(fabs((rq2[i] * rq2[i] - v1[i]) / v1[i]) >E);
    }
}