for i = 1:n
      for j = i+1:n
         println(r)
         sinA = r[j][i] / sqrt(r[i][i] * r[i][i] + r[j][i] * r[j][i])
         cosA = (sinA * r[i][i] / r[j][i])
         println(sinA)
         println(cosA)
         for k = i:n
            rj = r[j][k]
            r[j][k] = - cosA * r[j][k] + sinA * r[i][k]
            r[i][k] = cosA * r[i][k] + sinA * rj
         end
         for k = 1:n
            qj = q[j][k]
            q[j][k] = - cosA * q[j][k] + sinA * q[i][k]
            q[i][k] = cosA * q[i][k] + sinA * qj
         end
      end
   end