Если нет задачи мониторить давления, то надо: 1. убрать давление из параметров p вообще 2. в строке 153 заменить на P=get_partial_pressures!(P_total, u)
function RnC4H8(P::AbstractVector, EA::AbstractVector, EP::AbstractVector, ER::AbstractVector, kA::AbstractVector, kP::AbstractVector, kR::AbstractVector, T) -2*R1(P, EA, EP, ER, kA, kP, kR, T) - R4(P, EA, EP, ER, kA, kP, kR, T) - R6(P, EA, EP, ER, kA, kP, kR, T) end вот функция из 110 строчки