Я пытаюсь решить это уравнение ниже (см. Изображение) для «t» в R, но застрял в сценариях. Я пробовал некоторые функции, такие как uni root, оптимизировать и оптимизировать, но все еще мог получить что-то. Любая помощь будет приветствоваться. введите описание изображения здесь
Here are the equation and my lines of codes.
#soil parameters
Ks= 17000;
Beta= 1.3;
Eta= 4.5;
Theta_s= 0.33;
Theta_r= 0.03;
Psi_Water= -80;
Delta=Eta/Beta
#crop parameters
Tp= 5.5;
EC50= 7.5505;
p= 3;
Psi_Root=-6000;
b= 10;
Yr0= 0;
#Water Salinity (EC)
EC <- seq(0.5,6, 1)
#Irrigigation water
I <-seq(0.5, 12, by =.1)
#the function created for the equation
fct <- function (t){-t + (min(Tp,(Psi_Root-(Psi_Water/((I-t)/Ks)^(1/Eta)))(I-t)b)/
(1+((ECI(Theta_r+(Theta_s-Theta_r)((I-t)/Ks)^(1/Delta)))/(EC50(I-t)*Theta_s))^p))}
for (i in 1:length(EC)) {
for (j in 1:length(I)){
ans <- optim (fct( Tp, EC50, p, Psi_Root, b, Ks, Eta, Theta_s, Theta_r, Psi_Water, Delta,
I[j], EC[i],t), c(0.1, 13))
}
}