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 Thiersch

Biometria, Inventário e Manejo Florestal

Engenharia Florestal – UFSCar - Campus Sorocaba