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