DSolve[T'[t] = -k*T[t]^n] with T[0]=T0