<p dir="ltr">Olá pessoal obrigado pelas dicas, problema solucionado.</p>
<p dir="ltr">Abraço </p>
<div class="gmail_quote">Em 02/10/2015 17:58, "Jônatan" <<a href="mailto:jdtatsch@gmail.com">jdtatsch@gmail.com</a>> escreveu:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Wagner, veja a adaptação de seu código comentada abaixo.</div><div><br></div><div><div># função</div><div>Lmax_fun <- function(par){</div><div>  L <- par</div><div>  Ho <- 12</div><div>  Hvar <- 0.09061577</div><div>  So <- -0.03</div><div>  k <- 4.182488e-06</div><div>  ((Ho*Hvar/((-So-(k*L^1.75))*(1-Hvar)))-L)^2</div><div>}</div><div><br></div><div>L <- 1:300</div><div>Ho <- 12; Hvar <- 0.09061577; So <- -0.03; k <- 4.182488e-06</div><div>DEN_L <- ((-So-(k*L^1.75))*(1-Hvar))*L</div><div>plot(L,DEN_L,type = "l", ylim = c(0,max(DEN_L)*1.10))</div><div>abline(h=Ho*Hvar)</div><div><br></div><div># otimização unidimensional - optimize()</div><div># resultado depende do intervalo escolhido</div><div># pelo gráfico da função podemos ver os intervalos das raízes</div><div>optimize(f = Lmax_fun, interval = c(min(L), max(L)))</div><div>optimize(f = Lmax_fun, interval = c(1, 200))</div><div>optimize(f = Lmax_fun, interval = c(100, 300))</div><div><br></div><div>library(pracma)</div><div># findmins() divide o intervalo de busca n vezes (n = 100 por default)</div><div>mins <- findmins(f = Lmax_fun, min(L), max(L))</div><div>points(mins, rep(Ho*Hvar, 3), pch = 20)</div><div><br></div><div>#Condição para valer a escolha da função</div><div>mins[abs(So)/(k*mins^1.75) >= 1.75+1] </div></div></div><div class="gmail_extra"><br><div class="gmail_quote">2015-10-02 14:48 GMT-03:00 Wagner Wolff <span dir="ltr"><<a href="mailto:wwolff@usp.br" target="_blank">wwolff@usp.br</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div><div dir="ltr"><div><div dir="ltr"><div><font size="2" style="font-family:arial">###############################################################</font><br><span style="font-family:arial;font-size:small">##  </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Jônatan Dupont Tatsch</span></div><div><span style="font-family:arial;font-size:small">##  </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Professor do Departamento de Física</span><font face="arial" style="font-size:small"><br></font></div><div><span style="font-family:arial;font-size:small">## </span><span style="font-family:arial;font-size:small"> </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.9090909957886px;line-height:12.9971590042114px;text-align:center;background-color:rgb(239,239,239)">Coordenador Substituto do Programa de Pós-Graduação em Meteorologia (PPGMET)</span><br></div><div><span style="font-size:small;font-family:arial">##  </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Centro de Ciências Exatas e Naturais (CCNE)</span><br></div><div><span style="font-size:small;font-family:arial">##  </span><span style="color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;font-size:10.909090995788574px;line-height:12.997159004211426px;text-align:center;background-color:rgb(239,239,239)">Universidade Federal de Santa Maria - UFSM</span></div><div><span style="font-family:arial;font-size:small">##  </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Faixa de Camobi, Prédio 13 - Campus UFSM - </span><span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Santa Maria, RS, Brasil - 97105-900</span><br></div><div style="font-family:arial;font-size:small">##  <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px">Telefone: <a href="tel:%2B55%2855%2933012083" value="+555533012083" target="_blank">+55(55)33012083</a></span><br></div><div style="font-family:arial;font-size:small">##  <span style="font-size:11px;text-align:center;background-color:rgb(239,239,239);color:rgb(150,150,150);font-family:Verdana,Arial,Helvetica,sans-serif;line-height:13px"><a href="http://www.ufsm.br/meteorologia" target="_blank">www.ufsm.br/meteorologia</a></span></div><div style="font-family:arial;font-size:small">###############################################################</div></div></div></div></div>
</div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div>