#include ".cminpack.h"
#include <math.h>

int equation_sytem( void* ptr , int n , const double* estimations , double* errors, double *fjac, int ldfjac, int iflag )
{
            if ( iflag == 1 ) // вычислените правых частей системы упавнений
            {
                errors[0] = sqrt(estimations[0]) + estimations[1]*estimations[1] - estimations[2] - 3.0;
                errors[1] = (estimations[0] + estimations[1])/estimations[2] - 2.0;
                errors[2] = estimations[0] * estimations[1] * estimations[2] - 24.0;
            }
            if ( iflag == 2 ) // вычисление матрицы Якоби
            {
              //вычисление элементов матрицы Якоби
              fjac[0] = 0.5 * pow ( estimations[0] , -0.5 );   fjac[3] = 2*estimations[1];                   fjac[6] = -1;
                 fjac[1] = 1.0/estimations[2];                    fjac[4] = 1.0/estimations[2];                 fjac[7] = -(estimations[0] + estimations[1])/(estimations[2]*estimations[2]);
                 fjac[2] =  estimations[2] * estimations[1];    fjac[5] =  estimations[0] * estimations[0];  fjac[8] =  estimations[1] * estimations[2];
            }
    return 1;
}

int main()
   {
             double vecX[3] ; // вектор переменных
         double vecE[3] ; // вектор правых частей (ошибок пиближения)
         double matJ[3*3]; // матрца Якоби

         int n = 3; // порядок системы
         int ldfjac = n;
         int lwa = 330; // в доке написано что должно быть 10*(n*(3*n+13))/2;
         double wa[330];
         double    tol = 1e-12;  // точность вычислений
             // начальные условия
         vecX[0] = 1.0; vecX[1] = 1.0; vecX[2] = 1.0;
            // решение системы
         hybrj1( equation_sytem , (void*)NULL , 3 , vecX , vecE , matJ, ldfjac ,tol , wa , lwa );
            
             // корни уравнения в vecX

            return 0;
}