Решить систему нелинейных уравнений

HAlex спросил: 11 июня 2018 в 10:59 в: r

Я пытаюсь решить следующую систему из четырех уравнений. Я попытался использовать пакет "rootSolve", но мне кажется, что я не могу найти решение таким образом.

Код, который я использую, следующий:

model <- function(x) {
F1 <- sqrt(x[1]^2 + x[3]^2) -1
F2 <- sqrt(x[2]^2 + x[4]^2) -1
F3 <- x[1]*x[2] + x[3]*x[4]
F4 <- -0.58*x[2] - 0.19*x[3]
c(F1 = F1, F2 = F2, F3 = F3, F4 = F4)
}
(ss <- multiroot(f = model, start = c(0,0,0,0)))

Но это дает мне следующую ошибку:

Warning messages:
1: In stode(y, times, func, parms = parms, ...) :
error during factorisation of matrix (dgefa);         singular matrix
2: In stode(y, times, func, parms = parms, ...) : steady-state not reached

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

Спасибо вам!

0 ответов