##############################################################################
#                          FÍSICA COMPUTACIONAL II                           #
#                                   por                                      #
#                         Francisco  Carlos  Lavarda                         #
##############################################################################

Tutorial de MAXIMA - Equações Diferenciais Ordinárias 
                     (Problemas de Valores de Contorno)
=======================================================

** Digite a equação diferencial para a qual se quer a solução:
'diff(y,x,2)-4*'diff(y,x)+4*y=exp(2*x)/x^2; 

sendo que na equação digitada:
-  'diff(y,x,2) é a derivada de segunda ordem de y em relação a x;
-  'diff(y,x) é derivada de primeira ordem de y em relação a x.

** Use a função ODE2 para resolver a equação que foi digitada acima:
ode2(%,y,x);

Isto vai gerar uma solução geral, bem mais ampla que solução específica
deste problema, cuja solução analítica é y=-exp(2x).ln(x) .
No sistema norte-americano, a função ln escreve-se como log, fato necessário
para entender o que vem a seguir.

** Use a função BC2 para ajustar a solução a um problema de valores de
contorno em uma EDO de ordem 2:
bc2(%,x=1,y=0,x=%e,y=-%e**(2*%e));

O símbolo % no MAXIMA quer dizer que você se refere à última resposta dada
por ele.

** DICA: ao invés de usar o símbolo % nos comandos, que indica a resposta
imediatamente anterior dada pelo programa, você pode usar o símbolo que
aparece entre parênteses na frente da resposta que o programa dá, que é um %
mais uma letra e um número. A vantagem deste último método é o de poder se
referir a uma resposta dada pelo MAXIMA em qualquer instante  posterior em
que se precisar dela.