[R-br] Função Solve

[Cláudio Roberto Thiersch]

Walmes, Boa noite!
    Muito Obrigado pela dica sobre o CAS e o wxMaxima.

    Depois de sua dica e tentando ajudar o Alessandro  busquei mais informações sobre CAS e encontrei o pacote rSymPy.  Infelizmente para o modelo (taxa de variação do incremento médio) apresentado pelo Alessandro não consegui entender o erro que o sistema retornou, por isso, instalei o software SymPy e executei o seguinte 
código no R. 
cat('import sympy\nfrom sympy import *\n\ndef main():', file='ioc.py', append = F,'\n');cat(paste('  b0, b1, b2, b3, b4, b5, id1, ab1, s, x = symbols(',shQuote('b0 b1 b2 b3 b4 b5 id1 ab1 s x','csh'),')',sep=''), file='ioc.py', append = T,'\n');cat('  ioc = solve(exp(b0 + b1/s + b2/x + b3 * (id1/x) * log(ab1) + b4 * (1 - id1/x) + b5 * (1 - id1/x) * s) * (b4 * (id1/x**2) - (b3 * (id1/x**2) * log(ab1) + b2/x**2) + b5 * (id1/x**2) * s)/x - exp(b0 + b1/s + b2/x + b3 * (id1/x) * log(ab1) + b4 * (1 - id1/x) + b5 * (1 - id1/x) * s)/x**2,x)', file='ioc.py', append = T,'\n');cat(paste('\n  print\n  pprint(ioc)\n  print\n\nif __name__ ==', shQuote('__main__'),':\n  main()\n'), file='ioc.py', append = T);system('python ioc.py');


Resultado: [-b2 - b3*id1*log(ab1) + b4*id1 + b5*id1*s]
      Caso alguém queira reproduzir, apenas tomem cuidado com os espaços e tabulações porque o "python" é um pouco exigente.  
At.,
Cláudio ThierschBiometria,Inventário e Manejo Florestal

Engenharia Florestal– UFSCar - Campus Sorocaba



  
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20150501/238e4342/attachment.html>