3 ! iflag = 0 (RK-adaptativo), 1 (RK-fixo), 2 (Forward Euler), 3 (RK-adaptive-alternative)
2   1.	! n = numero de equacoes, tau0 = tempo inicial
1000.  0.1 1.d-4  ! tend = tempo final, dtau = passo externo de integracao para mostrar resultados, teta0 = valor inicial
1.d-8  ! tol = tolerancia para RK-adaptativo alternativo
1. ! alfa = parametro para RK-adaptativo alternativo
4. ! p = ordem do metodo usado para integrar
0.001 ! ht2 = passo interno inicial usado para integrar as EDOs
5. ! dtab = intervalo de tempo para construir a tabela	numerica