Алгоритм

}dir;

}

// длина шага в заданном направлении

double length (vector<double> dir, vector<double> grad, vector<vector<double> > A) {a = 0;<double> B (dir.sz, 0);(int i = 0; i < grad.sz; i++) {+= grad[i] * dir[i];

}(int i = 0; i < dir.sz; i++) {(int j = 0; j < A[i].sz; j++) {[i] += dir[j] * A[j] [i];

}

}a / inner_prod (B, dir);

}

// метод сопряженных градиентовConGrad (vector<vector<double> > A, vector<double> B, vector<double> &X) {<double> gP = B, d (B.sz, 0), gN = B, c;s;(conv(gN) > EPS) {= grad (A, B, X);(conv(gN) < EPS) break;= direction (d, gP, gN);= length (d, gN, A);= nextApp (X, s, d);

gP = gN;

}

}