а у меня же есть alglib::rmatrixinverse

#include <iostream>
#include <cstdio>
#include <time.h>

long long int binomial_C(int n, int k) {
    long long int result = 1;
    if(k > n-k) {
        k=n-k;
    }
    for(int i=0; i<k; i++) {
        result = result*(n-i)/(i+1);
    }
    return result;
}

int main()
{
    clock_t tStart = clock();
    int x,y;
    long long int binom_1, binom_2;
    for(x=1; x<100; x++) for(y=1; y<100; y++) {
        binom_1 = binomial_C(x-1, y);
        binom_2 = binomial_C(x, y-1);
        if(binom_1 == 10 && binom_2 == 5*x/2) {
            std::printf("Success\nx: %d, y: %d\n", x, y);
        }
    }
    std::printf("Time taken: %.5f seconds\n", (double)(clock() - tStart)/CLOCKS_PER_SEC);
    return 0;
}

Success
x: 6, y: 3
Time taken: 0.00324 seconds

Без компилятора слишком сложно. :<

У кого есть проц на 6ГГц

я сделал. только обратную матрицу для якобиана я нашел исключениями гаусса. не хочу прикручивать alglib из-за одной функции, потому что алглибские матрицы ниче  не умеют (складываться, умножаться)

пихайте в меня системы уравнений, насмехайтесь надо мной

А какого размера потянешь?))

на тебе систему

Первое уравнение - гипербола?

я чо-т таких гипербол не видал

там иррац-ть есть

а если её убрать, то 4-я степень

Ой да. Оно не похоже на 1/х

Хотя... Похоже...

-2,0 вроде подходит

Судя по всему, это еще и единственное решение)

как заставить ньютона не подходить к уже найденным корням? для 1-мерных ф-й надо разделить на бином (x-x_0). а с многомерными как?

init == {{-9.38216}, {-5.1402}}
sol == {{-2}, {1.44031e-08}}
nt.F(sol) == {{2.22045e-16}, {2.07448e-16}}
-----------
roots == {{{-2}, {1.44031e-08}}}

больше не могу найти корней

Сейчас сам поищу)

я делю исходную ф-ю на произведение (x-r_i) по всем координатам всех найденных корней. но такая редуцированная ф-я иногда выдает мусор, далекий от реальных корней.