<div dir="ltr"><div><div>Olá pessoal!<br></div>Estou precisando de uma ajuda, preciso encontrar uma solução para uma variável que terei que usar um método iterativo para estimá-la, entretanto, existe duas soluções (raízes) ou mais para esta variável e preciso da segunda solução pois é a que tem explicação física, sendo assim, a função optmim me retorna a primeira, gostaria que retornasse todas as soluções possíveis, alguém saberia como posso fazer isso? Ou indicar outra função. Segue o código e o gráfico mostrando as duas soluções, onde a linha horizontal intersepta o perfil.<br><br></div>Abraço<br><div><div><div><br>Ho <- 12<br>Hvar <- 0.09061577<br>k <- 4.182488e-06<br><br>Lmax_fun <- function(par,Ho,Hvar,So,k){<br>    L <- par<br>    ((Ho*Hvar/((-So-(k*L^1.75))*(1-Hvar)))-L)^2<br>}<br>So <- -0.03<br>(Lmax <- optim(par=100,fn=Lmax_fun,Ho=Ho,Hvar=Hvar,k=k,So=So,method ="BFGS"))<br><br>abs(So)/(k*Lmax^1.75) >= 1.75+1 #Condição para valer a escolha da função<br><br>L <- 1:300<br>DEN_L <- ((-So-(k*L^1.75))*(1-Hvar))*L<br><br>plot(L,DEN_L,type = "l",ylim = c(0,max(DEN_L)*1.10))<br>abline(h=Ho*Hvar)<br>points(x=Lmax,Ho*Hvar)<br>Lmax2 <- 129.3757<br>points(x=Lmax2,Ho*Hvar)<br><br><br><img src="cid:ii_if9y11ex0_15029aa56df7790f" height="400" width="442"><br>​<br></div></div></div></div>