Нахождение минимума квадратичной функции

#include <iostream>
#include <math.h>
#include <conio.h>
#include <iomanip>
#include <fstream>
#include <stdlib.h>
 
using namespace std;
//исходная функция
double f(double a11, double a12, double a13, double a14, double a22, double a23, double a24, double a33, double a34, double x1, double x2, double x3)
{
    return a11*x1*x1+2*a12*x1*x2+2*a13*x1*x3+a22*x2*x2+2*a23*x2*x3+a33*x3*x3+2*a14*x1+2*a24*x2+2*a34*x3;
}
//производная по x1
double ffx1(double a11, double a12, double a13, double a14, double a22, double a23, double a24, double a33, double a34, double x1, double x2, double x3)
{
    return 2*a11*x1+2*a12*x2+2*a13*x3+a14; 
}
//производная по x2
double ffx2(double a11, double a12, double a13, double a14, double a22, double a23, double a24, double a33, double a34, double x1, double x2, double x3)
{
    return 2*a12*x2+2*a22*x2+a23*x2+2*a24;
}
//производная по x3
double ffx3(double a11, double a12, double a13, double a14, double a22, double a23, double a24, double a33, double a34, double x1, double x2, double x3)
{
    return 2*a13*x1+a23*x2+2*a33*x3+2*a34;
}
int main()
{
    setlocale( LC_ALL,"Russian" );
    double x1;
    double x2;
    double x3;
    double a11=1.2475;
    double a12=0.847;
    double a13=0.4332;
    double a22=1.0135;
    double a23=-0.0556;
    double a33=0.739;
    double a14=-1.6227;
    double a24=-1.2001;
    double a34=-0.6323;
    double eps;
    cout << "Выберите метод решения квадратичной функции:\n" << endl;
    cout << "1) Градиентный метод\n" << endl;
    cout << "2) Метод Ньютона\n" << endl;
    cout << "3) Метод сопряжённых градиентов\n" << endl;
    int k;
    cin>>k;
    if ((k!=1)&&(k!=2)&&(k!=3))
    cout << "Неверный выбор!" << endl;
    else 
    {
        cout << "Введите точность eps" << endl;
        cin >> eps;
    }
    system("PAUSE");
}