sábado, 13 de janeiro de 2018

Resolvendo uma EDO de 2a. ordem com "Maxima"


Sobre o Maxima (não traduzido, mas o inglês é bem fácil de entender):

Maxima is a computer algebra system, implemented in Lisp.

Maxima is derived from the Macsyma system, developed at MIT in the years 1968 through 1982 as part of Project MAC. MIT turned over a copy of the Macsyma source code to the Department of Energy in 1982; that version is now known as DOE Macsyma. A copy of DOE Macsyma was maintained by Professor William F. Schelter of the University of Texas from 1982 until his death in 2001. In 1998, Schelter obtained permission from the Department of Energy to release the DOE Macsyma source code under the GNU Public License, and in 2000 he initiated the Maxima project at SourceForge to maintain and develop DOE Macsyma, now called Maxima.

Vejamos um exemplo simples:

f: 'diff(y,t,2) + 5*'diff(y,t) + 9*y = 5+3*cos(4*t);
g: ode2(f,y,t);
h: ic2(g,t=0,y=0,'diff(y,t)=2);

A primeira linha define uma EDO de 2a. ordem e armazena em f. A solução geral dessa equação é armazenada em g. Em h temos a solução particular para as condições iniciais y(0) = 0 e y'(0) = 2. Quando executado (Crtl+R) gera:


Que é a solução desejada. Um segundo exemplo:



*** Link para download do Maxima aqui.

Nenhum comentário:

Postar um comentário