Построение нелинейной функции в R

В настоящее время я пытаюсь построить следующую функцию в R.

(F1/F2^m * (R1 - F1)^((1 - s)/s)) - (m * (((R1 - F1))^(1/s) + ((R2 - F2))^ (1/с)) / (F1^m + F2^m)) = 0

Где я хочу построить график F2 как функцию F1. То есть F2 по оси Y и F1 по оси X. Значения R1,R2 являются фиксированными константами, а s,m также фиксированы.

Поскольку функция нелинейна, мне предложили использовать числовую аппроксимацию, чтобы найти значения F2 в ранее определенной области F1.

Я разработал следующий код, определяющий сначала константы, затем функцию и функцию поиска корня следующим образом:

rm(list=ls())

library(ggplot2)

R1_values <- c(100)  # Example values for R1
R2_values <- c(100)  # Example values for R2
s_values <- c(1)  # Example values for s
m_values <- c(1)  # Example values for m

# Define the original function to solve for F2
F2_function <- function(F1, F2, R1, R2, s, m) {
  (F1 / F2^m * (R1 - F1)^((1 - s) / s)) - (m * (((R1 - F1))^(1 / s) + ((R2 - F2))^(1 / s)) / (F1^m + F2^m))
}


F1_values <- seq(0.01, 100, by = 0.1)
# Solve for F2
F2_values <- sapply(F1_values, function(F1) {
  tryCatch({
    # 
    uniroot(F2_function, c(0.01, 100), F1 = F1, R1 = R1, R2 = R2, s = s, m = m, tol = 1e-8)$root
  }, error = function(e) NA)
})

data <- data.frame(F1 = F1_values, F2 = F2_values)

ggplot(data, aes(x = F1, y = F2)) +
  geom_line() +
  labs(title = "F2 as a Function of F1",
       x = "F1",
       y = "F2") +
  theme_minimal()

Однако все возвращает NA или чрезвычайно большие значения (случайно зная, когда и как). Даже если функция четко определена в другом программном обеспечении (Geogebra), изображение прилагается.

Следовательно, я хотел бы попросить помощи в том, как построить график (и восстановить значения F2 для каждой части области F1).

Уравнение в математическом плане выглядит так:

П.Д.: Я думаю, что у F2 есть два результата на каждую F1, но я тоже не могу разобраться в этом.