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