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