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
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
