2 ! iflag = 0 (RK-adaptativo), 1 (RK-fixo), 2 (Forward Euler)
2   1.	! n = numero de equacoes, tau0 = tempo inicial
3.  0.1 ! tend = tempo final, dtau = passo de integracao
2. 2. ! tp0(i) = valores iniciais (i=1,n)