Найти значение переменной с помощью функции Solver

Question

Я пытаюсь настроить функцию "Солвер", чтобы найти нормальную глубину канала (yn). Параметры приведены в коде ниже, где я могу оценить одну сторону уравнения. Все остальные параметры являются функцией yn. Мне нужно найти yn, который решает функцию A * (R ^ (2/3) = nQSo.

So=0.001                                       
n=0.013                                        
Q=30  
B=10 

nQSo=(n*Q)/(So^(1/2))

A=B*yn
P=B+2*yn
R=A/P

A*(R^(2/3)=nQSo

Answer 1

Вы можете взглянуть на optimize

So=0.001                                       
n=0.013                                        
Q=30  
B=10 

nQSo=(n*Q)/(So^(1/2))

error = function(yn,nQSo){
  A=B*yn
  P=B+2*yn
  R=A/P

  return(abs(A*(R^(2/3))-nQSo))
}

optimize(error,interval = c(0,2),nQSo = nQSo)

результат, как вы видите, yn = 1.239066