void SolutionSlau degreeMethod double lMax double lMax2 std vector dou

void SolutionSlau::degreeMethod(double &lMax, double &lMax2, std::vector<double> &xMax, std::vector<double> &xMax2)
{
    int size = A.size();

	std::vector<double> yCur(size, 1);
	std::vector<double> xCur(size, 1);
	std::vector<double> yNext(size);
	std::vector<double> xNext(size);
	std::vector<double> lambda(size, 0);
	std::vector<double> lambdaNext(size);


	for (int k = 0; k < 100000; ++k)
    {
        for (int i = 0; i < size; ++i)
        {
            yNext[i] = 0;
            for (int j = 0; j < size; j++){
                yNext[i] += A[i][j] * xCur[j];
            }
			lambdaNext[i] = yNext[i] / xCur[i];
        }
		double n = 0;
		for (int i = 0; i < size; ++i)
			n = std::max(n, fabs(yNext[i]));
		
		for (int i = 0; i < size; ++i)
			xNext[i] = yNext[i] / n;

		std::swap(xNext, xCur);
		std::swap(yNext, yCur);
		std::swap(lambda, lambdaNext);

		double t = 0;
		for (int i = 0; i < size; ++i)
			t = std::max(t, fabs(lambda[i] - lambdaNext[i]));
		if (t < eps)
			break;
    }
	xMax = xCur;
	lMax = 0;
	for (int i = 0; i < size; ++i)
		lMax += lambda[i];
	lMax /= size;

}